M>-A050  213  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CALIF  F/6  9/2 

CMS  HISTOORAM*  DENSITY  ESTIMATION  AND  PROBABILITY  PLOTTING  ROUT— ETC (U) 
dec  77  e I DANIKAS 

UNCLASSIFIED  HI- 


1 of3 

AD 

A0502I3 

a 

p 

' — — 

• 

k 

f 

1 ( 

1 

1 

■ 

• 

i 

■ 

G 

A 

i 

1 

J 

NAVAL 


POSTGRADUATE  SCHOOL 

Monterey,  California . 


i/ih  .t 


/he  SIS, 


fn 


FEO 


■’7^ 


CMS  Histogram,  Density  Estimation  and 
Probability  Plotting  Routines,  with  an 
I Application  to  the  Analysis  of  the  Output 


U 


of  a Simulation  of  a Correlated  Queue*  ^ 


leorgios  loanniy Danikas 


by 


-.N 


s7d 


'7  / 


Thesis  Advisor: 


P.  A.  W.  Lewis 


D 


Approved  for  public  release;  distribution  unlimited. 


SCCumTy  CLASSiriCATION  or  this  ^ACC  (Whm  Omm 


REPORT  DOCUMENTATION  PAGE  bbfome^cmple^o'form 


2.  OOVT  ACCCSSION  NOJ  I-  nCCI^)CHT*t  CATALOO  NUMMCH 


4.  TITLE  fold  SuA'KUJ  / *•  TYVC  OE  dlPOIIT  A PCmOO  COVEMCO 

CMS  Histograun,  Density  Estimation  and  / Master’s  Thesis; 

Probability  Plotting  Routines,  with  an  December  1977 

Application  to  the  Analysis  of  the  Output  ••  ^«"^om4iMa oa«.  aepoat  numee* 

r\f  a ^ 1 mil  1 ^-1  nn  0‘F  Por 7*0 1 a OiiAne 


7.  AgTHOUf*; 

Georgios  loannis  Danikas 


t.  AEEEOMMINC  OHOANIZATION  NAME  ANO  AODNESE 

Naval  Postgraduate  School 
Monterey,  California  93940 


II.  CONTHOLLINO  OPEICE  NAME  ANO  AOONESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


U.  MONITOHINC  AOENCV  name  a tiOOntSSal  dlltaiwnt  Inm  Cantnlllng  Olllcm)  IS.  SECUNITT  CLASS,  (al  OHa  riAen) 

Naval  Postgraduate  School  Unclassified 

' California  939  40  isa.  oeclassieication'  oowngnaoihc 

SCHEDULE 

IS.  OISTKIEuTION  statement  (el  title  Kepatt) 

Approved  for  public  release;  distribution  unlimited. 


12.  MCPOfRT  date 

/ 

December  1977 

1).  NUMBCR  OP  ^AGCS 

232 

I 17.  OlSTfllSUTION  STATCMCNT  (of  «A«rr«cf  mtf4  In  8tPek  20,  H dlUtmt  hpm  Hppoet) 


fltTsUliWLCMeNTAnY  MOTES 


1 19.  KEY  WOflOS  (Contifium  «n  it  npcmmmpfr  tdPPtiir  Ar  Moek  nymkPf) 


Histogram 

Plotting  subroutines 
Exponential  Moving  Average  (EMA) 


ao.  j^^^WACT  (Cpnanum  gw  rpvprmm  044m  tt  mcpppptt  PTd  idmtiifr  kf  ktopk  mmtkpr) 

The  object  of  this  thesis  has  been  twofold.  The  first 
object  was  to  develop  FORTRAN  versions  of  several  existing  APL 
programs  which  were  designed  to  analyze  univariate  data.  In 
particular  the  programs  were  designed  to  test  for  exponential- 
ity  and  normality  of  the  data  and,  by  sectioning  or  jacknifing, 
obtain  estimates  of  sampling  variances  of  sample  moments.  The-^ 


00 


EOHM 
JAN  71 


EDITION  OE  I NOV  EE  lE  OEtOLBTI 
S/M  0l02>9l4*«401  I 


tECUNITY  CLAEEIEICATION  OF  TNIE  NAEI  | 


1 


• 

second  object  of  the  thesis  was  to  use  these  programs  in  a 
simulation  study  of  first-come  first-served  queues  in  which 
the  service  times  and  the  inter-arrival  times  are  exponen- 
tially distributed  but  dependent.  The  dependence  is  intro- 
duced by  using  the  mixed  moving  average  autoregresive  struc- 
ture (EARMA  (p,  q) ) for  .exponential  sequence^introduced  by 
Lewis  and  co-workers.  ^Four  models  of  correlated  queues  are 
introduced,  giving  autocorrelated  and  cross-correlated 
service  and  arrival  times  in  various  degrees.  The  simulation 
study  gives  a quantitative  idea  of  the  effect  of  correlation 
on  the  mean  waiting  time  and  the  distrib’  i of  the  waiting 
time. 


DD  Fornix  14.3 
I Jan  >3 

S/N  0102-014-6601 


SICuAl'f''  CLAMI^lCATlOM  0»  THIS  Dim 


2 


Approved  for  public  release;  distribution  unlimited 


'f-. 


CMS  HISTOCaAM,  DENSITX  ESTIMATION  AND  PROEABILITY  PLOTTING 
BOOTINJS,  ilTH  AN  APPLICATION  TO  THE  ANALYSIS  OF  THE  COTPOT 
OF  A SIMOLATION  OF  A COSHELATEC  QOEOE 


by 


Georgios  loannis  Danikas 
Major  Hellenic  Aray 
3.S.,  Military  Academy  of  Greece,  1963 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degrees  of 


MASTEB  OF  SCIENCE  IN  COMPaTEfi  SCIENCE 
HASTES  OF  SCIENCE  IN  OPEfiATIONS  EESEAHCH 


from  the 

NAVAL  POSTGHADOATE  SCHOOL 
December  1977 


n D c 


1 ' 


fn  : - 1978  11 


Author: 


Approved  by: 


HU:.. 


-u  U LULH 


D 


yf  ^ A f ^ 


or 


GERALD  G.  GN^WN 


Second  Feader 


nd  Policy  Sciences 


1 


&ESTRACT 


The  coject  of  this  thesis  has  oeen  tuofold.  Ihe 
first  ctject  nas  to  develop  FOBTHAN  versions  of 
several  existing  APL  piograois  which  were  designed  tc 
analyze  univariate  data.  In  particular  the  programs 
were  designed  to  test  for  exponentiality  and  normality 
of  the  data  and,  by  secticnirg  or  jacknifing,  obtain 
estimates  of  sampling  variances  of  sample  moments. 
The  second  object  of  the  thesis  was  tc  use  these 
programs  in  a simulation  study  of  firsr-ccie 
first-served  queues  in  which  the  service  times  and  the 
inter-arrival  times  are  exponentially  distributed  tut 
dependent.  The  dependence  is  introduced  by  using  the 
mixed  moving  average  autoregresive  structure 
( EAai3A(p,g)  ) for  exponential  sequences  introduced  by 
Lewis  and  co-workers.  Pour  models  cf  correlated 
queues  are  introduced,  giving  autoccrr elated  and 
cross-correlated  service  and  arrival  times  in  various 
degrees.  The  simulation  study  gives  a quantitative 
idea  cf  the  effect  cf  correlation  on  the  mean  waiting 
time  and  the  distribution  of  the  waiting  time. 
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I . INTRODUCTIOM 


The  ocject  of  this  thesis  has  been  twofold. 

The  first  object  was  to  develop  FCEIBAN  v 
several  existing  API  programs  which  were  designed 
univariate  data,  in  particular  to  test  fcr  expo 
and  normalitj  of  the  data  and,  by  sectioning  or 
obtain  estimates  of  variances  of  sample  moments. 


ersicns  of 
to  analyze 
nentiality 
jackrifing 


The  second  object  was  to  use  these  programs 
simulation  study  of  first-come  first-served  gueues  in 
the  service  times  and  the  inter-arrival  times 
exponentially  distributed  but  dependent.  The  depends 
introduced  by  using  the  mixed  moving  average  autoregr 
structure  ( EAailA(p,g)  ) for  exponential  seg 
introduced  by  Gaver  and  Lewis  (1977),  Lawrance  and 
(1977),  Jacobs  and  Lewis  (1977)  and  Lawrence  and 
(1978).  Four  schemes  are  introduced  which 
autocorr elated  and  cross-ccrrelated  service  and  a 
times  in  various  degrees. 


in  a 
which 
are 
nee  is 
essive 
uences 
Lewis 
Lewis 
give 
rr ival 


Since  nc  analytical  properties  of  these  queues  can  be 
derived,  their  properties  are  studied  by  simulation. 
Because  the  EAEUA  gueues  are  neither  regenerative  nor 
Markovian  there  are  several  problems  in  the  simulation  study 
of  knowing  when  one  can  assume  that  a steady  state  has  been 
reached  in  the  simulation  and  of  obtaining  confidence 
intervals  for  estimated  parameters.  This  prcblem  is  handled 
by  generating  a number  of  independent  sample  paths  and 
cemparing  statistics  obtained  from  the  replications  at 
several  pcints  along  the  sample  paths.  Only  the  waiting 
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not  the  state 


1 


tiae 

space 


prccess  (W  ) is  studied  in  this  thesis 
prccess  M (t)  . 


9 


Unfortunately  the  simulaticn  study  c£  the  EAEHi-type 
gueue  uas  hanpered  by  the  size  and  speed  c£  the 
laM  Systen  360/67  systea  at  the  Naval  Postgraduate  School. 
Conseguently  the  results  were  not  as  extensive  as  had  keen 
hoped  fcr.  Fairly  detailed  results  are  given  for  the  case 
where  the  service  tine  process  is  autoregressive;  one  case 
where  the  traffic  intensity  is  0.99  and  the  correlation 
parameter  0.98,  was  intensively  investigated.  This  gives 
seme  idea  of  the  length  of  the  transient  in  the  gueue  and  of 
the  inflation  of  the  mean  waiting  time  because  of  the 
correlation  in  the  service  process.  Some  investigation  of 
the  case  where  the  service  times  and  inter-arrival  times  are 
cross-correlated  are  given. 

In  section  II  we  describe  the  basic  histograc  and 
statistics  computation  package  'HISIGS/F'  and  apply  it  to 
two  sets  of  telephone  error  data.  In  section  III  the 
plotting  method  and  the  plotting  Subroutines  'NCFNPl', 
'FXPLI*  are  described  and  we  test  them  with  various 
generated  data.  In  section  IV  we  describe  the  assessment  of 
variability  Subroutines  'S2CTN',  'JACK'  and  we  apply  them  tc 
telephone  data  set  1.  In  section  V the  Subreutine  'LIST'  is 
described  and  we  apply  it  to  two  sets  of  telephone  error 
data.  In  section  VI  we  describe  and  analyze  the  EAE!JA(p,g) 
model  as  well  the  program  which  we  developed  in  order  to 
simulate  it.  Computer  program  listings  are  provided  after 
section  VI. 


j 
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Here  the  Subtoutine  'HISTGS/HISTFS ' is  presented. 
' HISTGS/HISIfS ' is  used  for  estimating  the  profcahility 
density  function  from  a given  set  of  data  and  at  the  same 
time  computing  some  basic  statistics.  Basically  this 
Subroutine  is  the  library  FORTRAN-Subroutine  HISTF/G,  which 
was  created  at  N.P.S.  by  D.  H.  Robinson.  It  has  been 
fflcdified  by  adding  the  new  entry  point  'LIMIIS'. 

The  new  entry  point  may  be  used  for  wild  data  or  for 
data  having  a mixed  distribution.  Therefore  it  is  a useful 
tool  for  that  user  involved  with  this  kind  of  data.  Simply, 
the  user  has  to  define  those  values  in  the  range  of  data 
(sectioning  the  data  into  several  disjoint  sections)  that 
he/she  believes  are  useful  cut  points,  for  each  section  of 
data  the  user  may  have  the  probability  density  function  as 
well  as  the  basic  statistics  ( as  HISIF/G  does  ) . Alsc  by 
sectioning  the  data  over  its  range,  the  user  may  have  a 
histogram  of  the  random  variable  when  it  is  conditioned  to 
be  between  given  limits. 

The  number  of  sections  that  the  entry  point  'LIMITS'  can 
accept  is  as  many  as  5C. 


A complete  description  of  how  ' aiSTGS/HlSTFS ' operates 
is  given  in  the  subroutine.  Cn  the  other  hand,  a summary  of 
the  subroutine  is  given  by  typing  on  the  terminal  the 


J 
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ccaaand  C£5CEIE£  HISTGS  under  CMS  (Camtridce  monitor  System, 
fief.  4).  Typing  'DE5CBIBE  HISTGS*  the  following  response 
will  be  given  on  the  terainal: 


[ SOBfiCOIIBE  HISTGS/HISTFS 

t 

[ 

’ HISTGS/iilSTPj ' gives  you  a histogram  of  a set  of  data 
along  with  the  estimate  of  a set  of  basic  statistics.  Also 
it  gives  you  the  following  options: 

1.  Displaying  a smoothed  empirical  probability  density 
function  ever  the  histogram, 

2.  Sectioning  the  data  into  several  disjoint  sections 

(no  more  than  SO)  and  then  having  a histogram  and  the 

estimate  of  basic  statistics  for  each  section, 

3.  Scaling  the  histogram, 

4.  Displaying  just  the  histogram  (with  cr  withoct  the 

density  function)  and  no  statistics. 

The  esiimated  basic  statistics  are:  Hean , Hedian, 

Irimean,  Hidrange,  Geometric  Hean,  Harmonic  Kean,  Variance, 
Standard  Deviation,  Coefficient  of  Variation,  Mean 
Deviation,  fiange.  Midspread,  Third  and  Forth  Central  Moment, 
Skewness,  Kurtosis,  Beta),  Beta2,  Maximum,  Minimun,  and 
i^uantiles. 

CALLING  SEQUENCES 


CALL 

HISIGS 

( 

X,  N, 

NBAfi  ) 

CALL 

HISIf S 

( 

X,  N, 

NBAfi  ) 

CALI 

LIMITS 

( 

X,  N, 

XLIM, 

CALL 

NCSIAS 

CALL 

STATS 

CALL 

FIXS  ( 

SCALE  ) 

CALL 

NCFIXS 

CALL 

FfiCM AX 

( 

FSC  ) 

CALL 

NGFfiMX 
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ABGUHSNTS 

X Array  of  data  values 

M Muober  of  data  values 

M£AB  Number  of  bars  in  the  histogram 
Xim  Array  of  cut  points 
llflS  Number  of  cut  pcints 

SCALE  Vector  of  two  values  to  scale  the  data 
PSC  A real  variable  (between  0.  and  1.)  defining 
the  maximum  value  of  the  probability  axis. 

If  the  entry  point  NOSTAS  is  called  before  calling  any 
of  the  entry  points  HISTFS,  HISTGS,  or  LIHIIS  then  no 
statistics  are  printed,  otherwise  the  statistics  are  printed 
by  default. 

ay  calling  the  entry  point  FIXS  the  histogram  scale  may 
be  fixed  and  remains  set  unless  it  is  reset  ty  another  call 
to  FIXS  or  allowed  to  float  (the  default)  bj  calling  NCFIXS. 
The  reason  for  fixing  the  scales  is  so  that  when  comparing 
mere  than  one  batch  of  data,  we  can  get  comparably  scaled 
histograms. 

By  calling  PSOKAX,  the  maximum  value  of  the  probability 
axis  (Y)  can  be  set  to  PSC,  which  will  remain  ir  effect 
until  NOPHMX  is  called  (the  default). 

More  information  is  given  in  tne  subroutine. 

E.  USING  ' HISTGS/HISTFS*  WIIH  TELEPHONE  DATA  1 AND 
TELEPHONE  DATA  2 


Here  'HISIGS/HISTFS'  has  been  used  to  plot  the  set  of 
Telephone  Data  1 and  Telephone  Data  2.  The  algorithmic 
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procedure  for  both  sets  of  data  is  the  same.  That  is,  the 
entry  point  'HI5TFS*  is  called  first  and  then  the  entry 
point  'LIMITS’.  Pour  cut  points  have  been  used  (1,  2,  141, 

86000)  sectioning  the  data  into  three  disjoint  sections  as 
fellows: 

Section  1 : Consists  of  data  points  X , such  that 

<i> 

1 < X <2. 

<i> 

Section  2 : Includes  the  data  points  X , such  that 

2 < X < 141. 

Ci> 

Section  3 : Includes  the  data  points  X , such  that 

Ci) 

142  < X < 86000. 
tiJ 

No  print-out  for  section  1 is  given  because  of  the 
constant  value  of  the  data  points. 

Figures  1 through  6 give  the  histogram  and  the  basic 
statistics  for  both  sets  of  data  ( For  the  entire  set  and 
fer  sections  2 and  3 ) . 


Me  observe  that  both  of  data  sets  are  sc  unruly  that  the 
use  of  the  entry  pcint  'LIMITS'  of  'HISTGS'  is  suggested. 


Observing  Figures  1 and  4 we  may  conclude  that  both  data 
sets  appear  to  be  the  same,  having  seme  geometric  or 
exponential  distribution.  But  as  we  will  see  later  this  is 
net  true.  Thus  the  histogram,  in  that  case,  which  is 
obtained  by  using  just  the  entry  point  'HISTFS'  is  net 
helpful  and  leads  us  to  a wrong  conclusion. 

'LIMITS'  is  applied  now  to  give  us  a useful  answer  ; 

Observing  figures  2 and  5,  obtained  by  using  the  entry 
point  'LIMITS',  we  may  conclude  that  the  data  sets  dc  not 
have  the  same  distribution  ( at  least  in  that  interval  ) . 
In  addition  data  set  2 is  far  from  having  a geometric  cr  an 
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exponential  distribution,  as  it  can  be  seen  by  the  histogram 
of  the  entire  set  of  data.  This  is  a conclusion  which 
cannot  be  obtained  by  without  use  of  the  entry  point 
•LiaiTS' . 

Observing  also  Figures  3 and  6 , where  the  histogram  of 
the  3rd  section  from  each  data  set  is  presented,  we  may 
conclude  that  this  section  of  both  data  sets  may  have  an 
exponential  distribution,  but  not  the  same  one. 

The  mosT  important  facet  of  the  data  which  appears  is 
the  modes  at  about  124.  The  mode  is  much  more  apparect  in 
Telephone  Data  2.  In  fact  closer  investigation  reveals  that 
this  mode  was  due  to  leakage  of  dial  pulses;  this  leakage 
caused  bit  errors  at  multiples  of  124  bits- 

The  above  is  a quick  and  informal  analysis  of  the 
Telephone  Data  Sets  obtained  by  using  the  entry  point 
•LIMITS*  of  • HISTGS/HISTFS ' , and  the  conclusion  is  that  both 
of  data  sets  are  so  unruly  that  any  formal  characterization 
of  their  distribution  is  not  so  easily  obtained.  In 
addition  it  is  apparent  that  a mixed  distribution  governs 
both  sets  of  data. 

Having  quickly  analyzed  both  of  the  telephone  data  sets, 
the  use  of  entry  point  'LIMITS*  of  the  Subroutine 
• HISTGS/aiSIFS ' is  obvious. 
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Ill . PLOT II NG  SUBROUIINES 


A.  INTBOCUCIION 

Graphical  aethods  of  assessing  the  validity  cf  a 
probability  aodel,  and  of  estiaating  soae  basic  paraieters 
(especially  Iccation  and  scale  parameters)  have  been  widely 
applied  over  many  years.  A variety  of  different 
prescri ptions  have  been  advanced  for  the  plotting  positions, 
that  is,  for  the  set  of  values  at  which  the  ordered 
observations  in  a sample  should  be  plotted. 

Plotting  the  ordered  data  is  an  infcrmal  and  quick 
visual  method  for  getting  an  impression  of  how  tne  ! 

distribution  of  data  looks.  Depending  on  what  is  known 
about  the  distribution  of  the  data  there  are  many  \ 

probability  plotting  methods.  | 

j 

The  Plotting  Subroutines  presented  here,  'NOfiMFL'  and  j 

•EXPLI',  use  the  method  of  plotting  the  order  statistics  j 

( X ) versus  their  expected  value  ( EfX  ] ) called 
<iJ  ci) 

•SCOHES'.  This  method  does  not  require  knowledge  cf  the 

Cumulative  Distribution  Function  F (x)  (except  for  its 

assumed  continuity) . Therefore  it  may  be  very  attractive 

for  those  who  know  nothing  about  the  data  to  be  analyzed.  In 

addition  the  indication  of  the  graph  may  lead  the  analyst  to 

further  useful  results. 


f 
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Using  the  Plotting  subroutines  'NOEMPL'  and  'EXPLT'  one 
aay  analyze  the  resulting  graphs  obtained  tc  get  an  idea  of 
the  distribution  of  the  data. 


E.  DESCilPIICS  0?  THE  PLOTTING  HETHOD 


The  idea  behind  this  methcd,  as  above  stated,  is  tc  plot 
versus  E[ X . ]• 


It  is  used  where  data  X , X , X 

12  n 


arise  as 


independent  observations  of  a continuous  Bacdoo  Variable  X 


with  a distribution  function  which  is  believed  to  have  seme 


particular  form  F (x) 


Then  the  ordered  data 


X , X , X are  plotted  against  EfX  ]. 

Cl)  C2)  Ca)  ci) 


If  the  icdel  for  P(x)  used  to  get  E[X^  is  true  then 
the  resulting  plot  is  a linear  (regression)  relationship  and 

should  plot,  within  limits  of  sample  variations,  as  a 

straight  line.  Barnett  [2]  gives  the  value  of  E[X  ] as  a 

c i > 

linear  relation  of  the  plotting  position  That  is: 

E[X  1 = a + sx  , 

'■  C i)  ^ Ci> 


here  a,  and  s,  are  the  location  and  scale  parameters.  The 
Bean  and  Variance  of  X are  not  necessary.  As  we  see  the 
above  relation  uoiguely  defines  the  plotting  position  as  : 


X - B 

= EC ] . 

s 

Therefore,  if  the  aodel  is  correct  the  linearity  will 
probahly  be  ewident.  In  addition  if  the  resulting  plot  not 
being  linear  leads  us  to  reject  the  aodel,  we  may  cet  an 
indicaticn  of  the  type  of  departure  present  and  thus  an 
indication  of  what  action  tc  take  or  what  alternative  acdel 
to  use.  It  should  be  noted  that  the  linearity  of  the  graph 
suggests  informally  that  the  model  is  true,  but  rot  to 
accept  the  model.  Additional  formal  tests  must  be  done  in 
order  to  acoept  the  model. 

C.  SUBBCUTIKE  NCRMEl 


1 . Cescriptioa 

Subroutine  ’NCHMPL*,  is  used  to  test  a sample  of 
independent  observations  of  a continuous  Random  Variable  X 
having  a distribution  function  which  is  believed  to  have  the 
Ncraal  fora. 

'NORKEl'  sorts  the  data  into  increasing  order  to 
create  a set  of  Order  Statistics.  Then  the  Order  Statistics 
are  plotted  versus  either  Normal  Scores 

(if  n < 50,  n = number  of  data-points)  or  inverse  of  the 

- 1 

Standard  Normal  Cumulative  Distribution  Function  (F  (y)), 
(if  n > 50)  . 
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Evaluatioa  o£  Normal  Scores  is  tased  on  computing 

tb 

rhe  Expected  Value  of  the  i Order  Statistic  ( ) . 

Ibis  evaluation  follows; 

Let 

tb 

g(x  ) be  the  density  function  of  the  i Order 

Statistic  of  a sample  1 , 1 1 frcn  a population 

12  n 

having  the  Standard  Normal  Cistribution  f.  Then  g (x  ) is 

c i > 

given  by  : 

m i-1  n-i 

g (X  ) F (X  ) (1-F  (X  ))  f {X  ) 

<i>  <i-1)!  (n-i)!  X <i>  X <iJ  x ‘i) 

n*  i-1  n-i 

L p Q Z , 

(i-1)l  (n-i)! 

where : 

2 

-X  /2 
1 ii> 

2 = f (X  ) e , 

X <i>  1/2 

(2  n) 

2 

1 -u  /2 

P * P (X  ) e du  , 

X <i)  1 1/2 

<<^(2  n) 

Q = 1-E  . 

Now 

= ( X g (X)  dx 


The  value  of  expression  (1)  has  been  evaluated  and 
is  given  in  £19].  This  table  gives  the  values  of  Sormal 
Scores  frca  1 to  n/2  (for  n<50) . The  values  from  (n/2)+1  to 
n are  evaluated  using  the  symmetric  expressicn  : 


•NOEaEL*  uses  the  above  Normal  sccres  if  n < 50.  For 
n > 50  'NOaaEL*  computes  its  cun  Normal  Sccres  using  the 
program  function  'INVNBa',  which  gives  the  value  x cf  the 
inverse  Standard  Ncrmal  distributicn  (x  = F (y) , 

0 5 y i 1) . Justification  of  this  is  the  fact  that  if  the 
data  are  really  normally  distributed  then 

E[  X ] P(  i/ (n+ 1)  ).  Therefore  'NORaPL'  can  acconcdate 
< i ) 

any  sample  size. 


In  addition  to  the  plot  given  by  'NOBaEL',  the 
Hilk-Shapiro  test  for  normality  can  be  evaluated  by  the  same 
Subroutine.  A detailed  description  of  that  test  is  giver,  in 
[22].  However  a summary  of  that  test  is  presented  here  in 
order  to  help  us  to  describe  how  the  program  works.  The 
Hilk-Shapiro  test  is  based  on  the  statistic: 


where : 


k 

r a 

i=1  n-i+1 


(n-i  + 1) 


k=E/2,  n even 
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b = a (X  -X  ) ♦•  a (X  -X  } ♦ 

n tn>  t1>  n-1  <n-1>  <2) 


♦ a (X  -X  ) , 
k*2  <k+2}  <k> 


2 n 2 

S = ^ (X  - I) 

i=  1 1 


Ibe  coefficients  (a  ) are  defined  by: 

i 


k=  (n-1)  /2,  n odd 


a = S a u /C  f 
i j=1  j ij 


i — 1f2f.../n 


where : 


a . = E[X  ], 
3 ^3^ 


u = Cov[X  ,X  ], 
ij  <i)  <j) 


2 • -1  -1 
C = a y V a 


(That  is  the  'NORMAL  3CCHE3') 


i#j  - 1»2/.../n 


where: 


a - (a  ^a  /...a  ) 
12  n 


V = (u  ) , That  is  the  nXn  Covariance  Matrix, 
if  j 


There  are  some  approximations  associated  with  the 
evaluation  of  these  coefficients.  However  fcr  n < 5C  the 

(a  ) are  tabulated  in  the  above-mentioned  paper  ard  the 
i 

program  uses  them  by  FCBTHAN  DATA  initialization  stateients. 


’T 


Fee  n>50  the  (a  ) are  computed  by  the  program  using 
i 

the  follouing  approxioaxion  method: 


i i 

where  C is  given  by  the  following  Least-Square  equation  as  a 
function  of  n: 

2 

C = -2.722  ♦ 4.083n 


Also,  for  n=1 


2 _ + .5^ 

~ 1/2 _ 

X r(n/2  + 1) 

ilcte  that  a = a 

i n-i+1 


Some  properties  for  the  M Statistic  can  he  given 
here,  which  have  been  taken  from  the  above  mentioned  paper: 

(a)  . M is  scale  and  origin  invariant. 

(b)  . 'n  has  a distribution  which  depends  only  on  the 
sample  size  n,  for  samples  from  a Normal  Distribution. 

2 

(c)  . li  is  statistically  independent  of  S and  of  5 for 


samples  from  a Normal  Distribution. 


r 2r  2r 


(d) . 

E[M 

] = E[b  ]/ECS  ] 

for  any 

r . 

(6)  . 

The 

maximum  value  of 

W is  1 . 

(t)  . 

The 

sinimum  value  of 

2 

w is  na 

1 

/ (n-1)  . 

(g) . 

The 

half  and  first  moments  of 

H are  given  by: 

,.V^  H^^n/2  -.5) 

E[  W ] = 

1/2 

C y[n/2)  2 
2 2 


BCW] 


R (R  +1) 
2 


c (n-1) 


35 


where: 

2 • -t 

B - ■ 7 n 

(h)  . For  n=3  the  density  of  H is  given  as: 

3 -1/2  -1/2 

f^(w)  = (1-«)  w , .75  < w < 1. 


2.  Prog laa Structure 

'N0R2SPL'  is  a FOETSAN-callabie  Subroutine  with  each 
call  returning  (optionally)  either  a plot  fcx  a given  set  of 
data  or  a value  for  the  W Statistic  or  both. 

Ihe  prograa  is  divided  into  three  parts: 

The  first  part  is  the  ccctrol  prcgran  of  'NOEMEL'. 
In  that  part  the  'Noraal  Scores'  for  n = 2(1)25  and 
n = 26(2)50  are  stored  along  with  the  Sufcrcutine  ty  DATA 
stateaents;  on  the  other  hand,  the  Noraal  Scores  for  r > 50 
are  computed  ty  calling  the  program's  function  'INVNBM'.  In 
this  part  the  user  specifies  just  a plot  or  just  the  value 
of  U or  both.  The  Library-Subroutine  'PXSCET'  is  used  by 
the  program. 

The  second  part  of  the  program  ccnsisxs  of  the 
Subroutine  'PLOT'  which  accepts  any  data  set  to  be  tested 
for  Normalitj.  'PLOT'  itself  scales  the  data  according  to 
its  range  and  plots  the  scaled  data  along  110  equal  spaced 
positions  of  the  X-axis.  No  plot  is  given  if  the  data  has  a 
constant  value. 

The  third  part  of  the  program  is  the  'WILK' 

Subroutine.  In  this  Subroutine  the  coefficients  (a  ) 

i 

(for  n < 50)  of  the  N Statistic  are  stored  using  FCSTEAN 


36 


DATA  initialization  statements.  For  n>50  the  coefficients 

(a  ) are  comfuted  by  xhe  Subroutine  itself.  This  Subroutine 
i 

also  calls  the  'INVNEfl*  Function  to  get  the  ’Normal  Sccres' 

to  be  used  fcr  the  ccmfutaticn  cf  the  (a  ) 's  when  n > 50. 

i 

A complete  description  of  how  'NOBBEL’  operates  is 
given  in  the  Subroutine.  However  a summary  is  obtained  by 
typing  on  the  terminal  DESCRIBE  NOBHPL . When  the  user  types 
the  command  DESCRIBE  NCHHPL  under  the  CMS  environment  the 
following  response  is  printed  on  the  terminal: 


SOBBCUTINE 


NCRMPL 


test  Statistic 


'NOEMPl'  taJces  a set  of  data,  sorts  it  into 
ascending  order  and  uses  the  created  Order  Statistics  for; 

1.  Plotting  X versus  either  'Normal  Scores'  (if 

(i) 

n < 50,  n = Sample  size)  or  Inverse  cf  the  unit  Kcrmal 
Cumulative  Distribution  (that  is: 

-1 

X versus  Normal  scores  or  F { i/(n  + 1)  ) ,) 

(i) 

to  see  if  there  is  a linear  fit. 

2.  Computing  the  value  of  the  W 
(Silk-Shapixo  test  for  normality) . 

It  is  called  by: 

CALL  NCEMPL  ( X,  SCORES,  N,  K ) 

ARGUMENTS 

X Is  the  array  containing  the  data 

SCORES  Is  a worh  array  of  dimension  N 
N Is  the  number  of  data  values 

K User's  option 

for:  K = 1 A plot  only  , 
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K = 2 The  H-value  is  given  , 

K = 3 Both  of  the  above  are  given. 

Mere  informaticn  is  given  in  the  subroutine. 

3-  Interpreting  The  Output 

If  data  are  really  Normally  distributed  then  a 
straight  line  is  expected  from  the  output.  This  is  an 
indication  that  gives  the  user  a first  feeling  of  the 
distribution  of  data.  But  if  the  plot  is  not  linear  then 
the  user  has  to  reject  the  normality  assumption.  In 
addition  the  shape  of  the  plot  may  lead  the  user  to  make  an 
alternative  decision  of  what  model  ha  may  use.  For  example: 

(a)  . Shape  of  the  form  of  Figure  7a  suggests  that 
data  are  skewed  to  one  side  (specifically  positive  values) . 
In  that  case  an  Exponential  model  could  be  an  appropriate 
one.  The  histogram  of  the  data  is  given  in  Figure  7b. 

(b)  . Shape  of  the  form  of  Figure  8a  can  be 
interpreted  as  data  not  as  dispersed  as  the  Normal 
Distribution.  This  set  of  data  may  have  a Uniform  or  a 
Triangular  Distribution,  and  could,  for  instance,  arise  if 
measuring  tolerances  of  components  hand-picked  to  lie  within 
tolerances.  The  associated  histogram  is  given  in  Figure  8b. 

(c)  . Shape  of  the  form  of  Figure  9a  indicates  much 
more  dispersed  data  than  Normal,  but  the  data  is  symmetric, 
probably  ceming  from  a Double  Exponential  or  Cauchy 
Distribution.  The  associated  histogram  is  given  in  Figure 
9b. 

It  should  be  noted  here  that  a straight  line  does 
not  always  arise  from  Normal  data.  Figure  10a  for  example 
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I 


I 

i 

1 


gives  an  indication  that  data  are  Normal  (Straight  line)  . 
But  these  data  have  been  generated  by  computer  from  a 
Symmetric  Triangular  Distribution  in  the  range  (0  to  2).  On 
the  other  hand  figure  11a  plots  a data  set  cf  size  50  again 
from  a Triangular  (0  tc  2)  distribution.  But  in  this  case 
ue  observe  that  the  graph  is  not  a straight  line.  These 
results  demonstrate  the  informality  of  that  Probability 
Plotting  Method  and  in  these  cases  the  user  has  to  continue 
testing  the  data  by  other  formal  tests.  Ne vertheless,  since 
this  method  is  primarily  an  initial  screening  device,  ease 
cf  application  is  important.  Figure  11b  gives  the 
associated  histogram. 


i 


I 

I 

! 
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VUHbKB  OF  OBDeseO  PAIAS  » 200 


'SCAIB  : *•*  • 0.0S9B-01  UNITS 


* ccKCBCcsccKCCC««ccccccBrv:s:EEs;a:xrcc«:ct:i:ri:s;c«:x:t;ci:i:c¥:£- 


oeooooo 


I I 


I I 


009- 9^0^*^ 

r*  <N«tf 


f OOOOOO 


I I 


I I 


>1IUUUUM 

ASnOv'^OO 

9«c<N®r»® 

9*s«r^*no 

— 9'®**»-9‘ 

( III 

•,nn 

!/)M 

MV^ 

»*-«« 

««9X3b}{il 

ecv)»ao 


■ ■••11 

MWuaaM® 

<N(N'^0**«N 


«M»a  X 
«co  ^c> 
^ XCV> 
xo®«xa 

^tnu*"  ~ 


aooooao 

I I I I I I I 

UM  WMUUiM 
9 

^r« 

®<N^Or*<N^ 

^iN^o^ar* 


XMCcxBi: 

<OmOU(JX 

MtUXMMUj« 

CXHKSOX 


oo;  ■ stiff  aaiaoio  lo  taoBnt 


riccuiiciis  sinpie  size  • 2oo 


^OOCN*-0»» 

eoooeoo 

< I 

-o'^Oi^irtuir* 

000^004) 


t t t i 


V)  (M900<V9 

H 000009 

u uftiMMUioa 
S ^.9'£9<^9' 
9 ’OC*^*"'CsO 

B •“C*'*'l/'f*40 

■J  00“^»^nOP* 


U V)(A 


OB  Mtn 

au 


9 64aSt*C-> 

M ‘"i^Ksaiu 
s BCt/lBSO 


*-o«*o<-o 

oooooo 

MklhltalMtAJ 

<£Ou^t^Oio 

<N9>0«n»B 

oor»f*“ 

^>cor«<N^ 


a 

B a»  < 

a *kj»a  os 

« «o 

kl  M ft.s9'0 

m MOM^ata 

’A  xnocac 


^<N(N!NO 

ooooo 

till 

mwmuui 

m 9^m9' 

»» 

O 

a 

M • • . • • 

a 

a t I I I I 


M 

ftm(9 
« aacaca 

a w<« 

H a*"BRa 
a <c»Mua 
M 

(J  BK^-BB 
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L. 


47 


NUHDEft  or  OHOBEED  BAIRS 


HlLK'SHAriRO  TEST:  TEST  VALUE  » 9.6880-01  SEE  EIOBEtRlKA  (196SK  E.  S91,  FOR  CBItlCAL  VALUES 


rieCUlMClES  S&HPLE  SIZE 


M 


oeooeoe 

ill* 

MMKMMMM 

9»rNr*r>«ur>4f*^ 

QDm^Ou^^r^ 


u«u 


(St^o 

moM 


X UMMMM 


zassxz 

asaasss 

caocooE 


V)  xetnot/>OM 

a 


yj 

H OOOOOO 
a I I I I I I 
M MUMUMM 
u 

cj  'nr>*^tr'p*a 

C 

^ xA^O''e<« 
« 0<Nr«900 


£ 


U MM 

MM 

a MM 

w *0^<^ 

s aMM« 

a MSHH 

M "^.-vaaMM 

a ESMicea 


aeoaaa 

I I I I I 
ajMMMMM 

c 

-^<N\0<NM^ 


a 

M 


Ot 

M 


a 

tu  «»  a 

abi>a  a 
«a  aa 

M MaaM 

aaM<aa 


:»MU«aE 


OM— 

ooOQoeo 

t I I I I I 

MMMMWMM 

a'^r*4C>^<r«*» 

»•  Mf*4^o*o^r» 
U» 

a oaaairr«9* 

Ml  •••««•• 

a 


*j  saoMM 

« Ka<«aca 
a «uu]« 

H SMCcacM 
a <awow^ja 
u MuaMMua 
U CBHEEUa 


Figure  11b  - HISrOGRAM  OF  TH2  DATA  OF  FIGURE  11a 


**•  Using  * NORMPL*  with  Generated  Data  from  Various 
Distributions  . 

The  Subroutine  'NOHMEL'  has  been  used  here  to  test 
for  Horaality  some  data  sets  which  were  generated  ty  the 
computer  frcn  various  distributions.  The  plots  of  these 
data  sets  illustrate  the  advantages  and  disadvantages  of  the 
plotting  Subroutine  'NCRMPL'. 

Cbserving  Figures  7a-b,  8a-b,  9a-t  and  12a-c  we  may 
note  the  following; 


(a)  . Figure  12a  shows  a plot  of  a sample  of  5C  data 
points  from  a Normal  Distribution  ( N (1000,1)  ).  As  it  was 
expected  a very  straight  line  is  fitted.  In  addition  to 
that  plot,  the  Subroutine  ’NORMEL*  uses  the  option  to 
evaluate  the  ' WILK-SdAEIRO ' test  and  the  W-value  is  conputed 
and  printed  on  the  same  figure.  Comparing  this  value  with 
the  percentage  points  of  the  W-test  which  are  given  in  [22], 
for  sample  size  n=50  we  see  that  we  may  accept  the  normality 
assumption  with  a significanse  level  = .02  . Figure  12b 
gives  the  associated  histogram. 

(b)  . Figure  7a  plots  a sample  of  200  data  points 
of  the  Unit  Exponential  distribution.  Here  the  nonlinear 
graph  is  obvious  ( as  is  expected  ) and  suggests  tc 
reject  the  normality  assumption  of  the  data  without  doing 
any  further  formal  test. 


(c) . Figure  8a  gives  a plot  of  a sample  of  size 


200  from  Unirorm  ( 0,  1 ) 
that  plot  indicates  data  not 
therefore  the  user  is  told 
the  same  time  a suggestion 


generated  data.  The  shape  of 
as  dispersed  as  the  ncrcal, 
to  reject  normality,  getting  at 
that  data  may  ne  Uniformly 


50 


distrihuted  ( because  of  that  particular  shape  of  the 
graph  ) . 


(d) , Figure  9a  suggests  strongly  departure  frcm  a 
assuapticD.  Indeed  this  set  of  data  have  keen 
from  a Cauchy  Distribution  and  its  shape  really 
a symmetric  distribution. 


i 

i 

\ 


Normality 

generated 

indicates 


5-  Osin^  *HORMPL*  With  Cost  Overruns  Eat a 


•NOBflPL'  was  used  with  the  cost  overruns  data  [6], 
for  the  year  1950  to  see  if  there  is  a linear  fit.  Locking 
at  Figure  13  we  see  that  data  may  be  Normally  distributed 
since  the  plot  is  very  close  to  linear.  Of  course  some  data 
points  deviate  from  the  straight  line,  but  these  points  may 
be  considered  as  outliers. 


! 


54 


MUHBEB  or  OBDEBED 


-SCALE  : « 0.3SitB  00  UHITS 


D.  SOBaCOTINE  EXPIT 


1 . Description 

Subroutine  'EXPLT'  is  indended  tc  be  usee  for 

testing  a saiiple  of  independent  observations  of  a continuous 

random  variable  X with  a distribution  function  which  is 

believed  to  be  Exponential,  by  plotting  each  of  the  Order 

Statistics  X versus  its  Expected  Value,  and  by  estioating 

U) 

the  parameters  '^^^and  'the  observations  ( and 

introduced  later) . 

This  Subroutine  sorts  the  data  into  increasing 
order,  obtaining  the  Order  Statistics.  The  Order  Statistics 
are  then  plotted  versus  Exponential  Scores.  These  scores  are 
evaluated  by  the  Subroutine.  The  Exponential  Scores  are  the 
expected  value  of  the  Order  Statistics  ( E[ X ] ).  They 

(i) 

are  more  easily  derived  and  computed  than  the  Normal  Scores. 
That  is: 

Let  X be  a random  variable  Exponentially  distributed 
with  parameter  >=  1 
Then  it  is  known: 

--U 

f (X)  = V , X > 0 , 

A 

->x 

F^(x)  = 1-e  , ^ > 0* 
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Ihe  Density  cf  the  i Order  Statistic  is  given  by: 

n-1  i-1  n-i 

f (X  ) = n(  ) F (X  ) (1-F  (X  ))  f (X  ) 

I <i)  i-1  X <i>  X <i)  X <i) 

ti> 


Ihus  for  the  1 Order  Statistic  we  have: 


n-1  Cl) 

(x^^  = n(1-  (1-e  ) ) ;>e 

M) 


-njkx 


= nxe 


Therefore  X is  exponentially  distributed  with 
(1) 


paraaeter  ">  = n > . 
and  for  %=  1 we  have: 


1 1 

S[X  3 = — 

(1)  n ^ n 


It  can  be  proved  (see  Feller  [7])  that  the  n random 

variables  X , X -X  , X -X  -X  from 

(1)  (2)  (1)  (3)  (2)  (n)  (n-1) 

an  Exponential  random  variable  X with  density  function: 


f^  (X)  = ye  , 


X > 0 


are  independent  and  the  density  of  the  random  variable 

1=1  - 1 is  given  by: 

(k)  (k-1) 


f^(y)  = (n-k  + 1)  ^e 


- :\(n-k+1)  y 


y > 0 


That  is,  X is  an  Exponential  random  variable  with 
parameter  (c-k+l)  ^ . 
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Ibetefore  B|;  Y ] is  given  by; 


1 1 

2[  X -X  ] , 

(k)  (k-1)  :;^(n-k  + 1)  n-k+1 


for  = 1 


Ibis  fact  is  used  to  derive  EFX  ],  i = 2,3,. ..,n  . 

(i) 

Ibat  is: 

E[X  ) = EC  (X  -X  ) + (X  -X  )+...♦ 

(i>  (i)  (i-1)  (i-1)  a-2) 

(X  -X  ) ♦ X ] 

(2)  (1)  (1) 


= sex  ] + E[X  -X  ] + ...+ 

(1)  (2)  (1) 

E[X  -X  ] + E[  X -X  ] 

(i-1)  (i-2)  (i)  (i-1) 


1 1 1 
n n-2+1  n-3+1 


n-i+1 


2 

k=1  n-k+1  ' 


i * 2,...n 


•EXPLT'  derives  the  values  of  E[ X ] using  the 

(i) 

aoove  formula  and  then  it  uses  its  Subroutine  'SPLCZ'  to 
plot  the  Order  Statistics  versus  their  expected  value. 


In  addition  to  the  plot,  which  is  an  informal  and 
quick  test  for  exponentiality  the  estimates  of  the 
statistics  and  provide  another  informal  test  since  the 
values  of  and  have  a constant  value  for  any 
Exponential  distribution; 


Ey  definitioE  ^j^and  are  given  by: 


(X-r) 


ii  = Ef- 


(X-Y) 


3 - 3 


lb  US 


= ---  EC  (X-Y)^ 


1 3 2 2 3 

--  £[X  - 3X  y'  ♦ 3Xy^  - ] 


(1) 


1 3 2 2 3 

--  (E[X  ] - 3^E[X  ] + 3 y E[XJ  -y  ) 


Now 


E[X^] 


CO 


X f (X)  dx 
^ X 

«c 


r 3 -:\x 

= 3o  * dx 


r(«) 


•.4-1  4 ->-\X 

X € dx 


Ibe  above  integral  is  that  of  tbe  Gaiaaa  density  fucction 
with  parameter  oC  = 4 , and  since  its  value  equals  1 we 
have: 


3 31  6 

E[X  ] 


1 
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2 

Ccoputing  E[X  ] «e  have: 


2 

E[X  ] = 

e dz 

''O 

r(3) 

3-1  3 -;vx 

r X e 

1 r{3) 

= 

2 

• 

ECX]  = 

_L_ 

“oT 

• 

Aisc 

Substituting  in  (1)  we  have: 


3 


16  12  11 


3 3 

= y { 2/^  )=  2 . 


Iherefore  the  value  of  )|  ^ equals  2 for  any  Ex 
distributicc . 

Also  is  coaputing  as  follows: 


4 


4 3 2 2 3 4 4 

= EC  JI  - 4y*X  ♦6yX  - hy4X*y]/C‘  - 3 


4 3 2 2 3 4 

(2)  = (ECX  ] - 4yECX  ]*6yE[X]-4yE[X]4y 


on  en  tial 


4 

} /<7-3. 
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Therefore  the  value  of  equals  6 for  any 

Expcnential  distrihuticn. 

This  fact  is  used  by  the  Subroutine  'EXPLI'  tc  give 
the  user  cne  sore  informal  test  for  expcnentiality . It 
estimates  the  values  of  and  ^ 

and/or  ^ 6 then  the  data  may  not  have  an  Expcnential 
distriburior:. 

Estimating  the  values  of  and  , 'EXPLT*  uses 

3 a 

unbiased  estimator  fcr  c , E[  (X-^)  ] and  E[  (X-y)  ] using 

the  formulas: 

n 2 1/2 

a = (.r  (X.-X)  / (n-l)  ) 

1-1  1 

n _ 3 

3 1=1  1 

EC  (X-V)  ] = n 

' (n-1)  (d-2) 


EC  (X-y)  ] 


n _ 4 

(n(n-2)*3)  (,X  (X.-X)  ) 
1=1  1 


(n-1)  (n-2)  (n-3) 
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3 (T  (n-1)  (2n-3) 

n (n-2)  (n-3) 


1 


2.  2MEiIl  structure 

Ihe  Subroutine  'EXEIT'  is  also  a ECRTEAN  callable 
subroutine  with  each  call  producing  a plot  fcr  the  given  set 
of  data  and  at  rhe  same  time  calculating  and  printicg  the 
estimated  values  of  the  parameters  '^^and  cf  the  data. 

Basically  the  program  is  divided  intc  two  parts: 


Ihe  firsr  parr  is  involved  with  all  required 
computations  calculating  the  Exponential  Scores  and  the 
various  statistics  needed  to  compute  the  ^ and 
paramete  rs. 

Ihe  second  part  is  the  Subroutine  'EELOT*.  'EELQT* 
scales  tne  data  points  according  to  their  range  and  plots 
the  scaled  data  along  110  equal-spaced  positions  of  the 
X-axis.  No  plot  is  given  if  the  data  have  constant  value  or 
if  there  exists  a data  point  less  than  zero. 

A complete  description  of  how  'EXELI*  operates  is 
given  by  the  Subroutine  and  a summary  description  is  given 
on  the  terminal  by  typing  the  command  DESCRIES  EXPLT  under 
the  cms  environment.  When  the  user  types  DESCHI3E  EXELT 
the  following  response  is  printed  on  the  terminal: 

SDBfiCUTINE  EXPLT 

•EXEIT'  takes  a set  of  data,  scrts  it  into 
increasing  order  and  plots  the  created  order  statistics 
versus  exponential  scores  (expected  values  of  order 
statistics)  computed  by  : 


( 

L 
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E[  X ] = ( i/  (n~k+1)  )/  1 = 1/2f«a<#n 

(1)  ^ 
k = 1 


MsCf  it  computes  the  estimates  of  ^ and 
parameters. 


It  is  called  by  : 

CALL  EXELT  ( X,  SCCBES,  N ) 

ASGOEENTS 

X Is  the  array  containing  the  data 

SCOBES  Is  a work  array  of  dimension  N 
N Is  the  number  of  data  values 


acre  information  is  given  in  the  subroutine. 


^ • Internretinq  The Output 

If  data  have  an  Exponential  distribution  then  the 
plot  will  tend  to  a straight  line.  But  it  should  be  noted 
that  the  inverse  is  not  generally  true.  The  linearity  of  the 
graph  is  an  indication  only  and  gives  the  user  a first 
feeling  of  the  distribution  of  data.  If  the  plot  is  not 
linear,  however,  this  suggests  that  they  are  not 
Exponential.  In  addition  the  shape  of  the  plot  may  lead  the 
user  to  make  an  alternative  model  selection. 


Cn  the  other  hand  the  values  of  the 
statistics  jj  and  may  be  used  to  test  informally  tne 
distribution  of  data.  These  values  are  an  indication  only 
for  the  test  for  expcnen tiality . If  the  estimate  of  has 
the  value  about  2 and  rhe  estimate  of  has  the  value  about 
6 then  it  suggests  that  data  may  be  exponentially 
distributed.  But  a departure  of  these  values  suggests  that 
data  do  not  have  an  exponential  distribution. 
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I 


'*  • i]§iS9  ' EXPLT*  with  Data  Generated  from  Various 
listr ifcutiopg 


The  Subroutine  ’EXPLI*  has  been  used  here,  plotting 
data  generated  by  the  computer  from  various  distrituticns, 
in  order  to  get  a visual  sense  of  the  behavior  of  the  plot. 

figure  14a  gives  tbs  plot  of  ICC  Oniform  (0,1) 
variates.  Obviously  nonlinearity  governs  the  plot  (as  we 
should  expect)  and  thus  it  suggests  departure  from  the 

Exponenrial  distribution.  Besides  the  nonlinearity  cf  the 
plot  the  estimates  of  ^^and  '^2.  (O*!  -1.15  respectively) 

are  far  away  from  the  ccrres pending  values  (2  and  6)  of 
exponentially  distributed  data.  Figure  14b  gives  the 
associated  histogram. 

Figure  15a  plots  variates  generated  from  a 

Triangular  Symmetric  distribution  (0,2).  The  nonlinearity 
and  the  values  of  the  estimates  and  suggest  that  one 

reject  the  exponential  assumption. 

2 

Figure  16  shows  the  plot  of  X generated  variates 

with  10  degrees  of  freedom.  For  the  same  reasons  as  above 
the  exponential  model  is  rejected. 

Figure  17a  plots  Normally  (N  (1000,1))  distributed 
generated  data  and  the  nonlinearity  of  the  plot  is  apparent. 
Besides  the  nonlinearity,  the  estimates  of  and 

(.09  and  -0.63)  suggest  not  only  tc  reject  the 

Exponentiality  of  data  but  in  addition  lead  us  to  test  the 
data  for  Normality.  Figure  17b  gives  the  associated 
histogram. 


L 
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figure  18  plots  exponentially  generated  variates 
«ith  paraaeter  '^=  1.  As  we  see  the  graph  is  linear  anc  the 
user  gets  the  idea  that  the  data  possibly  are  exponertial. 
Besides  the  linearity  here,  the  values  or  the  estimates  of 
and  (2.06  and  5.4)  are  close  enough  to  their  true 

values  to  indicate  an  Exponential  model. 
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Figure  14b  - dlSrOGSAM  OF  THE  DATA  CF  FIGURE  14a 
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Figure  17a  - 


A SAMPLE  OF  SIZE  100  FROM  A NORMAL  (1000,1) 
DISTRIBUTION 
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BSlIBiTEO  PARAflETEBS  OP  DATA  : GAHHA1  « 2.061009E  00  CAaRA2  « S.ll0a789B  00 


Figure  18b  - HISTOGRAM  CF  THS  DATA  CJ  FIGURE  16a 


A.  DESCGIfllCM 


I 

t 

\ 


Another  Subroutine  presented  in  this  thesis  is  the 
Subroutine  'LIST'.  Its  purpose  is  to  list  a given  set  of 
data  in  ascending  order  and,  taking  advantage  of  like 
occurrences  in  the  data,  to  print  the  ordered  data  in  a 
compressed  form.  This  feature  becomes  highly  useful  uhen 
listing  a large  number  of  data  points  that  contain  many 
repeated  values.  It  is  also  a tccl  for  finding 
multiplicities  in  supposedly  continuous  data,  and  a 
probability  function  estimating  routine  fcr  data  which  is 
known  to  he  discrete. 

A complete  description  cf  how  'LIST'  cperares  is  given 
in  the  Subroutine.  However  a summary  is  given  by  typing 
CESCHIBE  LIST.  When  the  user  types  the  ccmmand 
EESCHIBE  LIST  under  the  CMS  environment  the  follcwing 
response  will  be  given  on  the  terminal: 

SUBaCUTIKZ  LIST 

•LIST'  scrts  a set  of  data  into  increasing  order  and 
gives  a 5'cclumn  print-out  as  follows: 

1st  column:  Serial  number  cf  the  first  occurrence  of 
this  value  in  the  ordered  List 

2Qd  column:  Value  of  ordered  data-value 

3rd  column:  frequency  cf  occurrence  cf  the  value 
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4th  cclumn;  Percent  for  the  value 


5th  cclusn:  Graphical  representation  cf  the  frequency 
for  each  value. 

It  is  called  by: 

CALI  LIST  ( X,  N ) 
where: 

X Is  the  array  of  data 
N Is  the  nuuijoer  of  data-values. 

More  information  is  given  in  the  subroutine. 

E.  INTEEE3EII1SIG  THE  OUTPUT 

The  print-out  is  a visual  rapresentaticn  of  the  data  and 
of  each  data  value  frequency.  For  data  points  having  the 
same  value  'IIST'  will  print  this  value  once  with  the  number 
cf  occurrences.  A 5-column  output  will  be  printed  ard  its  j 

interpretation  is:  i 

1.  First  column  gives  the  serial  number  of  the  first  i 

occurrence  of  this  value  in  the  ordered  List 

2.  Second  column  gives  the  value  cf  ordered  data-value 

3.  Third  column  gives  the  frequency  cf  occurrence  of 
the  value 

4.  Forth  column  gives  the  percent  for  the  value 

5.  Fifth  column  is  a graphical  representation  cf  the 
frequency  for  each  value. 


W 
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Zzaagle; 


Let  1,  1,  5,  3,  5,  1,  2,  2t  5,  6 be  a given  data  set. 

Then  the  crdered  data  will  be:  X =X  =X  = 1» 

(1)  (2)  (5) 

X =X  * 2,  X = 3,  X =X  =X  = 5,  X =6.  Then 
(4)  (5)  (6)  (7)  (8)  (9)  (10) 

the  print-out  will  have  the  following  form: 


SEH^A;.  MUMBEE OBDEBED  data FREOaEMCY P2BCENT PBOE.  , GBAPH 

1 1 3 .3  *** 

U 2 2 .2  ** 

6 3 1 ,1  * 

7 5 3 .3 

10  6 1 .1  » 

If  there  are  no  data-pcints  having  a coaocn  value  then 
•LIST'  gives  this  indication  and  prints  cnly  the  crdered 
data.  This  happens  when  data  have  a continuous 
distribution. 


C.  OSISG  ‘IISI'  WITH  TEL2PECNE  DATA  1 AND  TELEPHONE  LATA  2 

Subroutine  'LIST'  was  used  with  Telephone  Data  1 
(Figure  IS)  and  Telephone  Data  2 (Figure  20)  and  a brief 
analysis  of  the  output  follows. 

Looking  at  Figure  19  and  Figure  20  we  may  get  some 
infornaticn  of  each  data  set.  It  can  be  seen  that  both  of 
the  data  sets  contain  a large  number  of  multiple  occurrences 
of  the  data  value  one  and  data  value  two.  As  we  can  see 
reading  the  fourth  column  the  occurrence  of  ones  is  19!*  for 
Telephone  Data  1 and  24%  for  Telephone  Data  2. 
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Comparing  also  Figures  19a-c,  20d-d  ve  can  see  that 
multiple  occurrences  happen  in  the  range  1 to  24C  for 
Telephone  Data  1 and  in  the  range  1 to  132  for  Telephone 
Data  2. 

Furthermore,  guick  visual  informaticn  concerning  the 
range  where  ve  have  multiple  occurrences  can  be  obtained 
from  the  Probability  Graph.  Thus  we  can  see  that  for 
Telephone  Data  2 there  is  a region  from  113  to  132  where  the 
multiple  occurrences  of  values  is  larger  than  in  neighboring 
regions.  Therefore  'LIST*  gives  the  user  a useful  tccl  for 
analysis  cf  data. 
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03 

1 

0.00149 

537 

0.  15300000E 

03 

1 

0.00149 

538 

0.  15600000E 

03 

1 

0.00149 

539 

0.  ISdOOOOOE 

03 

1 

0.00149 

540 

0.  16100000B 

03 

1 

0.00149 

541 

0.  16500000E 

03 

1 

0.00149 

542 

0.  17500000E 

03 

1 

0.00149 

543 

0.  17600000E 

03 

1 

0.00 149 

544 

0.  17700000E 

03 

1 

0.00149 

545 

0.  18300000E 

03 

1 

0.00149 

546 

0.  laeoooooE 

03 

1 

0.00149 

547 

0.  19700000S 

03 

1 

0.00149 

548 

0.  1920000  OB 

03 

1 

0.00149 

549 

0.  19300000E 

03 

1 

0.00149 

550 

0. 20200000E 

03 

1 

0.00149 

551 

0. 21700000E 

03 

2 

0.00298 

553 

0. 22400000E 

03 

1 

0.00149 

554 

0. 22600000E 

03 

1 

0.00149 

555 

0. 22800000E 

03 

3 

0.00446 

558 

0, 23100000E 

03 

1 

0.00149 

559 

0.  23400000E 

03 

1 

0.00149 

560 

0. 2J700000B 

03 

1 

0.00149 

561 

0. 2J900000B 

03 

2 

0,00298 

563 

0. 24000000B 

03 

4 

0.00595 

567 

0. 24100000E 

03 

1 

0.00149 

568 

0. 24400000E 

03 

1 

0.00149 

569 

0. 24800000E 

03 

1 

0.00149 

570 

0. 24900000E 

03 

1 

0.00149 

571 

0. 25100000E 

03 

1 

0.00149 

572 

0. 252000008 

03 

1 

0.00149 

573 

0.  27000000S 

03 

1 

0.00149 

574 

0. 27900000E 

03 

1 

0.00  149 

575 

0.286QOOOOE 

03 

1 

0.00149 

576 

0. 2S700000E 

03 

1 

0.00149 

577 

0.30300000E 

03 

1 

0.00149 

578 

0. 31200000E 

03 

1 

0.00149 

579 

0. 31800000E 

03 

1 

0.00149 

580 

0. 34700000E 

03 

1 

0.00149 

581 

0. 35400000E 

03 

1 

0.00149 

582 

0. 3b000000E 

03 

1 

0.00149 

583 

0.  36400000S 

03 

1 

0.00149 

584 

0. 37Q00000E 

03 

1 

0.00149 

585 

0.39000000S 

03 

1 

0.  00  149 

586 

0. 39400000E 

03 

1 

0.00149 

587 

0.  435000006 

03 

1 

0.00149 

588 

0. 46000000E 

03 

1 

0.00149 

589 

0.  46500000B 

03 

1 

0.00149 

590 

0.47300000E 

03 

1 

0.00149 

591 

0. 48000000S 

03 

1 

0.00149 

592 

0. 491000008 

03 

1 

0.00149 

593 

0. 54900000E 

03 

1 

0.00149 

594 

0. 6COOOOOOE 

03 

2 

0.00298 

596 

0. 60900000E 

03 

1 

0.00  149 

597 

0. 62  100000E 

03 

1 

0,00149 

598 

0. 711000008 

03 

1 

0.00149 

599 

0.  8 1700000B 

03 

1 

0.00149 

600 

0. 83600000E 

03 

1 

0.00149 

601 

0. 92700000B 

03 

1 

0.00149 

602 

0.  1 1240000E 

04 

1 

0.00149 

603 

0.  1 1500000E 

04 

1 

0.00149 

604 

0.  12700000E 

04 

1 

0.00149 

605. 

0.  128900008 

04 

1 

0.00149 
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606 

0. 12980000E 

607 

0. 13050000E 

608 

0.  1328000 OE 

609 

0,  13480000B 

610 

0.  1355000CE 

611 

0,  14120000E 

612 

0.  14290000E 

613 

0.  14890000B 

614 

0.  14930000B 

615 

0.  15100000E 

616 

0.  15190000E 

617 

0.  15470000E 

618 

0.  163300U0E 

619 

0.  17370OU0S 

620 

0.2C720000S 

621 

0.  24830000E 

622 

0.  280600008 

623 

0.  29620000E 

6 24 

0.  3026000UE 

625 

0.  32010000S 

6 26 

0.  359300008 

627 

0.  36850000E 

6 28 

0. 3S5200C0E 

629 

0.  4 1570000E 

6 30 

0.  44690000E 

631 

0.  620600008 

632 

0.  76  140000E 

633 

0.  832200008 

634 

0.9C150000S 

635 

0.962500008 

636 

0.  980600008 

637 

0. 981600008 

6 38 

0.  10  154000B 

639 

0. 1039800DE 

640 

0.  104510008 

641 

0. 109390008 

642 

0.  1 12800008 

643 

0.  1 344700  08 

644 

0.  143850008 

645 

0,  151550008 

646 

0.  152940008 

647 

0.  15504000B 

648 

0. 156470008 

649 

0.  156580008 

650 

0.  162800008 

651 

0.  162990008 

652 

0.  163610008 

653 

0.  164030008 

6 54 

0.  168170008 

655 

0.  171740008 

656 

0.  176670008 

657 

0.  182160008 

658 

0.  186490008 

659 

0.  194610008 

660 

0. 216460008 

661 

0. 234990008 

662 

0. 246920008 

663 

0. 264430008 

664 

0. 3C9740U08 

665 

0. 356440008 

666 

0. 380030008 

667 

0. 40131000E 

668 

0. 471200008 

669 

0. 47592Q00S 

670 

0. 6 17100008 

671 

0,697750008 

672 

0. 85993C00B 

0« 


1 

1 

1 

1 

1 

1 

1 

t 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0-00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0-00149 
0. 00149 
0.00149 
0-00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00149 
0.00  149 
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SZaXil  SOBBBB 

OICBBBD  DBTi 

FBBQOBKCZBS 

PE8CSRT 

pBQBiam 

1 

0.  10000000B 

01 

178 

0. 24185 

179 

0. 20000000B 

01 

36 

0.04891 

215 

0.30000000B 

0 1 

n 

0.01495 

•• 

226 

9.400000008 

01 

6 

0.00315 

m 

2i2 

0.  500000006 

01 

6 

0.00815 

• 

238 

0.  600000 OOB 

01 

5 

0.00679 

• 

2'*J 

0.  700000006 

0 t 

5 

0.00679 

• 

240 

0.  800000006 

01 

4 

U. 00543 

• 

252 

0.  90000000B 

01 

4 

0.00543 

« 

256 

0. 100000006 

9 

0.01223 

• • 

Q. llQQQOOOe 

^ • 

2 

0.00272 

0. 120000006 

02 

3 

0.00408 

• 

270 

0.  13000000B 

02 

1 

0.00136 

271 

0. 140000006 

02 

1 

0.00136 

0. 150000006 

1 

0.00136 

273 

0. 160000006 

02 

4 

0.00543 

• 

277 

0.  190000006 

02 

1 

0.00136 

278 

Q. 2 10000006 

02 

1 

0.00  ? 36 

279 

0. 220000006 

02 

1 

0.001 J6 

200 

0. 240000006 

02 

3 

0.00408 

• 

283 

0. 250000006 

02 

3 

0.00408 

• 

286 

C. 260000006 

02 

2 

0.00272 

238 

0.  270000006 

02 

1 

0.00136 

289 

0. 3C0000006 

02 

3 

0.00408 

292 

0. 320000006 

02 

1 

0.001 36 

29J 

0.  J60000006 

02 

1 

0.00 136 

294 

0.  UOUOOUOOS 

02 

2 

0.00272 

296 

c. a2aououos 

02 

1 

0.001 36 

297 

a. 4ujojuoas 

02 

1 

0.00136 

298 

0. 470000000 

02 

4 

0.00543 

« 

302 

0. 430000006 

02 

1 

0.00136 

30J 

0.490000006 

02 

9 

0.01223 

• • 

0.500000J06 

02 

4 

0.00543 

• 

316 

0.  5 300COU0B 

02 

3 

0.00408 

• 

319 

0.550000006 

02 

2 

0.00272 

321 

0.560000002 

02 

4 

0.00543 

• 

325 

0.  570000005 
0. 560000006 

02 

3 

0.00408 

• 

328 

02 

1 

0.00 1 36 

329 

0.590000006 

02 

3 

0.00408 

« 

332 

0.600000006 

02 

2 

0,00272 

334 

0.620000006 

02 

3 

0.00400 

• 

3J7 

0.630000006 

02 

1 

0.001 36 

338 

C. 6UOOOG006 

02 

4 

0,0054  3 

« 

342 

0.O50000U06 

02 

3 

0.00400 

m 

345 

0.  660000006 

02 

4 

0.00543 

« 

349 

0.670000006 

02 

1 

0.001 36 

350 

0. 600000006 

02 

5 

0.00679 

• 

355 

O.lj'fOOOJJUl! 

02 

3 

0. 00400 

• 

358 

0. 700000005 

02 

3 

0.00408 

• 

361 

0.710000006 

02 

5 

0.00679 

• 

366 

C. 730000006 

02 

1 

0.00136 

367 

0. 740000005 

02 

1 

0.00136 

368 

0.  T50000J0E 

02 

1 

0.00136 

369 

Q. 760000005 

02 

1 

0.00136 

no 

0. 770000006 

1 

0.00136 

371 

0. 810000006 

02  - 

1 

0.001 36 

372 

0.  0bOOOOOOe 

02 

1 

O.UO 136 

373 

0. 380000006 

02 

1 

0.00136 

374 

0.890000006 

02 

1 

O.OOT36 

375 

0. 93000C006 

02 

1 

0.00136 

376 

0.950000006 

02 

1 

0.001 36 

377 

0. 960000006 

02 

1 

0.001 36 

378 

0. 970000006 

02 

1 

0.00136 

379 

0. 990000006 

02 

3 

0.00408 

■ 

382 

0.100000006 

03 

1 

0.00136 

383 

0. 101000006 

03 

1 

0.001  36 

304 

a. losooaooE 

03 

1 

0.001 36 

335 

0. 106000006 

03 

1 

0.001 36 

306 

0. 107000006 

03 

1 

0.001 jb 

307 

0. icsoooooe 

03 

1 

0.00 1 36 

388 

0. 1 10000006 

03 

1 

0.00136 

389 

0. 1 1000006 

03 

1 

0.001 36 

390 

0. 1 12000006 

03 

1 

0.001 36 

391 

0.  1 1300000E 

03 

4 

0.00543 

• 

395 

0. 1 14000006 

0 3 

3 

0.00408 

• 

398 

0. 1 15000005 

03 

1 

0.001 36 

399 

0. 1 16000005 

03 

4 

0.0054J 

« 

*03 

0. 1 17000006 

03 

1 

O.OOT36 
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r 

404 

0.  11800000E 

03 

5 

0.00679 

« 

409 

0.  1 1900000E 

03 

8 

0.01087 

• • 

417 

0,  12000000E 

03 

9 

0.01223 

426 

0.  12100000E 

03 

7 

0.  00951 

«« 

433 

0.  12200000E 

03 

12 

0.01630 

«*« 

443 

0.  12JOOOOOE 

03 

8 

0.01087 

453 

0.  12400000E 

03 

7 

0.00951 

460 

0.  12500000B 

03 

9 

0.01223 

«* 

469 

0.  12600000E 

03 

6 

0.00815 

» 

475 

0, 12700000E 

03 

3 

0.00408 

♦ 

478 

0.  12800000E 

03  ' 

2 

0.00272 

, 

480 

0.  12900000S 

03 

4 

0.00543 

; 

484 

0.  13000000E 

03 

3 

0.00408 

487 

0.  ’3200000E 

03 

3 

0.00403 

490 

0.  1 3300000E 

03 

2 

0.00272 

492 

0.  13700000E 

03 

1 

0.00136 

1 

493 

0.  14000000E 

03 

1 

0.00136 

494 

0.  1430000US 

03 

1 

0.001 36 

495 

0.  15200000B 

03 

1 

0.00136 

496 

0.  15800000E 

03 

1 

0.00136 

497 

0.  16600000E 

03 

1 

0.00136 

498 

0.  16900000E 

03 

2 

0.00272 

500 

0.  noooooop. 

03 

1 

0,00136 

501 

0.  17300000E 

03 

2 

0.00272 

503 

0.  17600000B 

03 

1 

0.00136 

504 

0.  178000005 

03 

1 

0.00136 

505 

0.  1800000QE 

03 

2 

0.00272 

507 

0.  18200000S 

03 

1 

0.001  36 

508 

0.  18500000E 

03 

1 

0.00136 

509 

0.  19700000E 

03 

1 

0.00136 

510 

Q.  19000000E 

03 

1 

0.00136 

511 

0,  19U00000E 

03 

1 

0.00136 

512 

0.  19500000E 

03 

1 

0.001 36 

513 

0.  19900000E 

03 

1 

0.00136 

514 

0. 20600000S 

03 

1 

0.00136 

515 

0. 20900000E 

03 

1 

0.00136 

0.21600000S 

03 

1 

0.00136 

517 

0. 22900000E 

03 

1 

0.00136 

518 

C. 23000000E 

03 

1 

0.00136 

519 

0. 23700000E 

03 

1 

0.00136 

520 

0. 2J900000E 

03 

3 

0,00408 

• 

523 

0.  24000000E 

03 

2 

0.00272 

525 

0,241000005 

0^ 

1 

0. 00 1 36 

526 

0. 24400000E 

03 

2 

0.00272 

528 

0. 24700000E 

03 

2 

0.00272 

530 

0.  24800000E 

03 

2 

0.00272 

532 

0. 25000000E 

03 

1 

0.00136 

533 

0. 25100000E 

03 

1 

0.00136 

534 

0. 25400000E 

03 

2 

0.00272 

536 

0.  25500000E 

03 

2 

0.00272 

538 

0. 2S60000UE 

03 

1 

0.00136 

539 

0. 25300000E 

03 

1 

0.00136 

540 

0. 27200000E 

03 

1 

0.00136 

541 

0. 27500000E 

03 

1 

0.00136 

542 

0. 28000000E 

03 

1 

0.001  36 

543 

0.29300000E 

03 

1 

0.00136 

544 

0.  3 1300000E 

03 

1 

0. 001 36 

545 

0. 32400000E 

03 

1 

0.00136 

546 

0. 34000000S 

03 

1 

0.001  36 

547 

0. 34400000E 

03 

1 

0.001  36 

548 

0. 34600000E 

03 

1 

0.001 36 

549 

0.35900000E 

03 

1 

0.00136 

550 

0. 36500000E 

03 

2 

0.00272 

552 

0. 36700000E 

03 

1 

0.001 36 

553 

0. 37700000E 

03 

1 

0.00136 

554 

0. 38000000E 

03 

1 

0.001 36 

555 

0, 38700000E 

03 

1 

0.00136 

556 

0, 41300000E 

03 

1 

0.001 36 

557 

0. 4 1800000E 

03 

1 

0.00136 

558 

0.42900000E 

03 

1 

0.00136 

559 

0.45300000E 

03 

1 

0.00136 

560 

0.  46600000E 

03 

1 

0.00136 

56  1 

0. 47200000E 

03 

1 

0.00 1 36 

562 

0. 50300000E 

03 

2 

0.00272 

564 

0. 5C4000U0E 

0 3 

1 

0.001 36 

565 

0. 50700000E 

03 

1 

0.00 136 

566 

0.5  1000000E 

03 

1 

0.00  1 36 

567 

0. 5140000UE 

03 

1 

0.001 36 

568 

0.  5 nOOOOOE 

03 

1 

0.  JO  1 36 

. 

569 

0. 528000005 

03 

1 

0.00136 
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570 

0. 57000000B 

03 

1 

0.00136 

571 

0.57200000E 

03 

1 

0.00136 

572 

0.58600000E 

03 

1 

0.00136 

573 

0. 58700000E 

03 

1 

0.001 36 

574 

0. 50100000E 

03 

1 

0.00136 

575 

0. 6 1400000B 

03 

1 

0. 00136 

576 

0. 62400000E 

03 

1 

0.00136 

577 

a.64Q00000E 

03 

1 

0.00136 

578 

0.663000U0E 

03 

1 

0.00136 

579 

0. 71400000E 

03 

1 

0.00 1 36 

580 

0. 71600000E 

03 

2 

0.00272 

582 

0.74000000E 

03 

1 

0.00136 

583 

0. 742000QOE 

03 

2 

0.00272 

585 

0.79100000E 

03 

1 

0.00136 

586 

0. 85300000E 

03 

1 

0.00136 

587 

0. 86200000S 

03 

1 

0.001 36 

588 

0. 86600000E 

03 

1 

0.001 36 

589 

0. 87500000E 

03 

1 

0.001 36 

590 

0. 88500000E 

03 

1 

0.00136 

591 

0. 89400000E 

03 

1 

0.001 36 

592 

0. 92600000E 

03 

1 

0.00136 

593 

0. 92800000E 

03 

1 

0.00136 

594 

0. 94 100000E 

03 

1 

0.  00136 

595 

0. 97600000E 

03 

1 

0.001 36 

596 

0.99100000E 

03 

1 

0.00  136 

597 

0-  10210000S 

04 

1 

0.001 36 

598 

0.  10230000S 

04 

1 

0.00136 

599 

0.  1C620000S 

04 

1 

0.  00 1 36 

600 

0.  108600002 

04 

2 

0.00272 

602 

0.  110700002 

04 

1 

0.001  36 

603 

0.  1 1 190000E 

04 

1 

0.00136 

604 

0.  1125000  02 

04 

1 

0.00136 

605 

0.  123  10000E 

04 

1 

0.00136 

606 

0.  125000002 

04 

1 

0.001  36 

607 

0.  12610000E 

04 

2 

0.00272 

609 

0.  129600002 

04 

1 

0.00136 

610 

0.  136700002 

04 

1 

0.00136 

6 11 

0.  141200002 

04 

1 

0.00 1 36 

612 

0.  1 41  300002 

04 

1 

0.001 36 

6 13 

0.  14620OJ0E 

04 

1 

0.00136 

614 

0.  148500002 

04 

1 

0.00136 

6 15 

0.  149700002 

04 

1 

0.00136 

6 16 

0-  150400002 

04 

1 

0.00136 

617 

0.  151700002 

04 

1 

0.001 36 

6 18 

0.  153400002 

04 

1 

0.00136 

619 

0.  1535000  CE 

04 

2 

0.00272 

621 

0.  156200002 

04 

1 

0.00136 

622 

0.  159400002 

04 

1 

0.00136 

623 

0.  159500002 

04 

1 

0.001  36 

624 

0.  160300002 

04 

1 

0.00136 

625 

0-  165200002 

04 

1 

0.00136 

626 

0.  165300002 

04 

1 

0.00136 

627 

0.  16930000B 

04 

1 

0.00136 

628 

0.  nOOOOOOE 

04 

1 

0.00136 

629 

0.  1715000CB 

04 

1 

0.00136 

630 

0.  17500000E 

04 

1 

0.001  36 

631 

0, 17550000E 

04 

1 

0.00136 

632 

0.  173600002 

04 

1 

0.00136 

633 

0.  182200002 

04 

1 

0.00136 

634 

0.  1824000  OB 

04 

1 

0.00 1 36 

635 

0.  185500002 

04 

1 

0.00136 

636 

0. 187700002 

04 

1 

0.001 36 

637 

0.  189300002 

04 

1 

0.00136 

638 

0.  189700002 

04 

1 

0.001 36 

639 

0. 19230000E 

04 

1 

0.00136 

640 

0.  1S640000E 

04 

1 

0.001 36 

641 

0.  19650000E 

04 

1 

0.001  36 

642 

0.  205  100002 

04 

1 

0.00136 

643- 

0.206300002 

04 

1 

0.  00136 

644 

0. 21670000B 

04 

1 

0.00136 

645 

0. 21850000E 

04 

2 

0.00272 

647 

0. 220300002 

04 

1 

0.001  36 

648 

0. 225300002 

04 

1 

0. 001 36 

649 

0. 226000002 

04 

1 

0.001 36 

650 

0.  227  100002 

04 

1 

0.00136 

651 

0. 228700002 

04 

1 

0.00 1 36 

652 

0.  23 1 600  002 

04 

1 

0.00  I 16 

653 

0.  232  100002 

04 

1 

0.001  16 

654 

0.242900002 

04 

1 

0.00 136 

655 

0. 243900002 

04 

1 

0.00136 

Figure  20c  - USING  LIST  WITH  TELEPHONE  DATA  SET 

(cent. ) 
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656 

0.2468000 os 

04 

1 

0.00136 

657 

0.24720000E 

04 

1 

0.00136 

658 

0.24900000E 

04 

1 

0.001 36 

659 

0. 26840000E 

04 

1 

0.00136 

660 

0. 27930000E 

04 

1 

0.00136 

661 

0. 28820000E 

04 

1 

0.00136 

662 

0. 28910000E 

04 

1 

0.00136 

‘ 663 

0. 29380000E 

04 

1 

0.001 36 

66U 

0.  308000(JOE 

04 

1 

O.OOl 36 

665 

0. 31170000E 

04 

1 

0.00136 

666 

0. 31b5000CE 

04 

1 

0.00136 

667 

0. 328600008 

04 

1 

0.00136 

668 

0. 33150000E 

04 

1 

0. 00 1 36 

669 

0. 33780000E 

04 

1 

0. 00 1 36 

670 

0. 35100000E 

04 

1 

0. 00136 

671 

0. 38070000E 

04 

1 

0. 00136 

672 

a.  3837000 OE 

04 

1 

0.00 1 36 

673 

0. 385000  OOE 

04 

1 

0.001 36 

674 

0. 38580000E 

04 

1 

0.001 36 

675 

0. 38790000E 

04 

1 

0.00136 

f 

6 76 

0. 38880000E 

04 

1 

0.00136 

! 

677 

0. 39440000S 

04 

1 

0. 00136 

1 

678 

0. 39780000B 

04 

1 

0. 00 1 36 

679 

0. 43730000B 

04 

1 

0.00136 

680 

0.43930000E 

04 

1 

0.001  36 

681 

0. 45120000E 

04 

1 

0.00136 

682 

0. 46370000E 

04 

1 

0.  00 1 36 

683 

0. 47590000E 

04 

1 

0.00136 

684 

0. 4835000CB 

04 

1 

0.  00 1 36 

685 

0.52160000E 

04 

1 

0.00 1 36 

686 

0. 54920000E 

04 

1 

0. 00 1 36 

687 

0- 56520000B 

04 

1 

0.00136 

688 

0. 57540000B 

04 

1 

0. 00 1 36 

689 

0.  6 1850000E 

04 

1 

0.00136 

690 

0.  624  10000E 

04 

1 

0.  00 1 36 

691 

0. 63720000E 

04 

1 

0.00136 

692 

0. 63850000E 

04 

1 

0.00136 

693 

0. 66590000E 

04 

1 

0.00 1 36 

694 

O.S8160000B 

04 

1 

0. 00 1 36 

695 

0. 68210000E 

04 

1 

0.  00 1 36 

0. 73070000S 

04 

1 

0.001 36 

697 

0. 73290000E 

04 

1 

0.00136 

698 

0. 77460000B 

04 

1 

0.00136 

699 

0. 79270000E 

04 

1 

0.00136 

700 

0. 80390000E 

04 

1 

0.00136 

701 

0.  80580000S 

04 

1 

0.  00  1 36 

702 

0. 83530000B 

04 

1 

0.00136 

703 

0. 884700008 

04 

1 

0.00136 

704 

0. 92060000E 

04 

1 

0.00136 

705 

0.  92560000E 

04 

1 

0.00136 

706 

0.95170000E 

04 

1 

0.001 36 

707 

0. 95410000E 

04 

1 

0.001 36 

708 

0. 95620000E 

04 

1 

0.001 36 

709 

0. 9965000CE 

04 

1 

0.00136 

7 10 

0.  1C020000E 

05 

1 

0.001 36 

1 

711 

0.  1037bO00E 

05 

1 

0.00136 

7 12 

0.  1C837000E 

05 

1 

0.00136 

i 

0.  10918000E 

05 

1 

0.00136 

! 

714 

0.  12042000E 

05 

1 

0.001 36 

0.  12793000B 

05 

1 

0.001 36 

716 

0.  13012000E 

05 

1 

0.00136 

717 

0.  1 3179000E 

05 

1 

0.00 136 

718 

0.  146920006 

05 

1 

0.001 36 

719 

0.  1592S000E 

05 

1 

0.00136 

720 

0.  168770008 

05 

1 

0.00136 

721 

0.  18885000E 

05 

1 

0.00136 

722 

0-  19273000B 

05 

1 

0.00136 

723 

0.  196960006 

05 

1 

0.00136 

724 

0.  19849000E 

05 

1 

0.00 1 36 

725 

0.  19851000E 

05 

1 

0.00136 

726 

0. 2C846000E 

05 

1 

0.00136 

in 

0. 21440000E 

05 

1 

0.00136 

728 

0. 2457300QE 

05 

1 

0.00136 

IV* 

0. 24770000E 

of 

1 

0.00136 

730 

0. 262780006 

05 

1 

0.00 1 36 

731 

0.27238000E 

05 

1 

0.00 1 36 

732 

0. 281 J3000E 

05 

1 

0. 00 1 36 

733 

0. 29913000S 

05 

1 

0.00 1 36 

7 34 

0. 39867000E 

05 

1 

0.00 1 36 

735 

0.  5 3 1 38000E 

05 

1 

0.00136 

7 36 

0.  6727100  0B 

05 

1 

0.001 36 

Figure  20d  - 
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Y.  ASSESSMENT  OF  7AE1ABILITY  BCUTIHES 


A.  INIECDQCTION 


In  th€  previous  sections  Prcbatility  Plotting 
Subroutines  were  presented  fcr  use  in  iatornial  estimaticn  of 
the  form  of  the  distribution  of  a set  cf  independent 
observations.  In  the  present  section  subroutines  computing 
the  variability  of  some  basic  statistics  are  described. 

He  know  that  I (sample  mean)  is  an  unbiased  estimator  of 

the  population  mean  ( = x f (x)  dx  ) , where  X is  a 

X 


population  with  unkncwn  distribution  (F  (x)  ) and  unknown 

X 

mean,  and  X , X , X is  a sample  cf  independent 

1 2 n 

Gbservaticns  from  X.  Furthermore  we  can  compute  the 


Variance  cf  J as: 


X 

VarCl]  = VarC.r-—-] 

1=1  n 


1 n 

Var[  ^ X ] 

2 i=  1 i 
n 


n (j 

= — - Yar[X]  = 

2 n 

n 
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Thus  w€  see  it  is  possible  to  estimate  the  7ar[X] 

2 

if  the  population  Variance  ( C ) is  unknown,  by 

2 

usiny  the  sample  variance  S as  follows:  Var[I]  = 
where  S = ^ (X  -X)  / (n-l)  (an  unbiased  estimator  o 

j.-i 


, even 
simply 
S /n. 


f (/)  . 


Therefore  S /n  is  an  unbiased  estimator  of  Var[X]  also. 


Furthermore,  especially  for  a Normal  population,  we  can 
obtain  a confidence  interval  for  the  aean  using  the 
t-statistic . 


£ut,  although  for  X there  exists  a direct  assessment  of 
variability,  for  estimates  of  other  population  parameters 
such  as  Skewness,  Kurtcsis,  Coefficient  of  Variation,  ard  so 
forth,  this  is  not  so  simply  obtained.  Thus  several  methods 
have  been  introduced  to  obtain  assessments  of  variability 
and  for  two  of  them  (Sectioning  of  data  and  the  Jacknife) 
the  Subroutines  'SECTN'  and  'JACK’  will  be  presented. 


E.  SUBBCQIISE  SECTN 


1 . Cescription 

'SECI.V’  Subroutine  is  used  for  assessing  variability 
of  estimates  from  data  based  on  the  Sectioning  Method. 

The  basic  idea  of  this  method  is: 

Assume  we  have  n independent  observations  X , X , 

1 2 

...,  X from  a population  with  unkrcwn  distribution 

function  ? |x).  Let  8 be  a parameter  of  F (x)  and 
X X 

8(n)  = (X  ,X  , . . . , X } be  a statistic  which  estimates  9 . 
1 2 n 
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Sow  we  want  to  estimate  the  variance  of  9(n)  which  is  a 
new  Bandca  Variable  and  we  wcrk  as  follows: 

1.  Divide  the  sample  into  r disjoint  sections,  of 
size  m (r  should  be  such  that  n=!nr.  If  this  is  not  possible 
seme  of  the  latter  data  is  discarded)  . 

2.  For  each  section  form  the  same  estimate  9i (■) » 

i-1^2/.../r 

3.  Compute  the  average  of  9it'»'),that  is; 


- 1 r ^ 

9 (n)  ^ 9 (m)  , 

r 1=  1 1 

doing  so,  we  have  an  estimate  of  9 and 


EC8(n)  ] = EC .2:  9.  (m)  ] 

r i=1  1 


= — rE[8(m)  ] 
r 


= 2C  9(ni)  ] , 


also  we  bave: 


Var[9(n)  ] = Var[ Z 9 (m)  ] 

r i=1  i 


1 r ~ 

VarC  Z 9 (m)  j 

2 i=l  i ^ 


= VarC  9(m)  ] 

r 

1 1 r 


^ ^ ^ 0 ft  2 

.-21,  (O,  (m)  -o{n)  ) 

r r-1  1=1  1 


1 2 
= — s- 
r 8(m) 
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Ihus  we  see  that  Var[  8(a)  ] and  Var[  9 (n)  ] can  be 


unbiasedl;  estimated  from  the  sample  variance  o£  ^-,(0),  the 
main  advantage  and  the  purpose  of  the  sectioning  of  data 
method. 

But  the  main  disadvantage  of  rhis  oethod  is  that  we 
would  like  r to  be  as  large  as  possible  in  order  tc  make 
the  variability  of  the  variance  estimate  S as  small  as 
possible.  Ihis  however  may  be  worse  for  the  bias 
properties  of  the  estimation  procedure.  Therefore  the 
choice  cf  r is  a factor  that  should  be  considered. 

It  should  be  noted  here  that  if  (a)  are 
approxima cel j normal  variates  (This  can  be  tested  by  the 
previous  described  Subroutine  'NORMPL')  then  confidence 
intervals  for  the  unknown  parameter  can  be  obtained  based 
on  the  t-Statistic , in  the  following  way; 


9(n) 


Lhv. 


(1-^/2) , (r-1) 


where : 

9 (n)  is  the  mean  cf  the  sectioned  data  statistics 
obtained  from  the  column  named  'mean'  of  the  second  table 
('estimated  parameters  of  the  sample  parameters')  cf  the 
program  output. 

1/2 

S = /r  IS  the  standard  deviation  cf  the 

9(n) 

sectioned  data  statistics  divided  by  the  square  root  of  the 
number  cf  sections,  obtained  from  the  last  column  of  the 
same  table  of  the  program  output. 

t is  the  (1-  <’^/2)  quantile  cf  the 

(1-  a/2)  , (r-l) 

t-distributicn  with  r-1  degrees  cf  freedom. 
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2,  groaill structure 

'33CI1I*,  using  the  Sectioning  method,  estimates  the 
following  statistics:  flean.  Median,  Variance,  Standard 
Ceviaticn,  Coefficient  of  Variation,  Skewness,  Kurtosis, 
Minimum,  and  Maximum.  (The  formulas  which  have  been  used  to 
compute  these  estimates  are  described  in  the  comments  of  the 
Subroutine.)  The  first  table  is  then  printed  by  the 
program,  containing  the  values  of  these  parameters.  Then 
•SECTN',  using  the  computed  estimates  of  all  sections  for 
each  parameter,  estimates  the  Mean,  Median,  Variance, 
Skewness,  Kurtosis,  and  Standard  Deviation  divided  by  the 
square  root  of  the  number  of  sections.  The  second  table  is 
printed  containing  these  values. 

There  are  three  restrictions  in  using  'SECTN': 

1.  The  number  of  sections  must  be  no  greater  that  100. 
If  it  is,  a diagnostic  message  is  printed  and  only  estimates 
from  unsectioned  data  will  be  given. 

2.  The  number  of  data  values  must  be  greater  than  3; 
otherwise  a diagnostic  message  is  printed  without  any 
calculation . 

3.  The  size  of  each  section  must  be  greater  than  3.  If 
it  is  less  than,  or  equal  tc  3,  then  the  program  gives 
estimates  fcr  the  entire  set  of  (unsectioned)  data.  A 
diagnostic  is  printed. 


A complete  description  of  how  ’SECIN'  operates  is 
given  in  the  Subroutine.  Furthermore,  a Summary  is  given  by 
typing  on  the  terminal  the  command  DESCaiBE  SECTN  under  the 
CMS  environment.  The  following  response  will  be  given  on 
the  terminal  when  the  user  types  the  above  oommand: 


■t 

i 


i 


1 


SOBHCUTINE  SECTN 

•SECTN*  is  intended  to  estiaate  a set  of  tasic 
statistics  c£  a given  set  cf  observaticns  using  the 
•Sectioning  cf  data  method'.  Also  for  each  estimated 
statisric,  estimates  of  some  basic  statistics  such  as  the 
mean,  standard  Deviation,  and  so  forth,  are  given. 

It  is  called  by  ; 

CALL  SECTN  ( X,  N,  K ) 

where: 

X Is  the  array  of  data 

N Is  the  number  cf  data  values  (must  be  greater 
than  3) 

K Is  the  number  of  desired  sections  (no  greater 
than  100) 

Note;  k should  be  a number  which  minimizes  the  number  of 
data  points  that  will  have  to  be  discarded.  'SECTN*  places 
the  data  into  the  equal  size  sections  discarding  any  data 
left  over. 

For  x<3  or  )c>1C0  or  (n/k)<3  only  estimates  from 
unsectioned  data  will  be  given  and  no  estimates  fcr  the 
estimated  statistics  will  be  computed. 

No  output  is  expected  if  n<3  . 

Cere  information  is  given  in  the  subroutine. 


! 
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Osina  * SEC  IN*  with  Islephone  Data  Set  _1 


•SECIN*  was  used  on  Telephone  Data  Set  1 to  assess 
the  variability  in  the  Mean,  Median,  Variance,  Standard 
Deviation,  Coefficient  of  Variation,  Sicewness  and  Kurtcsis. 

The  672  data-points  of  Telephone  Data  Set  1 were 
broken  down  into  16  disjoint  sections  with  42  data-points 
per  section.  Because  of  this  fcreak-down  no  data-pcints  were 
discarded.  Compare  the  values  of  row  'unsectioned'  (see 
figure  21)  with  the  values  of  the  corresponding  statistics 
computed  by  the  Subroutine  ' HISTGS/HISIFS ' (see  figure  1)  . 
The  values  axe  the  same. 

Now  if  we  want  to  assess  the  variability  in  ary  of 
the  above  parameters  we  proceed  as  follows; 

1.  He  take  the  mean  (0  ) of  the  parameter,  whose  the 
variability  we  want  to  assess,  from  the  table  under 
'ESTIMATED  EAHAMETEHS  OF  SAMPLE  PAfiAMETEBS'. 

2.  He  take,  from  the  same  table,  the  value  (S)  from  the 
last  column  (STD.  DEV. /NS**.  5)  for  the  same  parameter. 

3.  Using  the  t-statistic  with  k-1  d.f.  and  the  formula; 

9 ± St 

(1-^V2)  , (k-1) 

we  get  a (1-0')X  confidence  interval  for  the  parameter. 

Example 

In  order  to  find  a 9531  confidence  interval  for  the 
Skewness  of  Telephone  Data  1 we  take; 

8=  4.9343,  3 = .30312,  t = 2.131 

. 975, 15 


It  should  be  noted  here  that  the  use  of  the  Variance 
estimate  frcB  the  sectioned  data  in  order  to  get  confidence 
intervals  of  the  parameters  is  based  on  the  normality 
assumption  and  the  independence  of  the  estimates  from  the 
sections.  The  normality  will  depend  on  the  number  of 
data-points  in  each  section,  which  should  be  Icept  large. 
This  requirement,  however,  conflicts  with  the  need  to  make 
the  number  of  sections  large  to  reduce  the  variability  in 
the  estimate  of  the  variance  of  the  statistics.  The 
skewness  estimates  of  each  of  the  16  sections,  can  be  run 
through  the  normal  plotting  routine  to  see  whether  the  use 
of  the  t-statistic  confidence  interval  is  valid. 
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Figure  21  - OSING  SECTN  «ITH  TELEPHCSE  DATA  SET  1. 
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SQBEOUTISE  JACK 


"•  • C€gcri.ptioD 


This  routine  is  used  for  assessing  variability  of 
estimates  of  data  based  on  the  Jacknife  method,  and  for 
reducing  the  bias  in  estimates  also.  It  is  particularly 
useful  fcr  data  with  small  sample  size. 

A big  picture  cf  the  method  has  as  fellows: 

Lex  i , X , X be  a sample  cf  n independent, 

1 2 n 

identically  distributed  observations  from  a population  with 

unknown  distribution  function  F (x) . Alsc  let  ^ te  an 

X 

unknown  parameter  of  P (x)  to  be  estimated.  Furthermore 

X 

suppose  a method  (biased  or  unbiased)  is  available  for 

estimating  9 . Then  we  proceed  as  follows: 

1.  Divide  the  sample  size  into  r disjoint  groups,  each  of 
size  m (r  should  be  such  that  n = mr.  If  this  is  not 
possible  seme  of  the  last  data-points  will  te  discarded.) 

2.  Compute  the  estimator  8 (n)  of  9 based  on  all  n = mr 

ebservatiens . ^ 

3.  Compute  the  estimator  9 (n-m)  based  on  the  n-m 

i 

th 

ebservaxions , having  deleted  the  i group. 

4.  Compute  the  so  called  'PSEUDO  VALUES'  d (n)  : 

i 

9 (n)  = r 9 (n)  - (r-1)9  (n-m)  , i = 1,2,.  ..,r  . 

i i 
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Then  the  JacJcnife  estimator  is  defined  to  te  the  average  of 
Eseudo  Values: 


— 2 8.  (n) 

r 1=1  1 


---  ( r'8(n)  - 
r 


(r-1)  S’  9 (D-m)  ) 
1=1  1 


9(0)  - 


r 

Z 9 (n-m) 
i=1  i 


Havirg  so  defined  the  Jacknife  estimator  it  can  be 
proved  that  this  estimator  is  an  unbiased  estimator  cf  9, 
except  for  terms  n and  higher  order,  assuming  that  the 

* ^ Cv 

bias  has  tne  form  cf  S[  9(n)  ] = u a/n  + terms  of  higher 

order  in  n.  That  is,  the  Jacknife  estimator  eliminates  a 
-1 

n bias  term.  Namely  if 


EC  § (n)  ] = ^ + an  + 0 (n  ) , 


then  we  have: 

A - £-1  j;  " 

EC9(n)]  = E[r8(n) 3 


rE[8(n)  ] - -----  f E[9.  (n-m)  ] 
r 1=1  1 


Q -j  r- 1 A ra  r 

= r(v+an  •*-0(n  ))~ (rt7-«- * C (r.  )) 


n-x  n-m 


^ ra  Q a 

♦ ....  (r-1)  (V* ♦...) 

n m (r-1) 


ry  ♦ ♦ . 

n 


m (r-1) 
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^ zd  ~ zd  * * d *... 

mi  m 

= a. 

A special  case  of  Jacknife  estisator  cac  be 
■entioaed  here,  called  the  'Complete  Jacknife  Estimate!', 
where  r=a  and  m=1. 

Seme  properties  of  the  Jacknife  estis'  follcw; 

-2 

1.  It  is  unbiased  up  to  order  n 

2.  The  'Jseudo  values'  can  be  used  to  obtain  variance 

estimates  for  the  Jacknife  estimator  since  they  can  be 

considered  (see  [IQ])  as  approximately  independent  and 

identically  distributed.  Thus  (r(r-1))  57  (9  (n)-8(n)) 

i=1  i 

A. 

should  be  an  approximate  estimate  of  Var[9^  (c)  ] , 
and 

"■  1/2 

I_i£i5:ILL 

1~~~{  q"  2 1/2 

{,r  (8  (n)  - y(n))  } 

1=1  1 

should  be  approximatelly  distributed  as  a t-statistic  with 
r-1  d.f.  This  procedure  is  particularly  useful  if  the 
number  n cf  data-points  is  small,  but  it  must  be  used  with 
care. 

3.  ?or  large  n it  can  be  shown  that  fer  the  ccoplete 
Jacknife  estiaatcr,  under  very  general  conditions  we  have: 

A 

Var[  8 (n)  ] > Var[  8 (n)  ] . 
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2 


££SSid£-.§l£uct^e 


Subroutine  'JACK'  is  a FOElSAN-callable  Subroutine 
which  takes  a set  of  data,  groups  it  into  r disjoint  groups 

t h 

and  for  each  set  of  the  r-i  group  coaputes  and  prints  in  a 
table  the  statistics:  Hean,  Median,  Variance,  Standard 
Ceviation,  Coefficient  of  Variation,  Skewness  and  Kurtosis. 
The  same  statistics  are  also  computed  for  the  ungroupec  set 
of  data.  Then,  using  these  statistics,  the  Jacknife 

estimator  of  the  above  parameters  are  computed  and  printed 
in  another  table  along  with  its  Variance  and  Standard 
Deviation . 

The  program  is  divided  into  the  main  program  and  tne 
Subroutine  JACKES.  The  main  program  groups  the  data  and 
successively  calls  the  Subroutine  'JACKES'  in  order  to  get 
the  estimates  of  the  above  parameters  for  ungrouped  daxa  and 
for  each  group  as  well.  A table  is  then  printed  containing 
the  estimated  parameters.  The  'Pseudo-Values'  for  each 
parameter  are  also  computed  by  the  main  program,  and  are 
used  to  compute  the  Jacknife  estimator  for  each  parameter 
and  its  Variance  and  Standard  Deviation. 

To  avoid  division  by  zero  the  number  or  observations 
must  be  greater  than  three,  otherwise  no  output  is  given. 
For  the  same  reason  the  expression  (r-1) (n/r)  must  evaluate 
to  greater  than  three,  otherwise  the  program  will  give  only 
estimates  of  ungrouped  data. 

A complete  description  of  how  'JACK'  operates  is 
given  in  the  Subroutine.  A Summary  is  printed  on  the 
terminal  by  typing  DESCBIBE  JACK  under  the  CHS  environment, 
which  responds: 


99 


SUBBCOTINE  JACK 


Jack  is  intended  to  estimate  the  statistics  : Bean, 
Bedian,  Variance,  Standard  Deviation,  Coefficient  of 
Variation,  Skewness  and  Kurtosis  of  a given  set  of 
independent  chser vaticns,  using  the  Jacknife  method.  In 
addition  to  tae  above  parameters  which  are  given  for  each 
group,  the  Jacknife  estimator  along  with  its  variance  and 
standard  deviation  is  given  for  each  parameter. 


It  is  called  by: 


CAll  JACK 
where: 

X 

xs 

SIAT 

N 

IG 


( X,  XS,  STAT,  N,  IG  ) 

is  the  array  of  data  of  dimension  N 

is  a work  array  of  dimension  N (returns  ordered 

data) 

is  a two-dimensional  work  array,  dimensioned  by 
(1G,7) 

is  the  number  of  data-values 

is  the  number  of  groups  plus  one  (=r+1) 


Note  : r must,  be  such  a number  as  to  minimize  the 
number  of  data-points  that  will  have  to  be  discarded.  JACK 
places  the  data  into  the  equal  size  groups  discarding  any 
data  left  over  (last  observations  in  order  of  original 
input)  . 


No  output  is  expected  if  N < 3.  Furthermore  if 
IG  < 2 or  (IG-2)  (N/(IG-1))  < 3 only  estimates  for  ungrouped 
data  will  be  given. 

Bore  information  is  given  in  the  Subroutine. 
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1 


3 • Dsing  * JACK*  with  Telephone  Da  ta  1, 


Ic  assess  the  variability  in  the  Mean,  Median, 
Variance,  Standard  Deviation,  Coefficient  of  Variation, 
Slcewness  and  Kurtosis  of  Telepfone  Data  1 the  Subroutine 
'JACK*  is  now  used.  In  order  to  avoid  to  discarding  any  of 
the  data-points  we  use  16  groups  with  42  data-points  per 
group . 


Observing  figure  22  we  see  that  the  values  of  the 
paraaeters  for  ungrouped  data  are  exactly  the  same  as  the 
corresponding  values  of  the  same  parameters  of  figure  1 
which  was  produced  by  using  the  Subroutine  HISTGS/HISTFS. 
Now  in  order  to  assess  the  variability  in  any  of  the  abcve 
parameters  we  have  to  use  the  t-statistic  with  (r-1) 
d.f.  along  with  the  values  printed  under 
•ESlItJAlEE  JACKNIFS  PAEAMETEES*  using: 

9 ± St 


, (r-1) 

where  9"  is  the  Jacknife  estimator  and  S is  its  Standard 


Deviation.  For  example  to  assess  the  variability  of  the 
Skewness  with  a confidence  level  .05  we  have: 

9 = 7.37321 , S = .900118, 


t = 2.131  , 

.975,15 


Therefore  a 951  confidence  interval  of  the  Skewness  is  : 
7.37321  ± .900118*2.131  = 7.37321  ± 1.S1815 
*==>  [ 5.4551  , 9.2914  ] . 


This  should  be  compared  with  the  point  estimate  4.9343  and 
the  confidence  interval  estimate  [4.2884,  5.5802]  obtained 
trot  subroutine  SECTN.  The  data  is  so  skewed  that  one  will 
inevitably  have  trouble  here  with  the  Jacknife  procedure. 
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Figure  22  - USING  JACK  WITH  TELEPHCNI  DATA  SET  1. 


102 


KU6TOSS13  6.7079ie  01  H.bdOSbB  OJ  1.71UJ7e  01 


I 

i 

I 


VI.  SINGLE-SEBVEE  FIRST^OME  FIRSI-SiBVED  fiOEOES  SITH 

lE&RMAJ.  SIBUCTUBE 


A.  THE  EA£eA(P,Q)  ERCCESS 


The  starting  point  for  these  processes  is  the  defirition 

of  a first-crder  autcregressi ve  model  for  a stationary 

sequence  of  random  variables  {X  ) : 

i 


X = rx  + 6 ; i = 0,±1 ,±2,  . . . . 

i i-1  i 


If  the  marginal  distrifcution  of  the  X is  fixed  to  be 

i 

exponential  nirh  parameter  7^  for  all  i; 


P{X  < X)  = 1-e  , :^>  0,  x>  0 , 

i 


then  €.  is  zero  with  probability  r and  exponential (>)  with 


probability  "h-r.  Thus 


{rX  w.p.  r 

rX  + E w.p.  1-r, 

i-1  i 


(1) 


where  {E  J is  a sequence  of  i.i.d.  expcner tial  (2^)  random 
variables.  The  process  defined  by  (1)  is  the  exponential 
autoregressive  process  of  order  1,  the  EAR(1)  process.  The 
correlaticn  structure  is  r(j)  = r , j = 0,1,2,...  . 


i 

I 

J 

] 
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The  first  order  moving  average  exponential  process, 
£I1A(1)  is  defined  as 


{‘61.  w . p. 

I + E w.  p.  (1-1 

i i-1 


< ^ < 1;  i = 0,±1,±2,...] 


and  has  ccrreiations  y(1)  = (I-0)  and  p(1)  = 0,  j = 2,3,... 

Ibe  general  EbA  (g)  model  takes  the  form 


/E. 
q i 

i * i E 

g 1 g-1  1-1 


V . p . b 


w.  p . b 


X = 
i 


for  0 


r 2' 


-1  i-1 

* £ E 

1 i-q+1 

w.  p. 

E * . . . 

* £ E 

+ E w.p. 

-1  i-1 

1 i-q+1 

i-q 

..w  / ^ 1; 

q 

' 

i = 0/±2, 

. . . lihere 

i = q+1 f 

0-i)  ... 

q 

1 

1 

q > i > 2 

. (1-  ii')  ... 

(1-4) 

i = 1 . 

The  serial  correlations  are  given  by  the  eguation 


(j)  = corr(X.,X.  ) = 
1 i-r 


q-i+1 
± b b 

V=j  V v+j 

0 


(1  $ j ^ q) 


(g+1  < j <co  ) 
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Tc  ccnvert 


a 


mixed 

ciders 


this  EMA  (q)  process  into 

autoregressive  moving-average  exponential  prccess  of 

1 and  q,  called  Ef!A(1,g)  we  replace  E in  (2)  with 

i-g 


r A 

i-q- 1 

+ E 

i-q 

w.  p.  r 

w.p.  (1-r) 

For  further  results 

on  these 

processes 

and 

their 

properties  see  Gaver  and  Lewis  (1978), 

Lawrance 

and 

Lewis 

(1977),  Jaccbs  and  Lewis 

(1977)  , 

and  Lawrance 

and 

Lewis 

(1978)  . 


E.  USE  CF  AUTOREGRESSIVE  (EAH(P)),  MOVING  AVERAGE  {Ei!A(C)) 
ANC  MIXEI  STEUCTOHES  (EARMA  (P,Q) ) IN  MODELLING  QUEUES 


Consider  for  simplicity  a queue  with  a single  input 

stream  and  a single  server,  and  a f irst-ccme-first-served 

(FIFO)  service  discipline.  Let  S , i = 0,1,2,...  , denote 

i 

th 

the  service  time  for  the  i arrival,  and  let  X , 

i 

th 

i = 1,2,...  , denote  times  between  arrival  of  the  i and 


(i-1)  customers.  As  is  usual  we  assume  that  the  first 


customer  (with  service  time  S ) arrives  at  time  zerc  and 

0 

finds  Che  queue  empty. 


If  the  {S  } and  {X  } sequences  are  i.i.d.  exponential 
i i 

random  variables  with  parameters  ^ and  a respectively,  we 


have  the  M/M/1  queue. 
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Now  let 


0,  + 1 .. 


B 

i 


be  exponential  and  independent,  i 
be  exponential  (cO  an<3  independent,  i 


He  want  to  model  queues  with  correlated  (autocorrelated 
and/or  cross-correlated)  service  and  inter-arrival  times, 
the  service  and  inter-arrival  times  both  having  marginally 
exponential  distributions.  He  also  want  the  queing  model  to 
include  the  a/M/l  queue  as  a special  case. 


There  are  five  simple  possibilities  which  we  put  forward 
here,  based  on  the  use  of  EARMA(p,g)  processes,  giving  what 
we  call  an  EAfiMA  aD/MD/I  queue. 


In  what  follcws  let  p and  p be,  respective! j,  the 

S X 

order  of  the  autoregressive  components  of  the  service  time 

sequence  {S  } and  the  inter-arrival  time  sequence  {X  } , and 
i i 

let  q ,q  be,  respectively,  the  order  of  the  moving  average 
S X 

component  of  {S  } and  {X  ] . These  parameters  can  take 

i i 

values  C,1,...  . If  p = 0,  g = 0,  the  sequence  is 

independent;  if  p = 0,  the  process  is  purely  moving  average 
and  if  q = 0,  then  the  process  is  purely  autcregresi ve . 


The  five  possibilities  or  cases  are  as  fellows: 

1)  Let  (S  } be  EAHaA(p  ,g  ) over  {E  , S , ...}; 
i S S i i-1 

let  iXJ  = {£.}. 

1 1 

Thus  the  arrivals  are  a Poisson  process  and  the  service 
times  are  autocorrelated,  i.e.  a mixed  autoregressive 
moving  average  exponential  sequence. 
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2)  Let  S = E , i = 0,1,2,  ...  ; 
i i 

let  X te  EA2aA(p  ,q  ) over  [6  t 6 • •••}• 

i XX  i i-1 

Then  the  service  times  are  independent  and  the  arrival 
process  is  non-Poisson  because  of  the  dependency  between 
the  inter-arrival  times,  i.e.  the  arrival  process  is  a 
point  process  with  EAHI!A(p,g)  structures. 


Let 

{S.} 

1 

be 

EAB£lA(p  ,q  ) 

S S 

over 

{£./ 

1 

E , . • . } • 

i-  1 

Let 

{X.} 

1 

be 

EASSA  (p  , q ) 

X X 

over 

1 

5,  ...}. 

1-  1 

The 

service 

and  arrival 

processes 

are  autocorrelated. 

but 

the 

two  processes 

are 

independent.  The  marginal 

distributions  of  {S  } and  {X  } are  still  exponential. 

i i 


4)  To  couple  the  two  processes,  with  resultant  dependence 
in  the  arrival  and  service  processes,  the  sinplest 
procedure  seems  to  be  the  following. 

Let  JS  } be  EARilA(p  ,q  ) in  the  following  sense: 
i S S 

{S.j  is  ZAHMA(p  ,q  ) over  {E  . , » , a 5;.  , ...}  . 

1 SS 

Then  if  X = & , i = 0,1,2,...  we  have  that  {S  } is  an 
i i i 

autocorr elated  sequence  and  also  cross-ccrrelated  with 

{X  = & },  i.e.  £S  ,X  } is  a bivariate  dependent 
i i i i 

sequence  cf  random  variables  with  exponential  marginal 

distributions,  although  {X  } is  still  Poisson.  Note 

i 

that  the  S sequence  is  an  autccor related  sequence 
i 

because  of  the  cross-coupling,  but  not  a pure 
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EARMAip  ,q  ) sequence.  It  does  have  an  exponential 
S S 

marginal  distribution  however. 


As 

an 

example 

let  p = 
S 

1, 

= 1. 

S 

0 

= 

E 

0 

S 

1 

= 

£ E 

S 1 

V .p. 

A = a & 

1 > 1 

= 

♦ A 

1 

w .p . 

) 

s 

s 

2 

w .p . 

4 

A = rA 

2 1 

w.p 

= 

♦ A 

2 

w .p. 

(1-/ ) 
s 

= rA  + 

1 A 2 

w.p 

S = ^ E 
1 Si 


'. p . £ 


A = rA 
i i-1 


w.  p.  r 


= ''•P-  J = ''•P* 

^Si  1 S i-1.;^i 


In  these  equations  one  can  replace  6 with  X and  then 

i i 

the  cross-coupling  between  the  sequences  {S  } and  {X  } 

i i 

becomes  apparent. 

Interpretation.  We  have  positive  correlation  between  S 

i 

and  q previous  inter-arrival  times  (when  p =0).  If 
S S 
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the  6^  s (j  = i-Pg"!) 


are  short,  thee  S 


will  te  short  if  all  the  g previous  inter-arrival  times 

S 

are  short.  This  models  the  case  where  the  server  tends 
^2  Ji£  queue  qets  lonq.  Cf  cource  he  also 

slows  down  when  the  queue  gets  shert  and  it  is  not 
immediately  clear  what  the  effect  on  an  average  waiting 
time  will  he. 


4*)  Similar  to  4 but  {X  } is  EABllA(p  ,g  ) over 

i XX 

^ , IE  ,...  andS  ~E. 

i o.  i cx  i-1  i i 

Interpretation  is  net  clear,  but  it  would  have  the 

same  effect  as  balking  in  the  input  stream.  If  service 

times  get  long  (and  presumably  the  queue  gets  long)  then 

inter-arrival  times  get  long. 

5)  A mere  general  possibility  which  allows  one  to  model 
dependency  iu  tha  input  stream,  in  the  service  stream 
and  dependency  which  couples  them  is  as  fellows: 

Let  {£  } be  a unit  exponential  indepencent  seguerce,  a 
i 


driving 

sequence. 

Let 

a.) 

1 

be 

EA311A(p  ,q  } 
S S 

over  {E  , E ,. 

i i-1 

Let 

a.) 

1 

be 

EABSA  (p^, q^} 

ever  (<&. , e. 

1 1-1 

..}  . 

Let 

is.) 

1 

be 

SABHA(p  ,q 

C,  C 

) over  {E  , , IE  . , 

^Vi 

Let 

<X.} 

1 

be 

EABMA<p  ,q 

C^  C 

) over  {I  , ckc  , 
a i i 

OD 

i-1 
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This  is  a very  gaieral  scheme  which  reduces  tc  all  of 


the  previcus  schemes  as  special  cases.  Fcr  example  if 


F =p  =g  =0,  then  we  have  case  3.  If 

Cl  Cl  C^  Cj^ 

g =1,  p =0,  / =0,  p=g=0;p=0  (and  p 

C35_  Cj  Cj_  s S I Cl 

and  g are  not  equal  to  zero)  we  have  case  4. 

Cl 


In  all  cases  the  basic  equation  fcr  M/H/l  queue  waiting 

time  H still  holds: 
n 


a = 0 
0 


a =(a  +£  —x  ) n=i^2f...  . 

n+1  n n n+1 


These  equations  can  be  used  to  generate  successive  i s 

n 


in  a simulation.  We  have  not  touched  on  hew  the  correlated 
sequences  are  started.  This  is  somewhat  arbitrary  as  it  is 
not  Known  (except  for  the  Eai(1)  and  I3ii(l)  processes)  hew  to 
start  them  so  as  tc  produce  a stationary  expcnential 
sequence.  The  problem  is  aggravated,  as  in  cases  4,  4 and 


5,  when  we  are  dealing  with  bivariate  exponertial  sequences. 
Then  the  aarginal  processes  may  be  stationary  but  net  the 
bivariate  process. 


C.  PfiOGBAM  SIBOCTOEE 


1.  G6D€ral 


Tc  sifflulate  the  £AfiHA(p,g}  model,  a BORTRAti  Pic 
has  been  written  tc  take  care  of  all  the  existing  cases 
that  model  and  also  of  the  a/M/1  queue.  Eepending  cn 
input  values  which  the  parameters  (as  described  below) 
take,  the  program  can  be  used  for  creatirg  in  each  r 
particular  case  of  the  EARaA(F,q)  queue.  Eecause  of 
generality  cf  the  program  it  is  not  suggested  that  i 
used  for  those  EARMA  cases  where  either  the  arrival  tire 
the  service  times  cr  both  are  not  autccorrelat ed 
cross-correlated.  The  reason  for  this  suggestion 
efficiency.  It  is  suggested  for  use  with  the  EAEI1£( 
cases  where  both  arrival  and  service  times 
autocorrelated  and,  furthermore,  cross-ccrrelated  with 
crder.  The  program  is  divided  into  the  main  progran 
into  the  Subroutines  BETAS,  AUTOR  and  EARdA. 


2*  Main  Program 

The  main  program  performs  all  I/C  operations 
calculations  required  for  the  desired  statist 
Specifically  it  is  designed: 
a)  To  read 

1.  All  random  number  generator  seeds  require 
generate  the  sequences  of  the  Exponential  and  Uni 
variates  to  be  used  during  rhe  execution. 

2.  The  values  for  the  variables: 

N : Humber  of  arrivals  to  be  generated. 


gram 
of 
the 
will 
un  a 
that 
t be 
E or 
or 
is 

P/4) 

are 

any 

and 


and 

ics. 

d to 
form 
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w 


a : Sumber  of  replications. 

ax,  as  : Arrival  rate,  service  rare. 

CfX,  CES,  CPG  : Taking  the  values  1 cr  0 (1  fcr  the 

coupled  process) . 

KX1,  tHj  KS1,  KS  : The  order  of  movirg  average  for 

arrival  and  service  processes. 

QX1,  CX,  QS1,  QS  ; Taking  the  values  1 or  0 (1  for 

autoregressive  case) . 

£X,  BX1,  BS,  BS1  : Arrays  with  the  values  of  for 

moving  average  part. 

aflOXI,  aflox,  BHCSI,  BHCS  : Parameter  values  for  the 

autoregressive  parts. 

t)  To  calculate  and  print 


SU13X  (I,J)  ; 

SUES  (I,J)  ; 

H U,J)  ; 

«E  (I,J) ; 

D <I/J)  ; 

fcr  I = 1,2, 3, 4, 
where : 

(In  what  follows  k gets  the  values  N/4 
fcr  1=1, 2, 3, 4 respectively  and  J=1,2,. 


SOMXIl  (I  ,J)  ; 
SDMSH  (I,J)  ; 

HM  (I,J)  ; 

WME  {I,J)  ; 

DM {I,J)  ; 

J = 1,2, ...,M, 


M/2,  3N/4,  M 

rM) 


k (J)  (J)  th 

SUMX(I,J)  = 2:  X , X the  r arrival  for 

r=1  r r 


th 

the  J replication  in  tae  EAEUA  queue; 

k (J)  (J)  th 

SUaXM(I,J)  = Z (XM)  , (XM)  the  r arrival 
r=  1 r r 

th 

for  the  J replication  in  M/M/l  Cueue; 


k (J)  (J)  th 

soas  |I,J)  = S , S the  r service  for 
r=1  r r 

th 

the  J replication  in  the  EABMA  queue; 
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Ic  {J)  (J)  th 

SDI!SI!(I,J)  = z (SM)  , (SM)  the  r service 
r=  1 r r 

th 

for  the  J replication  in  M/M/1  Queue; 


(J)  (J)  (J)  (J) 

= W = max{  (H  ^ + S , - X 0) 

Ic  k-1  k-1  k 

th  th 

= waiting  time  of  k arrival  in  J replicaticn 
of  the  correlated  queue; 


«M(I,J)  = (WM) 


(J) 


The 


same  as 


but 


for  M/M/1  queue; 


k (J)  ^ 

»B (I,J)  = 37  W /k 
r=1  r 

J.  (J)  , 

aM£  (I,J)  = z (WM)  /k 
r=1  r 

D (I,J)  = {SUMS  (I,  J)  - SaMX(I,J)}/K  - 1/RS  + 1/FX 
DM(I,J)  = {SDMSM(I,J)  - SUMXM  (I,  J)  )/K  - l/SS  ♦1/HX 
k (J) 

S0flSM(I,J)  = Z (SM)  , 

r=1  r 


Mete:  Between  the  group- variables  (KX1,  KX,  KS1,  KS) , (CXI, 
QX,  QS1,  QS) , (CPX,  CPS,  CPG) , (HHOXl,  HBCX,  HHOS1,  EHCS) 
and  (BXl,  EX,  BS1,  3S)  there  are  the  fcllowing  relations: 

1.  The  number  of  elements  of  each  variable  (array)  from 
B-group  should  be  the  same  as  the  value  cf  the  corresponding 
variable  frea  K-group. 

2.  The  *0'  value  cf  a variable  from  Q-greup  implies  that 
the  corresponding  variable  from  H-group  is  net  needed. 

3.  The  value  '1'  of  CPG  dominates  any  value  of  CPX,  or 
CPS. 

Any  valid  combination  in  the  values  of  the  above 
variables  specially  cf  the  Q-group,  C-greup  and  K-greup  (0 
cr  not  0),  generates  a particular  EAfiMA(p,g)  case  and  is 
accepted  by  the  program. 
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Ibe  basic  approach  iu  the  main  program  is  to 
generate  at  once  and  for  each  replication  three  independent 
exponential  sequences  with  parameters  XL  (arrival  rate),  SI 
(service  rate)  and  1 (unit  exponential)  which  are  stored  in 
the  arrays  EXPXL,  EXFSI  and  EXP1  respectively.  Ther  the 
Subroutine  AQTOB  which  is  involved  with  the  autoregressive 
part  of  the  model  is  called  (if  it  is  needed  for  the 
particular  EABflA  case)  . A Icop  follows  which  is  executed  as 
many  times  as  the  value  of  N.  From  the  loop  the  EAEHA 
Subroutine  is  called  which  is  involved  with  the  scving 
average  part  of  the  model.  The  number  of  calls  of  Subroutine 
EABMA  depends  on  the  particular  case  which  is  being  rur.  The 
statistics  and  all  infcrmaticn  that  will  help  us  to  analyze 
the  model  are  also  computed  inside  the  Iccp  and  stored  in 
the  arrays  described  above.  The  program  continues  execution 
until  the  desired  statistics  for  all  M replications  have 
been  calculated  and  gathered.  Then  the  output  part  of  the 
program  fellows  which  gives  us  the  values  of  calculated 
statistics  on  paper  and  on  punched  cards  for  analysis. 

3*  S ubreutine  BETAS 

This  is  a simple  subroutine  and  its  purpose  is  to 
return  an  array  SUMBX1/X/S1/S  with  element-values  as 
follows: 

I 

SUflBX(I)  = X"  EX(r) , I = 1 , 2, . . . , value  of  K-group 
r*1 


to  be  used  for  the  choice  of  the  order  in  the  moving 
average.  It  is  called  at  most  4 times  during  the  execution 
and  only  if  the  K-group  variables  have  value  greater  than 

•O'. 
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Subgcutine  AOTCR 

'AOICB*  accepts  in  each  call  an  expcnential  secuence 
(generated  previously)  of  variates  and  transforms  it  into  an 
autocorr elated  sequence.  It  is  called  only  if  any  of  the 
Q-group  variables  equals  1.  For  each  Q-group  variable  having 
value  1 it  is  called  a times. 

5*  Subrcutine  SABKA 

This  is  the  main  Subrcutine  of  the  program  and  it 


has  been  created  so  that  it  can  be  used  for  all  £ABHA(p,q) 
cases.  Because  of  the  generality  quite  a feu  parameters  are 


transferred  and 

in 

each  call 

just  one  value  is  returned. 

I 

This  value  is  the 

th 

I 

arrival  or 

th 

the  I service  time  as  it 

> 

i 

has  been  modified 

because  of 

the  loving  average.  It  is 

1 

-» 

called  only  if 

the 

particular 

£ABhA(p,q)  case  requires 

i 

moving  average;  that  is  only  if  any  of  the  K-group  variables 
is  greater  than  0. 

^ requirements 

Seme  CPD-times  have  been  gathered  in  the  course  of 
the  simulation  of  some  EABtlA(p,q)  cases,  using  the  IBM 
360/67,  FCRTBAN  H compiler  as  follows: 


n 

K-GBOUP 

QSI 

C-G300P 

TiaE 

2000 

500 

0 

1 

0 

9 

Bin. 

51 

Sec . 

5000 

500 

0 

1 

0 

22 

Bin. 

25 

Sec . 

10000 

500 

0 

1 

0 

43 

Bin. 

30 

Sec . 

1000C 

500 

0 

1 

CPS=1 

55 

Bin. 

20 

Sec . 

from  these  results  we  see  that  the  value  cf  M is  the  main 
factor  that  affects  the  CPO-time  (linearly) . He  see  also 
that  one  cf  the  C-group  variables  aay  increase  the  CPD-time 
by  2531.  He  ignore  the  factor  of  M (it  is  another  sain 
factor  as  N is)  because  we  may  step  the  prograit  at  any 
replication,  getting  the  last  random  number  generator  seeds 
and  then  continuing  another  time  using  these  seeds. 

D.  AiNALXSIS  OF  SIMULATION  FESDLTS 


Since  nc  analytical  properties  of  the  EAHMA(p,g)  models 
can  be  derived,  a simulation  has  been  done  to  study  their 
properties.  Eut  as  we  have  noticed,  this  program  is  a 
general  program  accomodating  all  EASMA(p,g)  cases.  Since 
our  study  is  limited  to  only  the  two  cases  below,  in  erder 
to  save  CEO-time  a modified  program  has  been  usee  that 
simulates  these  particular  two  cases. 

• iQiigSI  with  Dependent  Service 

The  mcdel: 


Ihat 


.95, 


(X.} 

1 


= e. 


is.} 

is  EABMA(1,0)  over  { E 

...  } 

1 

i 

1-1 

is; 

r r A w . p.  r 

1 i-1 

where 

A = E 

[ rA  * E w.  p.  (1-r) 

i-1  i 

To  simulate  this  model  we  chose  the 

values 

0 0 

.25,  .50, 

.99  for  traffic  intensity  t and  the 

values 

.25,  .50, 

I 

I 

I 

( 
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.90,  .95,  .96  for  the  correlation  r,  in  order  to  co^er  a 
representative  range.  In  what  follows  the  string  S#tt#rr 
stands  for  the  run  with  t=tt  and  r=rr . That  is  the  run 
Si25#50  stands  for  the  run  with  t=.25  and  r=.50  For  each  t 
and  r we  ran  the  program  and  the  required  statistics  were 
gathered,  analyzed  and  plotted  using  the  Subroutines 
HISTGS/FS,  NCRdPL  and  EXPLT.  The  sample  statistics  from 
HISTGS/FS  were  also  tabulated  separarely  frca  the  figures  so 
that  one  could  obtain  an  overall  picture  of  whether  or  not 
the  waiting  times  had  converged  in  mean  and  in  distribution 
to  tne  limiting  distribution. 

From  the  analysis  of  all  runs  and  plots  we  have  the 
fcllcwing  results: 

1.  Simulating  the  type  S#25#rr  and  S#5Cirr  model  (with 
r =.25,  .50  and  .90)  it  was  possible  for  the  W and  i to 
reach  the  steady  state  for  U=2000.  On  the  ether  hand  for 
the  type  S#25irr  and  S#50#rr  (with  rr  = .95,  .98)  it  was 
necessary  to  go  up  to  11=10000,  in  order  for  the  M and  U to 
reach  the  steady  state.  Thus  we  see  that  the  high 
ccrrelaticn  affects  the  choice  cf  the  value  cf  N,  requiring 
N to  increase  as  the  correlation  is  increasing.  But  not 
only  the  correlation  affects  the  choice  of  N.  The  traffic 
intensity  affects  it  much  more,  since  for  the  type  S«95#rr 
and  S#99#rr  the  W and  U do  not  reach  the  steady  state  for 
N=10000  even  if  r is  low.  Because  cf  the  CFD-time 
requirements  we  restricted  ourselves  to  a detailed  study  of 
the  particular  case  S#99#98  increasing  successively  the 

value  of  N up  to  320,000,  where  the  W and  U appear  to  be 

close  to  the  steady  state.  Thus  ve  may  ccnclude  that  high 
traffic  intensity  and/cr  high  correlation  require  a large 
value  of  N in  order  to  achieve  steady  state.  Figures  23a 
through  23d,  24a-24c,  25  and  26a-26d  justify  the  abeve 
conclusions.  In  figures  23a-23d,  where  we  are  dealing  with 


I 


i 

■ I 


I 

i 


j 


\ 

1 


r 


117 


th 

the  waiting  time  of  the  2000  (»  ) arrival  for  the  case 

2000 

t=0.50,  r=0.25,  we  can  see  that  H has  converged  to  a 

2000 

value  of  about  6.3/  and  furthermore  a straight  line  appears 

in  the  plot  under  EXPLT.  He  also  see  approximate 

convergence  for  the  correlated  queue  for  the  i , which  is 

2000 

presented  cn  figures  28a-28c,  and  its  value  has  converged 
to  a value  of  about  3.03  (figure  23d) . Note  that  the 
introduction  of  correlation  into  the  service  time  has 
increased  the  average  waiting  time  from  2.5  tc  about  3.03, 
cr  by  approximately  20%. 

He  cannot  say  the  same  for  figures  24a-24c  where, 
because  cf  .the  high  correlation  (r=0.98)  we  do  not  have 
convergence  of  H and  TJ,  although  N=10,000,  and  for  figure 

25  where,  because  cf  the  high  traffic  intensity  (t=C.S9)  we 
do  not  have  convergence  of  H and  ¥.  Figures  26a-26d  also 
show  us  the  non-convergence  of  W,  ¥,  because  of  the  high 
correlaticn  and  the  high  traffic  intensity  (t=0.99, 
r=0.98).  This  is  the  case  for  which  we  eventually  had  to 
carry  the  siaulaticn  out  to  about  N=500,0C0  to  observe  a 
stationary  queue. 
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Figure  23a  - QUEOE  WITH  E&R1  AUTOREGRESSIVE  SERVICE  TIME 

SECUEHCE  AND  POISSON  INPUT.  HISTOGRAM  OF  THE  WAITING  TIMES 

W WITHOUT  ZEROS  PROM  THE  RUN  S#50#25;  in=500 

20  00 

REPLICATIONS,  RX=.2,  RS=.  4.  No  OF  ZEROS  = 255. 


119 


NUBBEB  Of  URDCBED  fAlBS  • 24B 


EsiiBAtto  i-AHAnnEio  ue  data  : oAnnAi  - f.iiiTnie  oo  oambaj  • ■t.&aasooi  oo 


Figure  23c  - QUEUE  WITH  EAS1  A OTCREG BESSI V 1 SERVICE  TIME 
SEQUENCE  ANt  POISSON  INPUT.  NORMAL  PLOT  (NORMPL)  CF  THE 
WAITING  TIMES  H WITHOUT  ZEROS  FROM  THE  RUN  S#50#25; 


WAITING  TIMES  H WITHOUT  ZEROS 

200  0 

id  = 500  replications,  EX=.2,  RS=.4. 


No  OF  ZEECS=255. 
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CCfiReLkTEIl  CUCUE 


figuie  23d  - flUEUE  EIIH  EEB1  EUTCSEGBESSIE  E SEHVICE  I1H2S  itID  ECISSCN  mPUI.  TEBUIETICH  Ot  S»(flE  STETISTICS  lOR  ThE 
CISlBlEUTIOEi  OE  CUBUIAIEO  INIEEEBBl  VAL  TldES  AT  B-^ldOO  AND  N>:2a00  (X  2 ANC  I U),  CUHULATED  SERVICE  TIMES  (S  2 AMD  S «)  , 
NAIJISG  IISES.  HTH  CS  BlIHOOT  ZSRCS  {MI  » 0,  HI,  I FOB  CASE  N-500,  lOOO.  1500,  2000),  CI'BDLAIEE  NAIlIliG  TIMES  (N1  SAB  etc.)  AND 
THE  AVEBAGEC  CIEEEBEDCES  BEiUEEN  S2  AND  X2  (02);  OtIAINEC  EBON  II^SCO  BEEIICAIIGKS  CE  THE  BUN  SIS0I25;  EX=.2:  BS-.H  . 
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igure  24a  - QUEUE  WITH  SAP 
EQUENCE  ANE  POISSON  INPUT. 

WITHOUT  ZEROS  FROM 

10000 

EPLICATIONS,  RX=2.5,  HS=10; 


1 AUTOREGRESSIVE  SERVICE 
HISTOGRAM  OF  THE  WAITING 
THE  RUN  S#25#98; 

No  OF  ZEROS=354. 


TIME 

TIMES 

in=500 
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Figure  24t  - QDEOE  BIIH  EAR1  AOTOSEGBESS I VE  SERVICE  TIME 

SEQUENCE  AND  POISSON  INPUT.  EXPONENTIAL  PLCI  (EXPLT)  CF  THE 

BAITING  TIMES  N WITHOUT  ZEROS  FROM  TEE  RUN  S#25#98; 

1 0000 

ffl=500  REPLICATIONS,  RX=2.5,  ES=10;  No  OF  ZES0S=354. 
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ElilirtTiU  rtmiltlEilS  or  UAT*  : CANRII  . &.U<ll2SU(  OO  CABRA2  • i.21fc«nc  01 


MUABBI  OP  OBDEBBO  FAIBS  • 146 


I 


Figure  24c  - QUEUE  WITH  EAH1  AUTOBSGfiESSIVE  SERVICE  TliJS 

SEQUENCE  AND  POISSON  INPUT.  NOSilAL  PLOT  (NOP.MPL)  CF  THE 

WAITING  TIMES  M WITHOUT  ZEROS  FROM  THE  RUN  S#25#98; 

100  00 

iB  = 500  REPLICATIONS,  RX=2.5,  RS=10;  No  OF  ZEROS  = 354. 
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figure  24(1  - QUEUZ  MIIH  E1R1  AUTCREGKESSIVE  SERVICE  USES  ANC  ECISSCE  INEUI.  TAbDlATICN  Of  SAEIIE  STATISTICS  fOB  THE 
ElSIfllUllCIIS  Of  COeUlAlED  I NT  £B  ABEI V A L IIHES  AT  M>5000  AND  N>  10000  (X  2 ANC  A 4),  CUnULATEO  SERVICE  TINES  (S  2 AKU  S 4|  , 
hAlTlNC  TINES.  AITH  CB  NIIHOUI  ZERCS  (MI  * 0,  HI.  I FOR  CASE  N^2S0C,  SOOO,  1500,  10000).  CUMULATED  hAIIIKG  TIMES  (N1  EAR  etc.) 
AND  THE  AVERAGED  DlFfEREACES  BETUEEN  S2  AND  12  (D2)i  OEIAIUED  FUCH  l-SOO  BffllCATICNS  CF  THE  RUN  SI2SA9E;  BI«2.S;  BS<10  . 
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Figure  24e  - QUEUE  WITH  EAR1  AUTOREGRESSI V£  SERVICE  TIMES 

AND  POISSON  INPUT.  HISTOGRAM  OF  THE  CUMULATED  AND  AVERAGED 

WAITING  TIMES  ¥ FROM  THE  RUN  S#25#96;  iii  = 500 

1 00  00 

REPLICATIONS,  aX=2.5,  ES= 1 0 . 
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figure  24f  - QUEUE  »1TH  EAH1  AUTOHEGHESSI VE  SEfiVICE  TIHES 

ASD  POISSON  INPUT.  EXPONENTIAL  PLOT  (EXPLT)  CP  THE 

COMOLATEE  ANE  AVEHAGED  WAITING  TIHES  Tl  FROH  THE  BUN 

10  0 0 0 


S»25#98;  fli=500  BEPLICAIIONS,  BX=2.5,  HS=1Q 


tS'JIDJtllO  PmiHITEBS  UF  UATA  : GAHHAI  • I.OIUIIIK  00  6ABBA2  • 0.62ll«ill  00 


HUHBKM  Cl  OBDENCO  IllBS  > bOO 


Figure  24g  - 
AND  POISSON 


QUEUE  WITH  EAF1  AOTOfiEGBESSI V£  SERVICE  HUES 
INPUT.  NORMAI  PLOT  (NOBflPL)  CF  THE  CUMULATED 


AND  AVEBAGEC  WAITING  TIMES 

100  00 

1=500  REPLICATIONS,  RX=2.5,  RS=10  . 


FROM  THE  RON  Si25#98; 
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CCKKELATED  CUEUE 


ElQuie  23  - QOEUE  DlIU  E1B1  EOIOEE6BESS1VE  SEBVICE  IIBES  ENC  ECISSCN  IMEUl.  lEEUlETICB  Of  SlEflE  STtllSTICS  FOB  THE 
ElSIfif OllCBS  Of  CVeULiltD  INIEEEBBIVil  TltlES  El  M>S000  END  1000U  (1  2 EMC  1 II),  CUBULEXEO  SEBVICE  TIBES  (S  2 EBO  S II), 
BEIUBG  USES,  ilTH  CB  BlIHOUT  ZEBCS  (HI  ♦ 0,  BE,  1 EOB  CESE  11-2500,  5000,  1500,  lOCCO)  , CUBULEIED  hEIimc  HUES  (SI  BEB  etc,) 
EBB  lUf  iVEBEGEO  OlfffBEfCES  BEIHEEII  S2  END  12  (G2)  ; OE1EJBEO  FBCB  1-500  BEEllCEllCBE  CF  THE  BUN  5199125:  FX-2.97;  BS-J  . 
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Figure  26a  - 
AND  POISSON 


QUSOS  WIT1  EAH1  AUTOREGRESSIVE  SERVICE  TIMES 


INPUT. 


OF  THE  WAITING  TIMES  tt 


HISTOGRAM 

1 0000 

WITHOUT  ZEROS  FROM  THE  St99#S8;  in  = 500  REPLICATIONS,  KX  = 2.97; 


SS  = 3 


No  OE  ZER0S=19, 
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NOHMt-H  UK  LIkDi-KKO  fAlMS  » 401 


t * 


Figure  26b  - QUEUE  WITH  EAR1  AUTOREGRESSIVE  SERVICE  TIMES 

AND  POISSON  INPUT.  EXPONENTIAL  PLOT  (EXPLI)  OF  THE  WAITING 

TIMES  W WITHOUT  ZEROS  FROM  THE  RUN  S#99#98;  in=500 

1 000  0 

REPLICATIONS,  RX=2.97;  RS=3  ; No  OF  ZER0S=19. 


EjUKATtl)  PAKAMbTtK^  ulK  DATA  ; CAHMAl  » 1.25633CE  00  UAHrA2  « 1.43e)77C  00 


Figure  26c  - QUEUE  WITH  EABI  AUTOREGRESSIVE  SERVICE  TIHES 

AMD  POISSON  INPUT.  NORHAI  PLOT  (NORMPL)  OF  THE  WilTING 

TIMES  W WITHOUT  ZEROS  FROM  THE  RUN  S#99#93;  in  = 500 

1 0000 

REPLICATIONS,  RX=2.97;  RS=3  ; No  OF  ZEH0S=19. 


133 


CUMMEIATEO  CUEUt 


1 

O 

1 

MB 

A 

9 

<•« 

O 

i 

I 

MB 

.3 

A 

9 

9 

A 

.3 

Mi 

«( 

s 

O 

o 

9 

4 

1 ^ 

Mi 

i> 

•B 

O 

3 

A 

Mi 

N 

■f 

o 

JB 

Q 

A 

A 

3 

O 

j -1 

« 

Mi 

MB 

M 

A 

' 

o 

1 

.X 

Ai 

M 

• 

■N 

■«#> 

A 

i 

O 

iv 

A 

A 

3 

A 

9 

o 

* 

! o 

•«B 

O 

V» 

Bto 

o 

Mi 

A 

O' 

(3 

UB 

u» 

Ito 

h- 

A 

3 

♦ 

•a 

■4) 

♦ 

9 

-r 

-3 

A 

3 

o 

Mi 

o 

' 

•O 

A 

O 

1 'T 

« 

• 

'V’ 

9B 

MB 

o 

MB 

M 

A 

o 

t 

o 

o 

X 

Ml 

MJ 

o 

! 

j 

4V 

j 

<v 

i> 

1 

A 

U 

o 

’.J 

iB 

o 

i 

lA 

r- 

A 

M. 

.3 

a 

• 

1 

it 

»• 

(M 

Ci 

A 

9 

a 

1 

«A 

tf 

Ml 

o 

A 

«» 

Ml 

a 

o 

MB 

W' 

MB 

ii 

Mi 

o 

A 

M» 

A 

a 

o 

o 

Mi 

m 

Vi 

j 

MJ 

A 

A 

MB 

a; 

'J 

t 

• 

• 

» 

• 

V* 

bT 

Mi 

MJ 

Vi 

3 

MB 

1 

4> 

M»- 

Mi 

Mi 

A 

— 

■T 

1 

A 

r» 

MB 

w 

M* 

O 

A 

O 

MB 

A 

MB 

a 

O 

A 

9 

9 

3 

M* 

Ml 

IT* 

9 

3 

o 

.3 

o 

Mi 

9 

MB 

3 

3 

A 

9 

A 

X\ 

o 

M- 

m 

IfB 

3 

•O 

A 

V 

r- 

o 

O 

A 

o 

1 ♦ 

■o 

o 

'9 

■T 

o 

3 

♦ 

•B 

A 

3 

A 

M 

M- 

A 

MJ 

s 

MB 

9 

O 

O' 

A 

A 

MB 

CO 

O 

3 

c 

A 

9 

I MB 

/M 

n 

Ml 

(A 

Ml 

3 

e 

3 

c 

a 

O 

O 

o 

MJ 

3 

O 

3 

O 

3 

c 

o 

Ml 

3b 

X 

Mi 

MJ 

3 

o 

o 

3 

A 

a 

o 

O 

i 

} 

3 

OI 

o 

<u 

A 

A 

h- 

«J 

3 

o 

A 

Mi 

«• 

9 

9 

3 

o 

U 

Ci 

«u 

U 

A 

O 

> 

Mi 

VB 

ii 

o 

MJ 

Mi 

<0 

O 

A 

l>» 

3 

Mi 

T 

o 

A 

A 

MB 

a 

Mm 

•*> 

3 

X 

A 

9 

9 

1 

MB 

MB 

<0 

<0 

«> 

3 

A 

O 

9 

Ml 

m 

A 

a 

9 

O' 

tfB 

MB 

M> 

<0 

A 

3 

A 

J> 

A 

A 

Mm 

9 

tt 

o 

O 

A 

rg 

A 

A 

9 

M» 

3 

3 

3 

9 

A 

Mi 

o 

O 

.f 

'O 

MJ 

3 

A 

9 

Ml 

O 

c 

3 

a 

a 

ai 

MB 

s 

•« 

o 

U 

•M 

3 

o 

A 

9 

A 

9 

3 

3 

C 

o 

Ml 

1 

JJ 

MJ 

MB 

3 

1 ♦ 

O 

o 

O 

3 

♦ 

fj 

A 

A 

A 

O 

a 

a 

a 

3 

4» 

O 

VB 

o 

1 

1 

Vi 

9 

!»• 

A 

A 

MB 

a 

A 

ilJ 

fH 

A| 

o 

3 

9 

9 

■J 

3 

s 

JJ 

X 

3 

M 

1 

3 

1 

O 

at 

M» 

«>i 

Vi 

<1 

VM 

•r 

UJ 

3 

9 

JB 

3 

Vi 

MJ 

a 

lit 

A 

<o 

oc 

© 

Mj 

9- 

MJ 

-» 

a 

A 

Ml 

MB 

■r 

o 

A 

cc 

A 

A 

o- 

r- 

Or 

< 

Mi 

c 

O 

A 

3 

A 

«J' 

A 

o 

3 

«* 

A 

A 

3 

Mi 

A 

3‘ 

« 

o 

•o 

rt 

<0 

.M 

3 

*- 

MB 

9 

3 

A 

3 

9 

M 

A 

A 

A 

r>« 

(T 

7* 

o 

MB 

lA 

9^ 

3 

< 

o 

A 

9 

3 

A 

A 

A 

9 

X 

MJ 

O 

MJ 

O 

t 

X 

• 

M» 

c 

*3 

.M 

A 

3! 

Ci 

<M 

o 

.f 

Mi 

-T 

9 

9 

o 

O 

M 

A 

O 

M 

O 

MB 

A 

r>> 

A 

MJ 

r*» 

A 

O 

O! 

X 

A 

• 

3 

• 

• 

o« 

>r> 

«V 

9 

Ai 

X 

A 

3 

MB 

» 

•r 

VM 

3 

3 

r>\ 

UB 

■T 

Vi 

j 

A 

Vi 

A 

Vi 

O 

A 

a 

a 

1 u 

Mi 

'i 

ar 

o 

(3 

tw 

u> 

MJ 

3 

MB 

*B 

o 

D 

■f 

9 

A 

O 

a 

3 

; ♦ 

M 

o 

Mi 

9 

A 

3 

o 

Ci 

9 

A 

3 

• 

• 

• 

MB 

MB 

CM 

3 

• 

MJ 

9 

MB 

•« 

V 

•4 

f 

9 

rB 

A 

o 

BJ 

— 

3 

3 

aC 

0 

A 

Ml 

O 

1 X 

1 

■r 

r- 

9 

■3 

3 

4 

A 

9 

■n 

:? 

JBl 

« 

/3 

A 

3 

9 

0 

Ml 

3 

4 

/ 

0 

3 

m 

M 

<B 

M 

3 

MB 

9 

a 

A 

A 

A 

4 

<M 

31 

MJ 

o 

rj» 

•B 

M 

0 

3 

Mi 

<3 

MB 

A 

A 

3 

A 

O 

9 

A 

3. 

! ^ 

.M 

Mi 

o 

• 

• 

■0 

O 

-T 

IB 

A 

Ml 

Mi 

< i 

X 

MB 

•i 

3 

A 

A 

C 

3 

O 

0 

O 

Ml 

N 

9 

3 

1 

o 

o 

BA 

MB 

A 

a 

o 

uB| 

Mi 

Vi 

MB 

r3 

3 

I 

j 

1 

MB 

MB 

V3 

1 

Oi 

<5 

•A 

O 

4- 

A 

A 

A 

3 

A 

9 

3 

A 

•B 

3l 

•r 

> 

MJ 

c* 

M) 

• M 

O 

3 

CO 

Mi 

3 

3 

sX 

ft» 

A 

A 

T 

A 

■-Bl 

1 

»/* 

>i* 

r» 

•a 

•7* 

3 

Mi 

■3 

3 

Mi 

3 

< 

MJ 

0 

A 

n 

*1 

• 

3 

A 

o 

3 

a 

C 

9 

CB 

Ml 

31 

! 

• 

UB 

WB 

O 

3 

A 

A 

9 

#B 

M 

! 

r\ 

o 

MB 

i 

© 

0 

O 

O 

3 

3 

o 

3 

o 

o 

M. 

3 

rj 

Ci 

Ml 

9 

M* 

C/ 

9 

rj 

o 

MJ 

9 

A 

• r 

Ci 

i 

MB 

.J 

.J 

rj 

3 

O 

9 

Ml 

oi 

i» 

3 

<M 

Mt 

s 

c? 

m 

Ml 

M 

a 

A 

•c 

it 

A 

9 

A 

3 

Vi 

Vi 

o 

•f 

0 

O 

J* 

E>» 

Mi 

(i 

A 

- 

o 

I ^ 

• 

• 

• 

X 

• 

OB 

>4 

Vi 

r» 

Ml 

<S 

r» 

3 

0 

M 

a 

•i 

'i 

9 

■V 

A 

M* 

CB 

•• 

f'" 

4< 

-M 

«3 

•B 

3 

2 

3 

9 

3' 

■ 

( 

lA 

< 

•J 

VM 

A 

MB 

A 

•* 

Ml 

A 

A 

31 

« 

<0 

tlB 

MB 

O 

o 

0 

A 

A 

O 

o 

9 

vM 

f 

3 

A 

A 

MJ 

A 

3 

.« 

... 

MB 

1 J 

0 

Mj 

•f 

rt 

A 

BJ 

M 

3 

»•, 

A 

MB 

3' 

O 

O 

« 

O 

9 

JB 

bA 

A 

JB 

3 

o 

Mi 

OI 

9 

0 

3 

r-. 

M 

A 

A 

ai 

o 

Mi 

o 

A 

'3 

3 

1 

MB 

» J 

JB 

Ml 

0 

.Mi 

3 

• 

Vi 

3 

<3 

O 

X 

•3 

n 

9 

• 

9 

TB 

MB 

X 

MB 

o< 

eo 

MB 

o 

u 

Vi 

A 

9 

— • 

Mt 

O 

*•> 

a 

BCN 

Mi 

r| 

0 

Bi 

3 

1 

•• 

A 

"• 

A 

oi 

W 

9 

ai 

M« 

A 

MJ 

A 

•J 

A 

a 

> 

«/) 

M 

X 

> 

A 

•m 

« 

> 

A 

• 

ae 

> 

A 

bm 

B/i 

A 

A 

U 

A 

A 

ft 

A 

Ai 

, « 

y 

«t 

2 

2 

X 

4 

3 

jf 

O 

c 

•t 

a 

s 

3 

X 

c 

»• 

X 

£ 

.C 

• 

3 

X 

9 

X 

X 

X 

a 

! ^ 

X 

X 

«t 

< 

0. 

J4 

x 

4 

a. 

X 

4 

4 

9 

r 

3 

X 

X 

X 

2 

X 

a. 

Oi 

4 

3 

a. 

Ji 

4 

3 

a.1 

X 

VI 

71 

Ml 

X 

X 

i/B 

A 

A 

& 

A 

A 

4 

A 

t. 

A 

A 

A 

4 

134 


figuie  i6d  - CUEUi  ulla  Kin  I AUTCBEG6ESSIVE  StaVICE  UtES  AND  ECISSCN  INIUI.  TAi.ULA'TlCN  0?  SAEflE  SIAIISTICS  EOB  IHE 
LISIAIAUIICNE  OE  CkCUlAIEO  INTEAABBIVAL  HUE:!  AT  M-eCOO  AND  N- 1UOOO  (J  2 ABE  X 4).  CUaULAIED  SEBVICE  TIBES  (S  2 AND  S 4), 
BAITIBG  IlafS,  Nllli  CS  BJIHOUI  ZEBCS  (UI  ♦ 0,  BI,  I EOE  CASE  B = 2500.  5000.  1500,  lOCOO)  , CUHUIATED  BAITING  IIHES  (B1  BAB  etc.) 
ABC  IBi  AVEBAuiD  OIEEEBEACES  BilUEEB  S2  AND  X2  |D2);  OBTAINED  EBCB  *‘500  BEIIICAIICBS  OE  TUE  BUN  SISSBSB;  BX=2.97:  BSO  . 


I 


I 

I 


2.  As  for  the  steady  state  of  the  value  cf  E[  W ] and  E[il] 
a large  N Has  required  for  high  t and/or  r,  also  for  the 

steady  state  of  the  distibution  cf  H and  H a large  value 

of  N Has  reguired.  Informally  testing  their  distributions, 
by  using  the  Plotting  Subroutines  and  the  values  of  their 
parameters  skenness  and  kurtosis,  we  may  conclude  that: 

(i) . The  distribution  of  H,  given  that  W>0,  has  an 
exponential  form,  if  the  plot  results  in  a straight  line 
under  EXPLT  and  the  '/^^and  ^^parameters  have  the  values  of 
about  2 and  6 respectively  (see  figures  23a-23d) . This  is 
true  only  if  the  steady  state  has  been  reached  while  it  is 
not  evident  what  their  distribution  is  if  the  steady  state 
has  not  been  reached.  Thus  observing  the  figures  24a-24d, 

where  the  W has  not  reached  the  steady  state  yet,  we 

1C000 

can  obtain  nothing  about  xhe  distribution  of  S. 


Note  that  for  the  M/M/1  queue  the  steady  state 
distribution  of  W,  given  that  W>0 , is  exponential.  The 
distribution  is  not  known  for  the  correlated  queue  and  this 
was  one  cf  the  objectives  of  this  study. 

(ii) . The  distribution  cf  U has  the  normal  form  if  we 

get  a straight  line  under  NORMPL  and  if  the  values  of 

skewness  and  kurtosis  are  around  'O'.  The  nornality 

assumption  holds  for  the  steady  state  ccly,  while  the 

non-steady  state  distribution  looks  like  an  exponent icnal. 

Figures  27a-27c  give  us  the  feeling  that  U tas  an 

10000 

exponential  form,  since  a straight  line  appears  under  EXPLT 

(see  figure  27b)  and  the  parameters  have  the  values 

2.3  and  7.3  respectively.  These  are  close  to  the  actual 
values  2 and  6 of  the  exponential  distribution.  We  can  see 


35 


from  the  flcts  that  as  the  U gets  close  to  the  steady 

state,  it  leaves  the  exponentional  form  and  comes  closer  and 

closer  to  the  normal  form.  Figures  28a-2dc,  23d  sbcN  us 

that  ¥ has  already  reached  the  steady  state  and 

10000 

furthermore  its  distribution  has  eventually  taken  the  ncrmal 

form,  since  a straight  line  appears  under  MOBMPL  (figure 

28c)  and  the  skewness  and  kurtosis  have  the  value  .7  and  .8 

respectively.  Compare  its  distribution  with  the 

distribution  of  ¥ of  figure  27b  where  ¥ is  not  in  the 

10000 

steady  state.  Again  we  may  state  here  that  convergence  to 

the  steady  state  of  the  distribution  of  «•  and  ¥ requires 

much  larger  N than  convergence  to  the  steady  state  of  their 
mean  values  for  the  same  t and  r.  And  also  the  higher  t 
and/or  r the  larger  N should  be.  As  an  example  we  can  see 

that  even  though  the  steady  state  (in  value)  of  ¥ for 

S#25#98  type  has  teen  reached  for  N=10000,  the  steady  state 
in  distribution  has  not  been  reached  yet.  See  figures 
24€-24g  where  we  can  observe  that  the  distribution  of  ¥ 

10000 

starts  to  leave  the  exponential  form  and  to  go  tc  the 
ncrmal  form. 
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NUlIbKil  Of  UIIUI.III!D  lAlRS  • SOO 


Figure  27t  - QUEUE 
AMD  EOISSOM  INPUT. 
CUMULATED  AMD 


liITH  EAR1  ADT08EGHESSIVE  SERVICE  TIMES 
EXPONENTIAL  PLCT  (EXPLI)  CF  THE 


AVERAGED  WAITING  TIMES 
S#5C#90;  ai=500  REPLICATIONS,  HX=.2;  8S  = .4 


i FROM  TEE  RUN 

2 0 00 
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e:>11NAIE0  t>AKAnt'tek;i  Ot  1>ATA  ; GAHHAI  • 00  GADOA^  « 7..)28S0Sl!  00 


NUNfiEB  Oi  OBUEBED  (BIBS  - SOO 


Figure  27c  - QUEUE  iillH  EAE1  ADTOHEGBESSI VE  SERVICE  TIMES 

AND  POISSCN  INPUT.  NORMAL  HOT  (NORMPL)  Cl  THE  CUMULATED 

AND  AVERAGED  WAITING  TIMES  T!  FROM  THE  BON  S#50#90;  ai  = 500 

2000 

REPLICATIONS,  RX=.2;  HS=. 4 . 
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-SCALE  : *•*  » 0.9o0t-ii1  UNITS 


rigure  Ida  - wUEUE  WITH  EAE1  AUTOREGRESSIVE  SERVICE  TIMES 


AND  POISSCN  IJiPUT. 
BAITING  IIHES  7 


HISTOGRAM  OF  THE  CUMULATED  AND  AVERAGED 


20  00 

REPLICATIONS,  RX=.2;  RS=.4 


FROM 


THE 


RUN 


S#50#25;  ni  = 500 


mo 


MUdBEB  or  OBDBBEO  UIBS  • $00 


Figure  28t  - QUEUE  BITH  EAE1  AUTOREGRESSIVE  SERVICE  TIMES 

AMD  POISSON  INPUT.  EXPONENTIAL  PLOT  (EXPLT)  OF  THE 

COMULATEC  AND  AVERAGED  WAITING  TIMES  m FROM  ThE  RUI' 

2000 

S#5C#25;  !B  = 5CC  REPLICATIONS,  RX=.2;  RS  = .4  . 
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ErillRlIED  rAHtntTEkS  Uf  OtT*  : tiASaAl  « 7.  C*B)<>2  • 7.92«eBU-01 


figure  28c  - QiJSOE  lilTH  EAR1  ADTOREGRESSIVE  SERVICE  TIMES 

AKD  RCISSCN  INPUT.  NCRHAL  ELOT  (NOHaPL)  CE  THE  CUHtLATED 

ANE  AVERAGEE  WAITING  TIMES  1 FROM  THE  RUN  S#50#25;  b=50G 

2000 

1 HEPLICATICNS,  ax=. 2;  RS=.4  . 


-sLALe  : *•*  * o.^bae*ot  units 


Continuing 

the  analysis  of  the  results  ue  present 

telou  the  re 

suits 

ue  got 

from  all  runs,  concerninc  the 

values  of  EC'S]  for 

various 

values  of  t. 

r,  N. 

N 

t 

r 

ECS](correlated)  E[S ] (M/M/l ) 

2000 

.25 

.25 

.3720 

.3332 

2000 

.25 

.50 

.4326 

.3332 

2000 

.25 

.90 

.7578 

.3332 

10000 

.25 

.95 

1.0576 

.3332 

10000 

.25 

.90 

1.7077 

.3332 

2000 

.50 

.25 

1.2106 

.9985 

2000 

.50 

.50 

1.5811 

.9985 

2000 

.50 

.90 

5.0785 

.9985 

10000 

..50 

.95 

9.3196 

.9985 

10000 

.50 

.98 

21.0542* 

.9985 

10000 

.95 

.25 

23.5651 

18.3528 

10000 

.95 

.50 

33.3367 

18.3528 

10000 

.95 

.90 

123.9350* 

18.3528 

10000 

.95 

.95 

196.2460* 

16.3528 

10000 

.95 

.98 

342.9240* 

18.3528 

10000 

.99 

.25 

62.3285* 

53.3371 

10000 

.99 

.50 

77.9370* 

53.3371 

10000 

.99 

.90 

193.9113* 

53.3371 

10000 

.99 

.95 

271.7713* 

53.3371 

10000 

.99 

.98 

426.6750* 

53.3371 

Sote  that 

these 

values 

have  been  scaled  from  the  values 

ettained  in 

the  siiulations 

by  multiplying  by  83.  Thus  for 

t^O.SO  and 

r»0.25  the  table  gives 

E[ 'H  ] <correlated)  as 

1.2106;  the 

value 

in  Figure 

23d  is  3.026x0.4=1.2106  . The 

values  aaiXed  * hav«c  net  converged. 
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A Study  on  these  results  gives  us  the  following 
chservaticns : 

1.  For  given  traffic  intensity  t the  value  of  E[¥]  is  not 
constant,  as  it  is  fcr  H/H/1  queues,  but  it  depends  cn  the 
value  of  the  correlation  paraaeter  r.  Thus  we  see  that  the 
value  of  ECiI  ] increase,  as  the  value  of  r becomes  larger 
and  larger.  Furthermore  from  figure  29a,  where  E[  T!  ] is 
plotted  versus  r we  can  see  that  the  rate  cf  increment  is 
not  constant  (therefore  is  net  a linear  function  of  r)  but 
it  increases  with  r and  eventually  gees  up  infinitely  as  r 
gees  to  1. 


2.  For  given  correlation  r,  we  see  that  the  value  cf  FCI] 
increases  with  t,  the  same  as  happens  in  the  a/M/1  queue 

(see  figure  29b).  Eut  regardless  of  the  value  of  r,  E[  ]!  ] for 
Il/a/1  queue  is  always  less  than  for  the  EARMA  model.  Thus 
we  may  state  that  the  weiiting  time  of  the  autocorr elated 
service  time  model  is  greater  than  the  corresponding  waiting 
time  of  the  H/M/1  model. 
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3.  The  values  of  £[17]  for  the  type  S#99#98  nodel 
increases  with  N following  the  curve  in  figure  30  a fact 
that  is  known  to  happen  in  the  tl/H/1  queue  also.  Me  can  see 
froB  that  figure  that  the  rate  of  increment  comes  dcwn  as  N 
increases  and  hopefully  should  come  to  zerc  as  ¥ reaches 

the  steady  state.  Ccmparing  then  the  M/H/l  queue  with  the 

correlated  queue  we  can  see  that  the  M/M/1  queue  converges 

faster  than  the  correlated  especially  when  r is  high.  Thus 

we  can  see  from  figures  32a-32c  that  both  the  value  and 

distribution  of  the  ¥ for  M/M/1  queue  have  reached  the 

10000 

steady  state,  but  the  ccrrespcnding  ¥ cf  the  correlated 

10000 

queue  (figures  33a-d)  has  not  (at  least  in  its 

distribution) . (See  the  skewness  and  kurtcsis  values  for  HI 

EAB,  H2  BAB,  H3  SAfi  and  M4  BAB  in  Figure  33d.  These  appear 

tc  be  decreasing  tc  zero,  but  have  certainly  not  reached 

zero  by  51=10,000.)  Looking  at  figure  31c,  we  can  cbtain 

the  result  that  neither  ¥ for  the  correlated  queue,  or 

10000 

¥ fcr  the  unccrrelated  queue  has  reached  equilibrium. 

10000 

This  is  because  the  skewness  and  kurtcsis  values  are 
definitely  net  close  tc  zero. 
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•4UMbLK  QF  UHiiLHEU 


OAMMAI  • I.UI^I^CE  00  GANHA2  « 00 


Figure  31c  - QUEUE  WITH  EAS1  AUTOREGRESSIVE  SERVICE  TIMES 

AND  POISSON  INPUT.  NORMAL  PLOT  (NORMPL)  OF  THE  CUMULATED 

AND  AVERAGED  WAITING  TIMES  U FROM  THE  RUN  S#S5#95; 

10  000 

01  = 500  REPLICATIONS,  RX=3.8;  RS=4  . 
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Figure  32a  - M/M/l  QUEOE.  HISTOGSAM  OF  IKE  CUMULATED  AND 
AVERAGED  WAITING  TIMES  ;»  FROM  THE  RUN  S#50#95;  m = 500 

I 000  0 

REPLICATIONS,  RX=4 ; RS=8  . 
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iiuaiit.il  01  oauKUKO  (tits  « SCO 


in 

■N 


Figure  32b  - M/a/l  QUEUE.  EXPONENTIAL  PLOT  (EXPLT)  0? 

THE  CUMULATED  AND  AVERAGED  WAITING  TIMES  ¥ FROM  THE  RUN 

10  000 

Si50#95;  in  = 500  REPLICATIONS,  RX  = 4;  RS=8  . 
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MUItUeU  Uf  OgUEBCD  (iIRS  > SOO 


Figure  32c  - H/M/1  QUEUE.  NORMAL  PLOT  (NCRKPL)  OF  THE 

CUMULATED  AND  AVERAGED  WAITING  TIMES  Ti  FROM  THE  RUN 

1 0000 


S#5C#95;  a=500  REPLICATIONS,  RX=4;  RS=3  . 
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WAITING 


3a  - QUEUE 
SON  INPUT. 
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HISTOGRAM 
FROM 

0 0 00 


AUTORE 
OF  THE 
THE 


GRESSIVE  SERVICE 
CUMULATED  AND  AV 
RUN  S#50#95; 


TIMES 
ERAGED 
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REPLICATIONS,  RX=4  ; 


RS=8  . 
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lUItGEII  Ot  ORbliSEU  lAIHS  • iOO 


Figure  33b  - QUEUE 
AND  POISSON  INPUT. 
CUHULATEE  AND 


WITH  EAR1  AUTOREGRESSIVE  SERVICE  TIMES 
EXPONENTIAL.  PLOT  (EXPLT)  OF  THE 


AVERAGED  WAITING 
S#50#95;  b=500  REPLICATIONS,  RX=U;  RS=8 


TIMES  i FROM  THE  RON 

10000 


157 


e:i1JBAIkO  up  data  : CAIIBAt  ■ 9,B7d092B-0l  GARflA2  ■ 00 


MUBBAi  cr  OBDEIED  (BIBS  • SOO 


Figure  33c  - QUEUE  KITH  EAH1  AUTOREGSESSIVS  SERVICE  TIMES 

AND  POISSOM  INPUT.  NORMAL  PLOT  (NOHMPL)  OF  THE  CUMULATED 

AND  AVERAGE!  WAITING  TIMES  H FROM  THE  RUN  S#50#95; 

10000 

1=500  REPLICATIONS,  RX=a;  RS=8  . 
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il9UCc  jjj  - CUEti  lllb  EtSI  AU'ICBtGSESSIVt  SEBIflCE  HIES  ANL  ECISSCK  INEUl.  TiEULATICN  OP  SAEFIE  STATISTICS  EOB  THE 
llSlilEUTICAS  CE  tUEUlATEU  I BTE  A ASA  I V At  TI.TtS  AT  ll="5000  AHb  N- 10000  (I  i ADD  X O),  COHULAIEO  SERVICE  TI8ES  (S  2 ARO  S 1)  , 

MAiiiiio  Tiais,  Kiia  oa  .ithooi  zesos  (ui  ♦ o,  21,  1 pos  case  ii^isot,  soco,  Tsoc,  locco),  cubocated  raitiac  tibes  (oi  bab  tic.j 

ABC  IBE  JVEBAGIC  CIEEEBEACES  bllBEEH  S2  ABO  X2  (02);  OEIAIBEO  FBCB  ■•SOO  SIEIICAIICNS  OF  THE  BUB  SASOfSS;  BX*4;  BS*e  . 


r- 


Another  result  obtained  iron  the  simulation 

concerns  the  cooparison  of  convergence  cf  tt  and  ¥.  From  the 

prograas  we  have  run  ue  observe  that  7 ccnverges  slightly 

faster  than  H (at  least  for  low  t and  r) . Also  the  standard 

deviation  of  D is  much  less  than  of  H.  In  figure  24d  ue  can 

see  that  1 has  already  reached  the  steady  state  but  not 

10C00 

tt  . Ccipare  also  the  values  of  their  standard 
10000 

deviations  (.005  versus  .2).  Figure  34  gives  the  results  of 
the  type  Sv50#98  aodel  for  N=10000.  We  can  see  that  while 

"S  is  already  close  tc  the  steady  state,  the  fact  that  K has 

not  reached  it  yet  forces  us  tc  continue  the  simulation  by 
increasing  the  value  of  N. 


i 


L 


2*  SiJSiiJ  with  Cross-Correlated  Service  lime 


Ihe  model  for  the  crcss-correlatad  .gueue  is; 

<X  .}  = g . 

i 1 

IS  } is  EAfiaA(1,1)  over  { E , fiX/aS,  HX/ES,  ...  } 

i i i i- 1 


lhat  is: 


S = 
i 


1 ^E.  + A. 


where 


f 


A = 
i 


i-1 


w.p.  ^ 

w.p  (1-^) 

w.p.  r 


rA  +6  SX/HS  w.p.  (1-r) 
i-1  i 


Again  we  chose  the  values  .25,  .50,  .95,  .99,  for  t 
and  .25,  .50,  .90,  .95,  .98  for  r and  we  run  a total  cf  20 

runs  to  ccver  all  ccahinaticns  of  t and  r,  naming  each  run 
as  SC#tt#rr. 

Analyzing  the  results  obtained  from  the  siiulaticn 
we  observe  the  fcllowirg: 

1.  The  value  of  the  traffic  intensity  t and  the 

correlation  r is  the  main  factor  for  convergence  of  R and  B. 

That  is,  high  traffic  intensity  and/or  high  correlation 

require  N to  be  large  in  order  for  W and  B to  reach  the 

steady  state.  Thus  simulating  the  models  with  t<.50  and 

rS.95  for  M=10000  it  was  possible  for  M and  i to  reach  the 

steady  state.  In  figures  35a-35d,  36a-36c  we  can  see  the 

justification  of  that  conclusion.  The  values  .1076  and 

.0957  of  R , and  B represent  the  steady  state  of 

10000  10000 
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the  SCf50#25  model.  On  the  other  hand,  regardless  of  the 
value  of  t the  models  with  i>.98  did  not  reach  the  steady 
state  for  N=10000  and  regardless  of  the  value  of  r the 
models  with  t>.95  did  net  reach  the  steady  state  for 
H=1C,000  . figure  37  gives  us  the  non-convergence  of  the 
model  SC#99#25.  Because  of  the  high  traffic  intensity  and 
in  order  to  get  the  steady  state  of  that  model  we  have  to 
continue  simulating  by  increasing  the  value  cf  S. 
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Figure  35a  - QUEUE  WITH  CEOSS-COHaELATIL  SERVICE  TIME 

SEQUENCE  AND  POISSON  INPUT.  HISTOGRAM  OF  THE  WAITING  TIMES 

W WITHOUT  ZEROS  FROM  THE  RUN  S0#50#25;  in=500 

1 0000 

REPLICATIONS,  RX=2;  R3=4  . No  OF  ZEEOS=279. 


NUHtlKIl  Cf  ORueSliC  ItlBS  • 22\ 


Figure  35fc  - QUEUE  HC  TH  CROSS-COH R ELAT E E SERVICE  TIME 

SEQUENCE  AND  POISSON  INPUT.  EXPONENTIAL  PLOT  (EXPLT)  CF  THE 

MAITING  TINES  W WITHOUT  ZEROS  FROM  THE  RUN  SC#50#25; 

1 00  00 

in  = 500  REPLICATIONS,  RX  = 2;  RS  = 4 . No  OF  ZERCS=279. 
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NUneiLil  0(  OUDERED  FAIRS  > 221 


Figure  35c  - QOEaE  WITH  CROSS-COSRELATZC  SERVICE  TIME 

SEQUENCE  AND  POISSON  INPUT.  NORMAL  PLOT  (NORMPL)  OF  THE 

WAITING  TIMES  H WITHOUT  ZEROS  FROM  THE  RUN  SC#50#25; 

1 00  00 

[D=500  REPLICATIONS,  RX  = 2;  RS=4  . No  OF  ZERCS  = 279. 
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tiquie  Jbd  - QUtUE  WITH  ChOSS-COBBELlTEO  SEfcVICE  IIHE  EiUUEECE  «kC  POISSCE  IMPUl.  IlfiUllllCN  Of  StEPlE  STillSTICS  FOH  THE 
EISlalEUllCkE  CF  CEtULlTEO  INTEFtBElVAL  Hcl£:>  il  N-ECQO  AMD  NclOOOU  (E  2 iNC  1 4),  CUMULATED  SECFICE  TINES  (S  2 AMO  S 4), 
kAITlMC  TINES.  KITH  CS  klTUOUr  EESOS  (MI  ♦ 0,  MI,  I FOB  CASE  N-2tOC,  5uOO,  <500,  10000),  CUMULATED  BAITIEG  TINES  (Ml  BAB  etc.) 
AMO  lUE  AMiiACEO  DIFIEBEACES  fllTktEM  S2  AMO  12  (D2) ; OBTAINED  FBCM  ■“SOO  BEEIICATIOMS  OF  THE  BUM  SCI50I2S;  il=2;  aS-4  . 


2.  High  trafic  intensity  and  correlation  affects  the 
convergence  of  the  distrituticn  of  M and  ¥ also. 
Furthermore  the  steady  state  of  their  distributions, 
requires  larger  N than  is  required  for  the  steady  state  of 
their  expected  values.  In  examining  that  conclusion,  look 
for  example  at  the  results  of  the  simulated  model  SC#50I98, 
where  although  H and  ¥ have  already  reached  the  steady 

state  for  N=  10000,  their  distribution  has  not  converged 

yet.  The  above  results  have  been  obtained  from  Figures 

38a-3ac,  3Sa-39d  where  W and  ¥ are  presented.  We 

10000  10000 

can  see  from  these  figures  (38a-38c)  that  although  the  value 

cf  W has  reached  the  steady  state  (see  Figure  39d) , 

10000 

its  distribution  cannot  be  characterized  since  the  plot 


under  EXfLT  (figure  38t)  is  not  a straight  line 
parameters  and  are  far  away 
desired  values  2 and  6 . The  same  results  can  be 
from  the  figures  39a-39c  which  deal  with  ¥ 

10C0O 


and  the 
frcm  the 
obtained 
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SEQUENCE  AND  POISSON  INPUT.  HISTOGHAM  OF  THE  CUMULATED  AND 


AVERAGED  -iAITING  TIMES  ¥ FROM  THE  RUN  SC#50#25;  in  = 500 

1 0 000 


NUBbEit  Cfr  OiDCREb  lAlRb  » 500 


Figure  36t  - QUEUE  WITH  CROSS-CORP.ELATEE  SERVICE  TIME 

SEQUENCE  AMD  POISSON  INPUT.  EXPONENTIAL  PICT  (EXPLT)  OF 

THE  CUMULATED  AND  AVERAGED  WAITING  TIMES  W FROM  THE  RUN 

10  000 

SC#50#25;  ai  = 500  REPLICATIONS,  RX=2;  RS=4  . 
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ESIlHATED  PABEBkTeBS  Ot  DATA  : CAHBAl  * 1.7bJ515E-01  GARBA2  « > 1 • 89S75«C-01 


Figure  36c  - QU20E  WITH  CHOSS-COfi HEL AT F I SERVICE  TIME 

SEQUENCE  AND  POISSON  INPUT.  NORMAL  PLOT  (NCRMPL)  OF  THE 

CUMULATED  AND  AVERAGED  WAITING  TIMES  "H  FROM  TEE  RUN 

1 0000 

sc#50#25;  ib  = 500  REPLICATIONS,  RX=2;  RS=4  . 


-scAie  : *•*  « o.96oe-oi  units 


CCHkPLATFt  CUtUt 
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J7  - QO£UE  KIIH  CEOSS-COaHELATED  SERVICE  TIME  SICOENCE  ini  ECISSOt  IMtUl.  TitOlAIICR  OF  SIEFIE  STATISTICS  FOB  TEE 
IISTilcUlICtS  CF  CUISIITEC  INIERARBIVAL  TINES  II  N-SCOU  AMD  10000  (I  2 AND  > V),  CUMULATED  SERVICE  TIMES  (S  2 AMO  S 0) , 
MAniAG  TIMES,  MI7H  CB  VITUOUT  EEROS  (MI  » 0,  HI.  I FOR  CASE  M^2S0C,  SOOO,  7S00,  10000),  CUMULATED  HAITIEG  TIRES  (HI  DAB  etc.) 
AMD  THE  AVEBA5ED  blEEEVEACES  NETHcEN  Si  AND  12  (02);  OfiTAIMED  FBCM  1>S00  BEEIICITIOMS  OF  THE  RUN  SCf99l2S:  91-2.97;  BS-j  . 
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Figure  38a  - QUEUE  WITH  CROSS-CORRELATED  SERVICE  TIME 
SEQUENCE  AND  POISSON  INPUT.  HISTOGRAM  OF  THE  WAITING  TIMES 
H ailHCOT  ZEROS  FROM  THE  RUN  SC450»98;  !n=500 


l 0000 

HEPLICATICNS , RX=2;  RS=4  . 


No  OF  ZER0S=245. 
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xuaaEH  cf  OBoaMEC  Kina  > 255 


Figure  38b  - QUEUE  Wir H CEOSS-CORB ELATED  SERVICE  TIME 

SEQUENCE  AMD  POISSON  INPUT.  EXPONENTIAL  PLOT  (EXPLT)  Cf  THE 

HAITIMG  TIMES  H WITHOUT  ZEROS  FROM  THE  HUN  SC#50#98; 

100  00 

18=500  REPLICATIONS,  HX  = 2;  RS=4  . No  OF  ZEHCS=245. 
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ESIIHATED  PAHABETEBS  OE  DATA  : GAHHA1  » 1.34115652  00  GAHHAl  « I.099609E  01 


Nunutn  Clr  OBUKfteO  IAIRS  > 2SS 


Figure  38c  - QUEUE  (illH  CROSS-CORBELATEE  SERVICE  TIME 

SEQOENCE  ANE  POISSON  INPOT.  NORMAL  PLOT  (NORMPL)  CF  THE 

HAIIING  TIMES  W WITHOOT  ZEROS  FROM  THE  RON  SC#50#98; 

1 00  00 

ib  = 500  REPLICATIONS,  RX=2;  RS=4  . No  OF  ZERCS=245  . 
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Figure  39a  - QD2UE  HIIH  CROSS-CORaELATEE  SERVICE  TIME 

SEQOENCE  AMD  POISSON  INPUT.  HISTOGRAM  OF  THE  CUMULATED  AtND 

AVERAGED  BAITING  TIMES  if  FROM  THE  SC#50#98;  ai  = 500 

10  000 

REPLICATIONS,  RX=2 ; RS=4  . 
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Figure  39t  - QUEUE  WIIH  CBOSS-COHRELATEC  SERVICE  TIME 

SEQUENCE  AHC  POISSON  INPUT.  EXPONENTIAL  PLOT  (EXPIT)  OF 

THE  CUMULATEE  AND  AVERAGED  WAITING  TIMES  TI  FROM  THE 

10000 

SC#50#9a;  a=500  REPLICATIONS,  HX=2;  RS=4  . 
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esiIRATED  PARAOeteitS  op  DAT*  : GAfltlAI  > ^.bQ220f^t  00  GAflB*2  « 4,912<l09E  00 
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figure  39c  - QUEUE  WITH  CROSS-CORRELATED  SERVICE  TIME 

JfwUEMCE  AND  POISSON  INPUT.  NORMAL  PLOT  (NORMPL)  CE  THE 

J.eJLATED  AND  AVERAGED  WAITING  TIMES  T!  FROM  THE 

10000 

• •9(1;  >«500  REPLICATIONS,  RX=2;  RS=4  . 
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figure  39il  - QUEUE  MITE  CEOSS-COEBELIT  ED  SERVICE  TIHE  SEQUENCE  A ED  EOISSCN  INPUT.  lAfiUIAIICN  OF  SIEEIE  STATISTICS  FOB  Iht 
EISTSIEUTIONS  Of  COBUIATED  INIERABRIVAl  HUES  AT  N^'SCOO  AND  11=10000  (I  2 ANI  X N),  CONOLATED  SERVICE  TINES  (S  2 AND  S U)  , 
RAITING  TISES.  NITH  CS  MJIHOOI  ZEROS  (VI  ♦ 0.  HI.  I FOB  CASE  N=250C,  5000,  7500.  10C00),  CUMULATED  BAITING  TIMES  (*1  DAB  etc.) 
AND  THE  AVERAGED  DIFEESEACES  DETNEIN  S2  AND  X2  (D2(j  OBTAINED  PRCM  »=500  EEEIICITICNS  CF  THE  RUN  SCI50A5e;  fcX=1.5;  BS=fc  . 


3.  lo  test  the  distribution  of  M and  ¥ we  also  usee  the 


r 


plotting  subroutines  and  the  values  of  their  skewness  and 
kurtosis  parameters  as  we  did  with  the  gueue  with  dependent 
service. 

Ihe  results  obtained  from  the  analysis  are  as 

fellows: 

(i)  . Given  x.hat  H>0,  we  may  say  that  W has  an 
exponential  form,  since  the  plots  under  EXPLT  Subroutine 
result  in  a straight  line  and  the  values  of  (^^i.and  ^j_are 
around  their  actual  values  2 and  6 respectively.  For  example 
the  SC#50#25  model  (which  is  represented  by  the  figures 
35a-35d)  , where  the  steady  state  has  been  reached,  gave  us  a 
straight  line  and  the  values  1.8  and  4.7 
for  and  2.  respectively.  Ke  have  to  state  here  that  the 
exponential  distribution  comes  only  if  W has  reached  the 
steady  state,  and  we  cannot  have  the  same  result  if  the  H is 
still  in  a transient  state. 

(ii)  For  the  distribution  of  "ff  we  can  say  that  the 
normality  assumption  is  an  appropriate  one,  since  the 
graphs  cbxained  under  NOBMEL  result  in  a straight  lire  and 
the  value  of  skewness  and  kurtosis  is  around  'O'.  The 
results,  for  example,  obtained  from  the  SC#25#25  model  gave 
us  a straight  line  under  NOBHEL  and  the  values  .17  and  -.05 
for  the  skewness  and  kurtosis  respectively  (see  figures 
40a-40c)  . The  normality  assumption  holds  if  the  1<  has 
reached  the  steady  state.  At  the  begining  of  the  transient 
state  it  is  not  possible  to  characterize  its  distribution, 

while  when  U is  in  a transient  state  we  can  see  that  its 

2 

distribution  takes  an  exponential  or  a X form.  This  is 
left  little  by  little  as  U comes  close  to  the  steady  state 


\ 
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and  then  eventually  it  takes  the  nocnal  icrm.  The  icdel 
SC#25#95  is  a representative  example  cf  that  result.  In 

that  model,  although  ¥ has  net  reached  the  steady  state  yet 

(£1=10000)  we  can  observe  the  values  of  skewness  having  as 
3.5,  1.8,  1.6  1.5  and  the  values  of  kurtosis  as  24.7,  7.2, 

5.7,  4.3  for  N=2500,  5000,  750C,  10000  r espectivelly . That 

is  the  values  start  frem  a high  level,  then  ceme  down,  they 
pass  from  the  actual  values  of  for  an  exponential 

model  and  then  continue  decreasing  and  hopefully  in  the 

steady  state  they  should  reach  a value  around  '0*  where  the 
value  cf  skewness  and  kurtosis  for  the  normal  distribution. 
Ihe  same  assumptions  are  obtained  analyzing  the  plots  of 
that  model  also.  That  is  for  IJ=2500  the  plots  are  not 
linear  under  both  EXPLT  and  NOEMPL,  for  N=5000  a straight 
line  appears  under  EXPIT,  for  N=7500  and  10000  the  plot 
leaves  gradually  the  straight  line  under  EXPII  and  begins  to 
give  us  a straight  line  under  NOfiMPL. 
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CROSS-CORRELATED  SERVICE  TIME 
HISTOGRAM  OF  IHE  CUMULATED  AND 
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1 0 000 
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- f IMAlCC-^IIRft«tttR5'Or  DATA  • I * * CARnAt  ■ — Ufe89te6E-01 CAMMAe  A *9664  1 8E-02 


Figure  40c  - QOEOE  WITH  CROSS-COSRELATEC  SERVICE  TIME 

SEQUENCE  ANE  POISSON  INPUT.  NORMAL  PLOT  (NORMPL)  CF  THE 

CUMULATED  AND  AVERAGED  WAITING  TIMES  tJ  FROM  THE 

1 0000 

SC#25#25;  id  = 500  REPLICATIONS,  RX  = 1.5;  RS=6  . 


184 


i 


I 

\ 


1 


o o -•  o o 
o o o o 


o ^ 


-i  3 ^3  «0 


Q O I O p 


^ ! > 

Jt  ' «« 
< C 


-J  a T b 


3 


I I 

i I 
I I 

■ i 

! 

I ' 


! = 


I 'i 


I 


I 


I 


I 

I 


l»  i 
5 i 

•T 

r>j  ri 

-5 


^ I — 


Q U 

I 1 

<JJ 


07  ci 


tA  t <*\  — 


O ; 

i 


fta  3}  3)  (U 


A 


•>*  3 

oi.; 


h.  VI  — wf*  J 

eo  ^ o — — p 

•*  a — A*  .f  p 


•«  M a 


r 


c 


« ' -r 

' A 

- |«v 
O I ■'*• 

V-  • J 

V ( I 

I 


■ o 


- ^ 


< a ^ 

I -• 

W V)  ' wl  ic 


I I 
i I 


a ,ij  o o‘  o 


I ;5 


|U  -o  eoi 


i -o  i'^ 
: o *«» 


O !P  P — 


I" 


1 P jO  O f<^^^  3 

p Ip  o oi  p 


A Ip  P M 


I**'  ' Q 'O  P 


’I  ? 


6 - 


i 


P P P O' 


5 2 ::l? 


185 


Figure  «Cd  - QU£0£  Klin  CBOSS-COmtliTli:  S£I)yiCE  llflf  IIOUFHCS  A»C  fOISSCN  IMEUI.  IIBUIXIICH  OF  SIFFIE  STATISTICS  FOB  TUB 
IISIFIEUIICIIS  Of  CUBULIIEO  INTEBABBIVAL  TIKES  AT  B^SOOO  ABE  U- 10000  (I  2 ABC  X 4),  CUKUEATED  SEBVICE  TIKES  (S  2 ABO  S 4). 
BAITIBC  TIKES,  BITH  CE  NITHOUT  ZEBCS  (BI  * 0,  BI,  1 FOB  CASE  N^2S0C,  SOOO,  ISOC,  10CC0),  CUHULATEO  lAITlNS  TIKES  (Bl  BAB  etc.) 
AND  THE  AVEBAGEO  DllEEBEICES  BETUEEU  S2  AND  X2  (E2)  ; OETAINED  FOLK  »‘500  IE f I IC AI IC AS  CF  THE  BOB  SCI2SI2S;  BS=6  . 


4.  The  followiag  values  for  17  have  been  obtained  frca  the 


runs  with 

various 

values  of  t 

and  r and 

for  1J  = 10000 

N 

t 

r 

Etl] 

E[-K](fl/M/1) 

1000C 

.25 

.25 

.1169 

.3332 

10000 

.25 

.50 

. 1724 

.3332 

10000 

.25 

.90 

.3538 

.3332 

10000 

.25 

.95 

.4075 

.3332 

10000 

.25 

.98 

.4679 

.3332 

10000 

.50 

.25 

.3831 

.9985 

10000 

.50 

.50 

.5648 

.9985 

10000 

.50 

.90 

1.4589 

.9985 

10000 

.50 

.95 

2.  1109 

.9985 

1000C 

.50 

.98 

3.7196* 

.9985 

10000 

.95 

.25 

10.4672 

18.3528 

10000 

.95 

.50 

13.4784* 

18.3528 

10000 

.95 

.90 

29.  1312* 

18.3528 

10000 

.95 

.95 

75.1180* 

18.3528 

10000 

.95 

.98 

145.34  84* 

18.3528 

10000 

.99 

.25 

35.5986* 

53.3371 

10000 

.99 

.50 

43.  1616* 

53.3371 

10000 

.99 

.90 

98.9556* 

53.3371 

10000 

.99 

.95 

139.0167* 

53.3371 

10000 

.99 

.98 

219.7920* 

53.3371 

An 

uljzing 

the  above  results  we  may  state  conclusions 

as  follCHS 

(i)  . 

Fcr  a given  traffic 

; intensity 

t the  value  of  U is 

net  constant,  but 

depends  on 

the  value 

of  ccrrelaticn  r,  as 

the  gueue  with  dependent  service.  He  observe 

also  that 

13  increases  as 

r increases,  but  the  rate  of 

increaent 

is  larger  as  r goes 

to  1 , and 

very  small  for  low 
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values  cf  r;  therefore  TI  is  not  a linear  function  of  r. 
Figure  41  gives  the  flats  of  ¥ versus  r for  t = .25,  .50, 
.95,  .98  and  also  for  the  M/M/l  queue  (r=0)  . Furthermcre  we 
can  observe  that  the  value  of  1 for  the  correlated  queue, 
ccmparing  it  with  the  result  for  the  a/U/l  queue,  is  not 
always  greater  as  in  the  case  cf  dependent  service,  but 
depends  cn  the  value  of  r.  It  is  less  for  r<.50  and  greater 
for  r>.90 

(ii) . Fcr  a given  correlation  the  value  of  V increases 
as  t increases.  Furthermore  high  values  of  t cause  high  rate 
of  increment,  which  becomes  larger  as  the  correlation 
increases.  Figure  42  gives  the  plots  cf  1?  versus  t for 
r = 0,  .25,  .50,  .90,  .95,  and  .98  . Observing  that 
figure  we  can  see  also  that  fcr  a given  t the  order  of 
increment  of  U is  for  r = .25,  .50,  0,  .90,  .95,  .98 
respectively.  Thus  we  can  conclude  that  the  waiting  time 
fcr  the  cross-correlated  queue  is  generally  greater  than  the 
waiting  time  for  M/M/l  queue  when  the  correlation  is  high. 

5.  The  last  result  obtained  from  the  simulation  of  that 
model  was  the  convergence  of  W and  TX . (These  statistics 
alsc  converge  in  value  and  in  distribution  slower  then  the 
same  statistics  in  a/li/1  queue)  . The  mcdel  SC#25#98  reveals 
that  i converges  faster  than  W (see  figure  43) . Eut  the 
model  SC«S9#50  gave  us  the  conclusion  that  W ccnverges 
faster  than  U (see  figure  44) . Thus  we  can  state  here  that 
fcr  low  values  of  t U converges  faster  than  B and  fcr  high 
values  of  t the  opposite  is  true. 
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figure  UJ  - QUEUE  ma  CB0SS-C0BEEL»TED  SEBVICE  IIBE  SIQUENCE  tNC  tCISSCS  INEUt.  TiEUliTICB  OF  SiEFlE  ST4TISTICS  f08  THE 
EiSlElEUlICIIS  Of  CUEULiTED  INXEBiaBIVil  TIMES  »T  H«S000  AND  N=  10000  (X  2 ANE  I 4),  CUMULATED  SEEVICI  TIMES  (S  2 AND  S 4)  , 
NAITING  IIBIS,  NITH  CB  BIIHOUT  ZEBCS  (HI  ♦ 0.  HI,  I FOB  CASE  N='250I,  5000,  "iSOO,  lOCOO)  , CUMULATED  BAIIIBG  TIMES  <B1  BAB  etc.) 
AND  1M£  iVEBAGLO  OlffEBftCES  BETBEEN  S2  AND  12  (02);  OETAINED  FBCM  i^SOO  BEEIICATICAS  Cf  THE  EON  SCI25IS8;  BJI'l.S;  BS=6  . 
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Elguce  44  - glEUE  KITH  CBOSS-COHEEL  AT  E 0 SERVICE  TINE  SEgUERCE  EMC  ECISSCR  IHEUT.  TEBUIATICII  OF  SEEEIE  STATISTICS  FOB  TEE 
CISISIEOTIONS  Or  CVEULAIED  INTEAANhl  VIL  TIMES  AT  N-SOOO  AND  N==10000  (X  2 ANC  X 4),  CUMULATED  SERVICE  TINES  (S  2 AND  S lil, 
RAITIHG  TIMES,  AITH  C8  NITHOUT  ZERCS  (N1  * 0,  NI,  I FOB  CASE  N=2S0C,  SOOO,  iSOC,  10C0C)  , CUMULATED  MAITIKG  TIMES  (lil  BAR  etc.) 
ANC  THE  AVEBAGEC  DIFEEBEACES  BEINEEN  S2  AND  12  (D2)  ; OBTAINED  FRCM  ■•SOO  EEElICAlIOliS  OF  THE  EON  SCISSISO;  11*2.97;  ES*j  . 
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PROGRAMMER  : GFORIOS  J.  DANIKAS 

UNDER  DIRECTIONS  CF  PROF.  P.  A.  W.  LEWIS 


SLBROUTINE  JACK 


• » •>  » •»* 


••'OUJ  ></)«/) 


I 

(M 


- 2Q 

ai(/7 

^4 

» • • 

• • • 

♦ 

• 

IM 

» 

» 

a 

+ 

• • • • 

• • • 

A ^ 

c\i 

2< 

•zo 

o 

• • « • 

» w • 

fN4  to 

z 

•» 

O -1  X 

w 

^ » 

-H 

z 

to 

m m m 

4*  ^4* 

X 

r-  z»  i-» 

< 

Z 

-5  O 

► • «»•  • 

m m m 

1- 

♦ 

a 

«» 

II  fVi 

m m m 

o 

u 

• a 

• 

m m m 

a 

X -J 

(SJ 

K < •‘tfOUJUJClC 

Q 

UJ 

o 

< • • • • 

m m m 

LU 

a 

Q mC^KTi 

z 

K Q 

\£\ 

a. 

u 

oc 

in  tf> 

Z K-IX 

(/)  » ,1 

m m 

3 

o 

CM  (Ni 

►I  < * 

—lo 

• f»«  m 

• • • 

O 

X 

* LU-  <K0S«:=50 

IX 

m 

I 

o o 

(O^CsJ 

z* 

^ • aim  m 

• • • H- 

II 

o 

h- 

u 

1-  H- 

•J 

,>11* 

a m 

z 

< 

< 

— • • *0 

xvi 

z—  ••••••. 

m m 

O' 

3 

K 

LU 

o o 

— t/lioo 

^z 

• • • 

• • • o 

WO 

u o 

•-4XXXZ 

o 

LL 

U. 

_|UJ  1 

►c 

m m m 

* • • 

o 

•» 

o 

C/5Z 

cc 

X 

<\i  tz  — 

X 1 

•+- 

UO 

to 

A M 

*-•  -J 

UJ 

LU 

<M  o • O— 'VIh-'— 

Oa: 

IT* 

1- 

1- 

<\i<M  rn 

•O-Z^M”?  •XLUOC 

A 

Z'« 

<r 

< 

•>  o 

<Mi-i*  5^0  II  Z 

M 

— Ci-Ql 

• 

• 

z 

to 

z 

Z 'Z  • 

-J  •—It  II  “1  (I  U</1  •• 

X »- 

Ui 

h* 

Ui 

M 

LU 

<X-^i-<Xa.  n 

z 

•J 

1- 

X 

h» 

J 

o 1 ”?x».— -jacM 

c 

»-*0^ 

• 

*CL 

(/) 

u 

to 

O •—*  •<MO  II  >-•  ai  o •II 

II  LO— 

z*  < 

1 '-•OC 

LU 

< 

Ui 

CVIZ  1 

N.Z  ••  c -»  >-• 

X(M 

o a. 

a 1 

QCeciNiCj  •-  (MZU.</> 

•-^-J 

20«->0* 

LU 

LU 

z II  >-«)M_l_|  II  II  (\J^>-«X_J_JI- 

— o 

00< 

M 

I zu 

►- 

1- 

tn 

^ |M4 

zu 

— z 

Z3 

II  f\J— . w 

II  (M<-i^O  <<Krot/i 

UJ,-< 

II 

II 

CL 

•J 

O. 

II  Z II 

-J  — _IQ  UOt/^ff^X-lh- 

xot- 

II 

z 

z 

(MO  « 

MO 

'-a< 

U.OIU.C’ 

o 

u 

l''OCJU.<'iU.OCJO  O 

<h- 

Q-IO 

'►-zz 

u 

u 

u 

_lZZ'-'U>->^ZO  o 

Uto 

****** 


CVJ  m 

(\»  -H  m 


UUO  UUU 


O 


WRITE  (6,200 


X 

r~ 


It  r- 

•• 


••II 


— -7  (/I 


< 

£ 


<o 

i/)-^ 

'-K 
< I- 
to 

NH(/0  UJ 


z 

— • X 

O-H  5C 

oo  o 

••  ••  "> 

•O'O 

UJ 

H 

UJUJ  Z3 
a. 
£ 

accc  u 

Z3  o 


m 

«b 

o 

•b 

•» 

• 

UJ 

X 

LU 

K 

II 

-o 

O 

< 

(M 

2 

2 

z 

•b 

< 

9* 

c 

• 

h- 

A 

o 

QC 

a 

L/) 

-5 

z: 

< 

< 

LU 

•» 

1 

> 

N-4 

CM 

«b 

• • 

2 

• • 

t/) 

LU 

H- 

zu. 

QC  •> 

LU 

< 

<UJ 

UJX* 

— CD 

CM 

N-O 

H-f^a 

•b 

«• 

t/) 

OO 

UJ  u. 

* VilU 

* 

LLi* 

2 

— 

2 

<1*  LU 

* <*. 

— 

«b 

q:-iu2 

* • X 

“J 

•• 

- X 

a^QC 

•Jl^  <. 

«-n 

f\l 

ax 

*r>  2< 

— *• 

O 

o 

UJO  •‘V. 

• UJXO. 

< < • 

K—  — 

ir>u.x 

SUJOir^ 

<— oc— 

-5 

LL  • 

L/^^4  »-4 

I-I-—  1 

»• 

on— 

X*  >• 

ai2<x 

l-H  UJ 

l/><>.OC 

Q. 

X 

<ZUitO 

Ouoor^ 

■H-h-*-*— 

a<OM 

Q.X*  < 

a ojo. 

— c/)l-— • 

M 

♦ ^ 

^LU  (/) 

»-<o  'n 

-J*  * QC 

CL2  ^ J-*  3 j 

1 

M 

• OU'V  -TX 

LUOC 

QC—  "X 

•> 

OO 

ron 

Q •+-!- 

X P-  -UMVIOCOCX 

— 1 1 — 

good 

a 

UJX«/)QClf'-<C3  X 

— OCOCQC-J 

2 

H-O-  3 

• LU  •••O'Oc/)  Z 

1 — — — •• 

1/1 


II 


— — — "VCSi 
-J I I PJ— 
•— — »!- 
tSTK  < 


£ 

QC. 

< 

a. 


LUUJ 


— </i-> 

X • 

z^-o 


»<Mw 

•HZ»- 

Q«l 

kOH* 

^/1£^/) 

X I 


p+-«*-a.—  oca: 

•-<w  I — jES 

a:h-»/x.«»/:oc  ■-* 

»-C— >-»-Q.C7  + 

+ I/ll/)— uo  O — 

f^l/)(\J  + ••JO  fOIV  — 

+l/)^II  II  IIP--II.  O Z—  — 

• II  1-10.1/)  "-“OiO  • —I—  — 

-j~J  h-  Q.— ,-i-H  • o OUJH-^XOt  II  fM 

- • •—  It  II  I •M-'O  O 0-J_l  O O 

II  O *00  II  • ••  • II  OO  CM  »•  II  II  KO-vOO 

—O  (M»-1*-tC>JCO  — •• 

II  II  •0-Jl/)i/>— ~)  II  UJ  -o 
II  <0)-a.Q-KI-l-  II 


z •• 

>-i»  X UJQ 
••p-uj 
oOO  -XO-O. 
Luo«  r— JD 
. ocuj  a 
ijij  - -oc 

z«  xo 
X <i/)inz 
<0  M— ■</)  o 
•4-XQCUJ  •• 
c/i<zm 
X >3:—  < 
X • UJ  X 
» X "i-i 

•J/)X 
V.X-  fo 
voo 


■etac  UJUJ—  - 
£UJ>*rri_i_j  Q 
a.o.<£UJ 
HUJ>»  —£>OQ. 
l/)£  <Z0CZ> 

LJ<U.  1/1  — U.C 

- o:0  OL 

< -JJHH/)0 

-a  •<•!»••«  ocz 

X-  >0*  * LUO 

UJ  — • — • 

-H  -C  - LU 

X "X  ••£  - 

• x.i-x<x 

O — rJ-lPIQC  — 

'^•T< 


< — IM 

•tH-— »/)</)0 

4’«/)»-Q.Q.Q 

u 

a 


<« 


o»-iz  z 


i/>  I/) «/)  r- m 3 o I- 3 ►- 


u 

c 


•o 

>o 


ujacuJa:u.Li.u 
ac3ac3— o 

- ui  in  c 

- — (M  c 


XCM  — 

— — >»v. 

CX.  ••'O 

X 

o— 

fM 

1 — 2t- 

l-l- 

I wUJQC< 

<< 

<<< 

< 

_ll-l-3£ 

££ 

£££ 

£ 

_j«i-i»-(ac 

a;o: 

OCQCQC 

QC 

<ii“a:ujo 

oo 

OOC3 

O 

LJI/)3QCLL 

LULL 

LLLLU. 

LL 

•tt 

* * «•  * * 

* 

i o 

o— 

OO  — 

CM 

o 

oo 

OOO 

O 

CM 

<0  0000 

00 

oou 


216 


u 

— • 

• r\u 

QXK 

uj  ra< 

1 

< O" 

z 

X'vt— 

X 

•-•>..1/) 

w 

h-  lu 

—* 

</)• 

«— 

Ol*  Z 

— »M 

* lU 

ro  1 

LU«  1X1 

1 z 

CO  CO 

zx 

—1 

— x— 

yt 

o<i> 

a. 

— * z 

zx««. 

z? 

-*  —X 

Uj 

► (M 

< 

to  1 — 

uc<x 

z 

p» 

•—  z 

ot- 

«*  x>«. 

l/)l—  < 

•> 

s 

1-  — 

a c 

o 

3 

< * — 

UJ(/) 

t/) 

1-  —CO 

H-Q. 

««</1  <-t  1 

UJ30 

•» 

— 1 z 

XCUJ 

. i/y 

zx  • 

4(X.Q. 

PM 

a. 

OW>  01 

QC.03 

maj«— (M  1 

< a 

— 1 < ■^rncM 

a.u.ac 

K 

UJ  —■ »•  * * 

OO 

<. 

•» 

— •'V(M  -*•!«■* 

uu  Z 

1- 

m ^ 1 

ti.aC3 

00 

£ 

X — Z 1 ^Ol 

mUJ 

3 

oO  •'tx— Z X X 

Zco 

•» 

—•4/3 

0^ 

■*XXH  > Xl/H/) 

xxz 

X 

h- 

xH—  0</1#— •—x-i 

ii:oa 

^ «> 

1 1-  »x— <M-» 

oZo: 

0-4- 

UOCdfO^ 

Z<  — x<  1 IMI-I- 

< u. 

—a: 

XI-  -*l-  1 Z*<«I 

"5* 

'-'X 

w*  * 

w</^  —<;ZX*I-|- 

• » (/) 

</y 

K 

z 

KU-U-U. 

—<  I- X «-—.(/)(/) 

* OC 

UJ 

< •> 

•Vw  xc/1— .*  I'O 

M ILi 

l-t-oo 

00  z x's.'». 

X t- 

u 

002: 

t/7 

OJI—  X-  llwXOO 

^ »UJ 

< 

X 

XX 

+ + + 

XQC  W +x,rr,^ 

XZ 

n 

+ 3 

1 

30  •I-— '*.00— XX 

Xl-I< 

•» 

Z«-4U0 

ooooO<»rs  l-XX 

ZLL 

-X 

—XXX 

h-  x<r)_< 

< 

— t » 

II 

-<•-•333 

II  II  II  O0V1X  1—  II  II 

XQ 

O >/1  OOO 

—»/)(/)  </1 

— 3*1/1 

Uj 

Oil  —COO 

II  X 

X— —1/1^*  — — * 

z 

Z 

• II  t-l  • • t 

II  II  II 

oisj-inooi-*  xcovor- 

•M 

o 

o*-*  xOOO 

*-*  II 

•»  •»  •‘tO^Z3  •*  •“ 

K 

MOO 

^1/) 

ivirn^t' 

O04/3O0<HX</)  1 004/) 

h-  Z 

3 

OO^f 

Z II  X"  II  II  II 

U.XXX 

•-M*— ♦-I  U3 

< a£ 

O 

z 

^3— 

<M*-333 

MPwwMP  n II  -a: 

X 3 

0^ 

UJU 

II 

<MQV0(/H/) 

h-H*-  I-K3 

QC  t-O 

03 

x< 

X <xxx 

<<^  <<H- 

a ixiz 

PMUJ 

Z3(_i  1-333 

O 

ZZ  K^-liJ 

U.  OCUJ 

00 

a<x 

X(/)U  i/)«/li/)</) 

o 

4/)4/)00«—  XX  </34/>q^i 

« « 

•M-  * 

ro 

C«J 

o 

0^ 

fM 

00 


217 


URN 


•k 

o 

•k 

D. 

UO 

UO 

l/J 

U 

u. 

a. 

Q. 

M 

•» 

u 

X 

9k 

c_; 

p‘ 

UO 

o o 

U7 

.-<o 

u. 

oo*^  ^ o 

tn 

u. 

X 

X 

PM 

00<n  ^ in 

OJ 

o 

3 * 

OO'-'  'T  - 

X 

UO 

00>-«  ^ 'T 

3 

•» 

••Q. 

P"t 

a.  ^ 

1/7 

«kH 

X 

XU 

X 

ww(QO  O CL 

X 

3 

3 

PM 

X^Xi-O  X 

•» 

•k 

3X30  • X 

•k 

tp4 

U-V)0 

f>i4 

X 

•-♦X 

m4 

► ^ o •* 

l/J 

LJ 

U70 

UO 

CO 

CO 

3 

1. 

3X 

N-* 

UJ 

N-* 

LX. 

Qd 

ooo«-<xc>  ► o 

Z 

•k 

•k 

•k 

o.ooroo.iO'T  o 

3's. 

•» 

UO 

•k 

UO 

ooo'— x^—  in 

1/7 

o 

X 

3X 

UJ 

-tOXOi-«C  *•  •“ 

liJ 

a 

O 

PM 

ill 

•— X —• 4- 

►O 

•-H 

ffk 

I 

*Z 

X#..*— 'X  ►CD  ••  >-» 

^ • 

X 

•-<cd 

m 

< 

Q.XX3*'  —OX 

X 

X 

X 

u 

O3u.l/)0  *003 

X 

X 

•k 

X *> 

3 

3 

f^ao 

3 

a. 

» fi  ^ tiOOUN  ^ ^ 

X 

•k 

UO 

PM  ^ 

</) 

Q- O O*-* 

•» 

X 

I 

idX  ^>4/7^  XXf^ 

M 

• 

OOOm.—  X'-xO 

Qd 

Qt  com—  XXOO  XX—  XXX 

X 

a:o 

00~— ItDXl/OO 

X 

•k 

» • *<n  i£^  «.  *en 

3 

►oo-4<o  xxm 

3 

«k 

^ » ^uom ' ^ •►xuo^^ 

M 

GC. 

r-ti-«i-)</)D.  wXO  ► 

PM 

UO 

uo 

mm  — ao(M  II  II  *in— r.kixi.4*  ii  ii  * 

o 

X— —CDX— 3</)'J’ 

Xw 

X 

o 

XOin  It  11  «i_^ro  II  II  — — ^ 

CQ 

3i/>*H2:ujoc/> 

CDf-i 

•> 

z 

z *0  m — ►->  *o  X 

uo: 

OCUJ  • 

3</>3  ro  fO 

X 

lOX  ► 

CC 

*q£C  </'I—  * ►XO  Xh-  * ►X 

31 

0.1/)  10 

► ItlO  vw  *— 3 

•“CO 

30</)0 

mi—*  *x  — 1—  * *x  — 

Ml 

301 

X •• 

—X 

X «k  MM  fk  MM  O^*** 

— </> 

O—  ► *Ol/>UO  ► 

CJV'> 

cOO^^ 

CJUJ«» 

O— fOcDOO  — 

^CO 

**  •• 

Cd 

»a:o*'**i/)  — imo— —x  — i-ix 

Ml 

•<Oo 

*00t-40  inir 

*^*-X 

s mmcD  mmffi  xxco  xxcd 

i/>: 

a 

^OOcf^O  ^ 

CO  f 

uox 

X ~ 

* asffiz  CDOZ  CDCQZ  corns 

X! 

wZ«. 

1/1 0 O— — '<•  .t  o 

CD^ 

3U, 

*X  3—  3 3—  3 

Ml 

<<Ot 

o«-*Oi/)»*o— '>^in 

fH 

-ZoQt 

ZQC— "m  .«.wm— w— >m  —— m 

xo  * 

— •f-*0D»HCOX(/)  ► 

(X  a. 

►—•'-XQ-— XX  4- 

— X 

M w>f— ^ 

-> 

< — 

lOlOl/7ZXlOOO  — 

^co 

• • 

UJX< 

3 3 U.  </)  UJ  CD  1/1  lO  3 

'-X 

PM 

<13L 

too 

O 

CJ 

O ►-»-4f\JUO,J#-4CMUO^MCMUO,J*-«<M</7 

M 

-JO^Ct 

•k 

Zaato 

o 

o 

o 

O 300<  •00<300«t  •00< 

o 

<o< 

o 

^4 

m 

i\j  •fon.ii-  coniK  .fDfMh-  incMh- 

rsj 

a:ooi 

z 

•> 

•k 

•k  » HJJ«-*4  » ftUJ  P>  ^UJm-4  » MJJ 

LUcC 

o 

m 

m 

m 

^ com>ocDu7U0kOC0X4nkO(SXin'O^ 

'O 

zaoi 

OC^ 

w 

.-  Xw«.  s^ww 

W 

lo  z 

UJi/> 

z 

o *- 

oz  < 

o 

CUk^UJ'^UJwUJwUJ 

OliJ 

OiUJ 

x< 

Q 

a 

Q 

V-  01-3  Ql-3  QI-3  QI-3 

K 

)-X 

x»- 

< 

< 

< 

MM  <C^J 

II  M 

xz^ 

Z"- 

a< 

lU 

UJ 

UJ 

Od  cLUjQC<.u.uJod^cLUJQd<u.ujod< 

X 

►-»-< 

— Q 

uo 

Od 

oc 

(X 

3 — X3U*-'0C3O  — 0C3U*-0£3U 

~»3 

******* 

* 

* 

m ii>  uo 

in 

<M  fO 

uou 

\ 

u 

< 

218 


KS,  KSlt  KXt  KXlt  QS,  CS1«  QXt  QXl 


— .00 


«/)x 

m.  *U.U. 


adZZ— Z-»Z 

z 

» .-z  ••z  •• 

••-JJ  ^ ^ 

a.a.a.x3v)Z> 
XXX3  =) 
lUUJUJ  «■  > 

»VH  •«*H 

•XXV/IV) 

-<^^3333 


•» 

fk 

U.U.  ^ » 

to 

u> 

X 

X 

vox 

u» 

IL 

u. 

u. 

••  XDO 
VOXII 

ouaac 

zx 

aoc.  » » 

cc 

cr 

QC 

tx 

• • 

0 

0 

0 

a 

f 

i 

1 

h- 

1— 

K 

=5 

3 

3 

(/>X  » ^ 

< 

< 

< 

«»  «<^>x 

^f-4Q.a 

a.u.xx 

1.1 

XXlULLi 

•J 

UJLLi 

< 

< 

< 

< 

*•  •• 

u 

0 

OJ 

0 

OOXVOX 

ILU.0.0 

— V-I-W 

0 

0 

ocaaoc 

• 

• 

• 

• 

oaao 

0 

0 

a 

0 

l-KI-l- 

lJj«i»UJ^UJ<««UJ< 

3333 

• 

• 

• 

• 

in<<<< 

• 

^ • 

M • 

M 

<0 

00 

oo 

X 

X 

Q. 

«>o. 

»o. 

o<<<< 

oouuu 


O rl>-« 


t/1  10  X X 

UL  X u.  U. 

Ci  o o o 

z ••z  -z  -z  - 

•«/)  •(/)  (X  tX 

o o u u 

Z X I I 

Mota-rio.  f-io.  f^a: 


-«•  ><■ 

O0Xi/)O0 

>* 

0^  0:0:0 

o.<-njzx 

X'v'v*  ♦ 

• •x  • 

<ZZfriZ 

Z 

II  II  II  II 

II 

-HOsIfTlX 

It  II  II  II 

^ot/)x 

xzzzz 

I/^XOO 

zzzz 

zzzoouo 

COOQQQQ 

o.a.a.zzzz 

xxx<<<< 

ujanucLa^ac.oL 


“>«-><JOOOOU 


(/IX 

0.0. 

-(XX 

1-40101 


O • 400 
0iC7OU-0j 
•ujoj  • • 

• • 

(.0  ill) 

a.(/)X(/)xo 

uoooorin 

— wwO 


U.U.U.U.U.O 

— •t-if.d  o 


or  »<3<ytO  •ox 

01 1-HOI  OOi<-4  oj  o 

• • • < 

^ ^ • • f 

• *00  • »0  • 

(ooovixxxxi/)  00  00  O 

00.0a  00.  OQl  II  II  II 

X X X X Mil  II  II  II 

a^OJ— 'Ol—OI-aOl  ZZ  Z 


•-•Ota-aQCi-iOl*- 

u.  u.  o u.  a£oia'oi<z<i 

a— I <(/Xl/)0(/)0 

■»•■»•»■* 

in  o 

aO  in 


weAR 


Mf'BAR 


p» 

p 

• 4/) 

► X 

pH 

Q. 

•» 

a 

4/) 

X 

X 

pH 

X 

o 

<3 

40 

UJ 

X 

lii 

A ^ 

03 

03 

W 

p 

«>i» 

p»  M 

40 

p»  ^ 

X 

<•12:4/) 

N>« 

<f-42:x 

• ^ 

(fi 

• 

03 

2:XCC3 

mm 

<x:3< 

rH 

(x:3< 

pH 

P pp 

< UJ 

40 

< UJ  •> 

X 

► 40 

-J  « 

•*  X 

-J  *— 

UJ  » f-l 

U. 

UJ  •»  pH 

Li. 

PH  U, 

40  «pp 

Li. 

X pp 

0 

0 

40 

X 

It 

H 

•*  It 

^ Li. 

► II 

^ u. 

-J 

*H 

pH 

pp 

:s 

40  3 

II 

X 3 

II 

<^<> 

uu 

Uj 

X UJ 

•» 

X LJJ 

p 

z 

4.J>0 

z 

z 

40  X 

z 

^ 3 

4/1 

«> 

► 4/) 

X LL; 

► X 

X UJ 

» fM 

Z 

> nmm 

Z 

^ 40 

f-#  X 

40 

^ U) 

X 

•*  X 

40 

X (S 

3--» 

40 

0-^ 

X 

— (D—r 

A 

H«  mm 

3 

M .>» 

3 

r 3 

X r> 

mmmm 

ODOcO 

0 

0D040 

0 

^ 0 

•H  0 

^ pp 

40 

40 

40 

X 

• • •» 

• 

• t •» 

• 

0 • 

0 

0 • 

0 

h- 

3 

}-  •‘>-•-4 

3 

X -J3 

< 

X -13 

U)^<:;4o 

UJ 

(J^OX 

UJ 

OC  P.H&UJ 

X • 

QC  »<UJ 

X • 

•40  •a. 

f 

• X •4X. 

f 

<o>  • 

QL^  3 

<o>  • 

Xp»  3 

40 

u. 

U. 

UJ 

< UJ 

UJ 

^ UJ 

f^ 

^ ^ • 

pH 

UJp^  • 

LiJpp  • 

lO  -^O-J 

40 

X oX-l 

X 

.J  40 

X 

pp 

o 

— 

X— 

3^— JXOO 

X 

«pp  40 

-i<-^X 

^ X 

h- 

40  40  CL 

X 

X Q. 

-J>40 

40  .JX 

-l>x 

X^X 

UJ  a.  X 

UJ 

Cl  X 

< UJ 

< UJ 

O 

t-i  »X  (UJ 

f 

tX  tUJ 

• 

4J  • 

-1  > 

u • 

-J  > 

o 

«>CUJO 

Ow 

OUJO 

0 

pH^O 

< » • 

pHwO 

^ » • 

-IZ  z •> 

Z 2 * 

z 

40  .JZ 

0^4»Q 

X JZ 

4jPH#C 

gO< 

<x 

< ►<•-« 

< 

CD40< 

Z 

a3x< 

z 

Q.  •«/) 

•Q. 

• »-<  •X 

• 

XO.  • 

xa  • 

ppp< 

X t<0  03 

X 

X 03 

— DX 

-JP-*  • 

— IDX 

^pH  • 

LU  CQ  -J 

1 UJ 

03  -J 

1 

40  UJ 

10  Uj 

^<’-p 

0 •-< 

0 ^ 

0 ^ 

>40 

0 PH 

>x 

•»  ••  <0>  * 

» x> 

0 CDfH 

P# 

0 03ph 

• 

• • *0 

• 

• • •u 

• 

•pH  • 

•pH  • 

p 

a 

— O-iO  +0— 

a^to  +0 

►-40pp*ar 

•40#  • 

►-X— 0 

• X#  • 

1- 

UJ 

•^Uj  UJ  » 

UJfM 

UJ  UJ 

Uj 

OLi.^UJ 

H-CD  3 

OU.PHUJ 

3 

Z2l-<-*+  + 

• 

w • ••  •1—2  •w 

; • 

• V • 

OXSLU 

• w • 

4r>X3UJ 

^ ^ a. 

^4 

»-«  UJ 

»pai« 

• 3UJ  • 

• ^LiJ  • 

4040  40 

X X -z-* 

fH^40pH 

40Z 

^.^XpH 

402 

t02  »/)X 

•a 

rOO.  •>Q.^4O4OG0 

•a.  »a.i-«xxo 

4040  0340 

•40  4040 

•XXCDX 

• X XX 

II  It  XXOU 

0«-'0«/3 

<3  1 

ou-<ux 

at- 

x<3.  OOX  > 00X0.  oox  * 0 

II  II 

II  (/I  V >1 

II 

X ^11 

X M 

40  It 

X It 

X It 

HM  O II  It 

II  w 

•> 

II  -X'- 

— 0 

--pUJ  ^ 

M «- 

It  hpUJ  w 

M ^i'U. 

II  00_l  2 22  2K2  2 2 2 

lU  lU  UJ_J  LU  UJ-iliJ<^UJ^UJ 

rJOOZUJUJ<U.  zu-  LL  ZU.Z<U.  U.  ZU.O  O.  ZU.<U.  ZU.<Uu  Zu.<u.  ZU. 

«MV)22i^*^>>-<  V11-*  >-•  <->  Xi—O  •-*  Xm>m  Xn 

u 

O 

**  #*  W * « * 


u 


u 


u 


u 


220 


221 


» » » • i>  11  «<  M M II  i(\i^rg^<\js^i  II  II  II  i>^»^i 

rv4f»>.i-rNj>»"— «»— *■»••»•>  r » 
-•w»,w— ms— — w^www i-KsjfT.^— sxs2:<->c>ifr>*»-^c\jrisr 
OQXX«>0«/)  eD01<n(OZ5EXX«/5t/l———OOXX«/’KO— 

COCDCB  SZXZZXXS£SZZOOZSZ:S  SZSZ  CQaSfiCdDCO 

333  333033333333  3333  CQCDcQcr)  3333ZZZSSZZZ 

</1V)V7V0  l/></)«<0t/733333  333  «/)V7tOV/133333333 


csim'«r<ocDr\t'>r(\j'>r^rom>r— rvjm>r<ooooc3uooooououuo«jooc)0oooooooouo 
m m — fH  <M  cvifo  it>  m fo 't  •4' in  m m uMf>  m -o  >o  >0  >0  o >0  «o  <0  r- r~  in  uMn  tn  >o  <0  <0  o <o  vO  ^ o r>- 


MUiiiiiiijliiH  mill  ji|iii  HI  HI  Mill!  ii|  uinnii  fii  i^ljii  HI  III  ijiii^  til  III  m ulMiiii  fii  111  ii|  m iinimlJli  III  1^1  m 111  1^1  Ilf  1 


oc  oc  oc  a QC IX  (X  oc  oc  (X  QC  a a a ix  oc  oc  Qc  a:  a oc  Qc  o:  c£  a£  a a:  oc  a oc  oc  oc  oc  a o:  o:  a oc  QC  oc  ac  a:  oc  a£  a:  IX 

33333333333333333333333333333333333333333333333 


222 


1 


<NiCNJ<MCg<MCM<NjrM<NJ<MCNJC\JOJOa«NJ<MCsiCMrNJrsJ<NJ«\J<VJCM<NiCSJCMfSjrs| 

(iL  a.  u- u*  u.  u»  u.  u.  u.  a.  a.  u.  u.  u.  u.  u.  bu  u.  Oi  u.  u.  a.  u.  u.  u. 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

■ <■  «||  li  1 1 li  1 1 II  1 1 IP  Ili  II  Ip  Ii  <■  ||  <1  <1  *■■  *1  li**!  <»  l»  1 1 )i  I 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
m lo  ir\  irt  uMn  trt  irttn  irMcin  UMO  tn  uMr\  kP  in  IT.  u>m  tr  tfMr\  iDirt  IT  ID 


i-HrsirO'r 

cco^iXaC  ocacac.(x. 

<«<  <!<<< 

—tf\Jir>>fiDcoeoto  fM>»’f\i>^i-»njfosfcDcQcDaarM^ 


(\i>4’ZZS£ZZXX 

XXXl/)l/1  XXV1«/1  XZXXXX 

33333<<<<<<<<  3333<<<<<<<< 


i<0«/)(/1(/)k/)3333  32  33  33333  32 


XX  •■  ••  •XXXXXXXXXX  ••  ••  ► fXXXXXXXXXX 


^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ >1*  >0  ^4^  ^ ^ ^ *0 

• • • 

cococo 

U.U.lLf 


•-•II 

XM 

•» 

(Mi-n 

» fM 

• 

X 

in  » 

II  :7 

u.* 

m 

^ • 

x^  • 

• II 

•M 

inxx 

in 

xin 

II  X 

>.  MNiN.XO.  X 
NifclO  » X» 


•*ff^  fKim  » 

II  • 

h-«  f»  mxm  m 

• Xi^U)  v-4 

•»  or  * ^ 

• II  » KNJCJ  «C 


II  XX  in  ^ 
^tn  u\x  ^ • 

• MNi</)*  •• 

Xi-iOC  II 
»(n  »<•  II  1/7 

X •«■  «/)  a 

om  uj 

II  XkOQ.  LU 

• inoo 

*••  1^  !•  • 

>o  l/>  OOC  <»  — • 

M II  QlfMs  •>  * 

mxtn  in 


• i/)  »ti.m  * •» 

XX  • *0 

II  - in  • 

►•  • • ►— o 

2 *<Nj  mm  * »•-* 

• Xt—  II  U.LL»  - U- 

m •■  *•  II  (M 

•■  ••  •■X.  . II  ,H  ^ 


xm»  ot.  II  3 

- II  o«  •— m 
•4‘  II  t-  _J  ••x  • 

«o  “ ► i-i/»xaxmmm 
.-i»  ■ mxooom  '-•u. 
•—i-i  <3mxx«  ^ *>m>o 
m.  II  • ocac 


<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 

Z£XZXSXZZZZZ2£X2r2ZZZX£2S;ZZ£ZZ£2 

oc  oc  oc  oc  a QC  a (X  ac  a:  o:  ae  o:  oc  Of  o:  IT  ac  <£  O' o:  a:  ix  oc  oc  cx;  o:  <x  o;  o:  a£  (X 
oouoaooooooaouatjouaaoaoooooooooa 

U.  U.  UL  Ik  U.  U.  li.  U.  U.  U.  U.  11.  U.  kJL  U.  u- U.  U.  U.  U.  U.  U.  U.  U.  U.  U.  U.  U.  li.  U.  U.  Li. 


oo  C3  ^ ni  .r  PK  ^ «M  fTK  »i^<M  m >ii  00  <M  .f  m <M  in  (M  m *T  ^ <*> 

m .o  r~  cj  fM  IT)  m m r<i 'T -r  < m m i-<  .H  r\j  <M  ro  fo  in  fn  m m 

«M  »M  cij  nj  fM  iM  nj  <\J  fM  cgrunj  <M<M  nj  (vi  fsiiM 


'</)  •>.  ► •x— » 
xx>.xxm 
m m •mm 


Hi— H'H’ 

<< 

<<<< 

XX 

XXXXQ. 

ococ 

acoeacacoo 

oo 

OUOO»-2 

11.  li. 

iLikikikmui 

****** 

oo 

l-«M(-ll-4 

oo 

oooo 

•HtM 

riKMinm 

o 


223 


< 

> 


£ 

3 

CO 


fk 

CM 

a 

w 

m 

CL 

a. 

•• 

•m. 

X 

X 

a 

X 

o 

at 

aj 

X 

m 

o 

# 

# 

UJ 

VO 

o 

a 

o 

m 

•• 

o 

00 

+ 

GO 

— 

o 

3 

o 

o 

«> 

M 

o 

H- 

oc. 

u 

— 

a 

a 

t» 

z 

U 

to 

CL 

o 

OD 

•^00 

-J 

o 

M 

It 

HN 

a 

3 

1 

ro— . 

H- 

(h  » 

0-* 

A 

a 

*-< 

a 

o— 

A 

a: 

A 

m 

O 

X 

•— o. 

o 

ir\ 

•» 

— o— 

ai 

— XZU. 

a 

•» 

fO 

A 

X 

-*Q£» 

o.aj 

a 

o 

0^ 

D 

o 

O 

•-10.0 

•k 

X 

«> 

M 

o 

o 

K 

1 o 

a. 

UJ*  3 

1- 

w 

H^CC 

X 

CO 

O 

— o. 

A 

0. 

o 

o 

o 

o— 

— 40ac 

0^ 

X 

o 

VO 

•» 

1 -»♦ 

<•—«/) 

< 

a 

X 

<✓) 

• 1 

cr 

O 1 

^ • 

M • 

X 

• 

< 

1 11  1 ^ 

o 

0-4 

>- 

a 

» 

to 

#i4 

a 

*- 

K 

O'*.' 

-»o 

1 o 

< 

a 

•j 

a 

*S  — 

oo: 

• 

• 

UJ 

• 

X 

• 

CO 

r 

^4 

< 

-4« 

fO 

a 

» 

-OcOO 

r<j 

</)^ 

• (MOC  (/} 

UL^i 

0,.hcn»-#  *- 

CO 

• 

• 

^ O Q.  H 


—O'—  ££'— 


a 

lO  UJ  11  II 

UJ 

Z3  II 

3 

a 

a 

o 

z 

2 • II  — 

z 

ZZ»4#  II 

z 

z 

^ t 

o^- 

« 

O II  — 4 —4  > 1^ 

— 

C JOC^— 4 

— UJ 

-JO  *3 

K 

—4  ^ O — <-4 

1- 

_J  1_  — — 

3 

K 

•-O 

XX 

II 

3 

1/1— II  ox  + z 

3 

Zco  II  *-• 

zz 

3 

(/) 

O 

Zi-i  w4-taC3Xa£ 

O 

oz  — ^— 

U.>-4(£ 

O 

Z II  — 

M 

fHQca' 

GG 

UJ— O 1-10.1/)— 3 

Q£ 

-JXUI— _J  '-U.-<H'3 

o: 

-JUJ 

S 

xro  xi-o 

(O 

<XX-4-J 

r»-o 

<x-l 

•—OOCULU  3UJ^ 

UiO—— <U.O 

-3 

ai-«<aaa 

U' 

VO 

ocoa.— o cooeui 

(/) 

a£oau.o<-<o 

uoGa 

to 

0GO>^X^ 

o 

224 


11  CONTINUE 


1C  ISl  = IS+Il+KPl 


225 


o 

CJ 

in 


00 


•z 

9* 

— s 

X 

X 

o 

> 

;► 

o 

QC 

QC 

••in-* 

UJ 

UJ 

—— o 

VO 

-J 

-JZO 

< 

omm 

*» 

ac 

QC 

OZ-' 

X 

X 

UJ 

< 

03X 

> 

> 

ccz 

l-l 

»«* 

<UJ 

z 

*— 

q: 

QC 

-iU 

3 

1/)  » - 

< 

< 

3 

UJ»-* 

l_JUJ 

UCJ 

•» 

Z 

-oo 

a 

CC 

t-t- 

> 

—mm 

<.  «<* 

a 

cc 

o— — 

cc 

CD 

<11. 

UJ 

ozco 

X 

X X 

0.0 

I/O 

0X3! 

Ui 

or 

•» 

z 

Uiz 

m 

—3  ► 

H* 

K VO 

xo 

z 

—1/5  — 

>-►- 

> 

X O 

< 

< ^ 

t- 

UJ  »o 

socsxx 

a:< 

oc 

— m 

« < 

oo 

< 

-o— 

^CD  ^CO  ^ 

u.»-< 

—os 

3XSXLL 

u. 

m 

om 

oo 

lO 

UJSUJ3— 1/5 

><- 

cc 

o—  •• 

XX 

UJ 

-* 

z 

Z UJ 

«— 

->o 

< 

om— 

aa: 

-* 

V/) 

—1/5  — — 

— 

zo 

CO 

0X0 

</) 

X 

e xm 

OZ 

—30 

» » UJ 

» 

zz 

-I-  -t-o  - 

X 

— i/>m 

(/)(/)  M 

X 

>— >—ZX 

Ul 

o< 

u.  — 

Q^cx: 

UJ 

••  •• 

x<o£<a:uj 

•*  — 

UJ 

> 

-X 

•» 

xm 

UiSUJS  — 

— — - 1 

mco 

a: 

— *o 

• • X 

UJUJ 

to 

vo  •* 

3^- 

UJ 

— o 

XX  UJ 

-f-J 

— xm  a 

«/) 

CO  - 

C£QC  — 

u. 

— — 

c/^X 

-3  -SUJU. 

— U'  UJ  -J 

UJf-  • 

om— 

X 

xm 

>UJ>UJ  — 

X—  o 

(Ot-iZ 

•* 

0—0 

— — — 

UJUJ 

•»-  ♦ 

—z— z ► 

uiu.  ^ ^ m 

< 

> 

0X0 

zz  -u. 

9» 

-III 

axxxu.  - 

z 

— xm 

o— 

-> 

M A 

<i/)<m 

XS-f-fXX—  + 

<ujO 

a 

-3— 

— — 

*— 

f-Ul  >>  — 

OinO 

<O3V0QZ 

Z3f-— Z»— ot— O 

<a: 

M/)X 

— — 

-UJ— < X— XX  lum— -j 

XUQ> 

ouoo 

ZZ 

P-tOlUi 

•f  X — 

— z <3mxuj<  uoui— u 

< 

oooo 

o 

oo 

c 

o oo 

m3  <m  U.3 

CM  vox 

ao. 

W H M 

fO 

m h-<M 

II  II  II  II  11  II 

ooci 

z 

—XX 

•» 

- 3 - 

II  II  II  II  II  II 

o « - 

o 

>0«fNfV 

UJUJ 

ITS 

sO  <<o 

— QQ  XX3 

000£*-« 

•mm  • • 

•f 

— w 

QQ  J-J>>>>UJ3»- 

«ujm  z 

Z^«iN 

-J-JOQ—X— XZUJ— • 

« 

oz  o 

UJ  OUi 

^COX 

UJ  UJ 

(mccz3:xujxujxz< 

«nzz 

UUJUJ  xokoo 

K M II  M 

_J-J 

fMUJUJ 

o 

1-  -Jl- 

(Mm  3 3m<(/5<mm</5  3 

Moeoc 

II  — 

II  -J-J 

< 

— «J— 

i« 

UJ^I-I  LJUJQ^UJUJ 

acx-i-i 

<< 

u 

UJ 

ac.  <x 

u 

1— aiuj 

oc>-Q  uccxoconazinx 

ouo 

o 

ac 

3 U3 

O 

* * If  «■ 

o 


uuu 


u u u 


u 


226 


•» 

• 

u.  «• 

•» 

••  "X  • 

X 

II 

X*  in* 

> 

in 

a 

V) 

M •»!! 

X 

UJ 

o 

> 

V) 

X 

mo:  90: 

•<to<o 

UJ 

•k 

• 

OcQ'HX 

</) 

X 

> 

•» 

U.*  • 

X 

•k 

X 

at 

m 

•k  •>»  fk 

> 

< 

o 

- X X 

•» 

(M 

mum 

QC 

«> 

V. 

II 

< 

at 

m 

•*a.  ^ 

< 

«» 

UJ 

QL  •CC.  • 

a 

3 

•• 

• 

UJCT'VtQ' 

<i 

ro 

V)U.»  u. 

o 

«k 

« 

•• 

• 

3 

h- 

►o 

•k  «»  Ck 

•J 

XU. 

V) 

X. 

X » 

< 

tn 

a 

X m 

• 

H 

3 

•» 

UJ 

m II  II 

o 

_!  + •-< 

•k» 

OJ 

•k 

• 

X < 

•k 

V) 

•o-i^x 

O II 

>k  it 

• 

n)M 

^(/)  <UJ< 

X. 

• <o«i 

mm 

II 

»COZCQ 

V) 

•» 

03^3 

<-4^  II  11  11  H II  11  |l  II 

“t- 

-J  1 3Z</)Z 

►oc 

X 

U.U. 

(/>  IJJ3 

: 3 

v0« 

in 

u. 

«— 

OZ 

•»  m 

M 

^ fk 

«IUJ3 

1 •vl/)  •>  •—» 

fk 

►X-  X 

•• 

s 

Z XoiX 

1/1  »'X 

• in  m 

XX 

1 

O •+-UUI-0 

UJ  >0 

im* 

II 

M f*4 

II  ► f 

o z 

z a:< 

< 

II 

rtixm 

H- 

.013 

; -“30 

«k  ^ 

• 

> •>  • 

k'NXm  cc  XX 

o • ^ 

vx</) 

> tr 

X zro 

MOOttO) 

^J■  XX330Xt/1X 

• ZOH-< 

> • 

>M  » 

Ui  • • 

11 

atu.<u. 

m 

• "OO  XX3 

Oi/)  m2 

o— *-•  «.3a:3 

«-  o 

< - 

m>j-(/i</»  30 

• < 

V>C£>UJ<tiJ 

^U. 

-> 

• tk  m> 

o 

i-u  • mm 

• +UJi  + 

a zt-<»-z 

: z 

^ X M 

>• 

• ’ 

>k^— kU.4n 

UJ  .J 

V OC  OCX 

•■X-* 

u,  •* 

» X 'su. 

m ww 

-JO  •♦oc 

ZXX  'V<  1 <l/T^</)0 

«^X  II  O II 

rom  "O. 

••J  < 

>> 

o 

X 

kO  nJ 

•k 

•xm 

OXo^cO 

>>  >s»-a  at  1 X 

O •ll  • 

•mop'ix  3f<i  •oi'J 

.Q£  QCtijZ<  >. 

,n 

' (M 

1 OJ  -z 

aiu  Q^^t/ii-ajjxxm 

rn  • 

(T|*-»-IX  • 

> Ck^ 

' •^ZV.SI^MU.rOXlNJCNJfVJCMkUnKVJfM 

< <3 

.XcCUjO 

o o 

Oro»  i/i<-i»Hr<iMinf-«<Mn-4’U0a-4<\iro 

3 II  3 II 

t-icQ  II  II  <3uj(/)rsi 

cn  UJ 

<M  • U. 

• ^ 1 

LL 

(M  rj  CM  <M  «M 

II 

II  II  <3  3 1-0—  f 

•»  • 

»^«k^«k»«k«k 

WZ<.>  OCUJ 

3 — • II  «*  -Jj 

P-  o 

r- 

UJ 

1-  q:o 

II  II  II  — 

1-^  <CQZ 

-J-)  It  — M 

h-H- 

t-i-p- 

l-l 

-HI- 

U.<lLCOX>-< 

— —XX— -J  -'UJ 

UJ 

UK.< 

<<< 

<4<<UJOIUJUJUJUJUJUJ 

i-*3*-*33l- 

X I/O  . -J  X </n  w— *7 1- 

h-  h- 

HXX 

XXX 

x: 

Z 

xx-j'^xx—oa-j— <- 

— oeot 

ofotot 

ot  at  at  a£  »-< k- 1-1 . . 1-* « 

a 

o oj -1  o X a. ->  x QC 

OC  U.U 

OtUQ 

ouo 

u ui-i  a at  at  at  a at  at  at  at 

u 

</)  VI 3 2 (/)  (/)  3 2 o o 3 

3 •-'O 

30.0. 

U.O.U. 

U.U.U.U.333333  33 

* 

* 

* * 

(M 

moo 

O.^fPI'T 

<M 

^oo 

ooo 

OOOO 

«-«CSl 

<Mrsj<M 

u u 


227  j 

Ji 


^ ^ -4- 'O  ^ 'T  ^ '4- >0 
• ••••••••• 

(Ni(NJ(MCMr>JC\irsJOJCsl<Ni 

l<p  l>  ll  I 

U.  U.  IL  U.  U.  kU  11.  IL.  U.  Lk 
OOOOOOOOOO 


i|  II  II  II  II  II  II  II  II  II  II  II 


OO  X(/)  CO  zx  s 
rxxx3:*cix«/irxx 
3030  XX330 

»/>v1  30 

»/)«/> 


xxxxxxxxxx 

in  m ID  u>  (n  (D  lo  irt  ID  IT) 


(Ni  CM 

(X  a. 

< < 

DICO  CMCMIMCDIM 
(Ni(M  (MXXX 

X</>  X</)  XZ 

t-KO  (->-0 

XXh.!—  XX^i*-* 
33«  33<< 

<D</533  CDID33 


» ^ ^ ^ 

X »xxx  ••  »xxx 
>i-x>4-'t>i-xx>j'><-.a- 
^ .O  ■4" -O  >0 -O  ^ .O 
<0  'O^ 


X 


♦ 

~ X 

X 

o « 

UJ 

i/) 

O X 

O 111 

•k 

ox 

-<  + 

«» 

3 — 

w 

LJ 

1/5  ^ 

kl 

M 

—■»  u. 

M 

•• 

cc 

O—  II 

ox 

O 1 

X 

1^1  1 

U. 

HI 

X zu. 

» 

Ill 

u. 

+ •• 

••  3 

OX  • 

O Q— 

fM 

a 

o*.-. 

1 • 

c 

O X 

K 

3 

• 

u.  • 

< 

ZU.—  1- 

z 

xxo 

lua  • 

CM  »-< 

o 

UJ 

*Z-. 

tl  3 

Z 

z <-< 

11 

OXlX— 

^ wl 

1- 

— 3 

A 

3 

ZiD  M 

O 

oz  — 

a: 

XUJ'~>_J 

03 

XX^-J 

3 

a— «iu. 

O ' 

1/) 

OOlLO— 

o • 

u 


t^fOID 


— U.  U.  ' 

(Vjv 

D ID  ID  >0  >0  ID  ID  >0  m3  ^ ^ ^ • 

4"  <M<M<M<M<M  DJCMDIDJCU 

•o  >«  r>- r»  r- r^r- r- h- r- 

► 

UJ  lU  lU  111  111  ILI UJ  UJ  UJ  HI  uu  U,' <<  <<<<<<<<<<  < 

MH  ►-•oc  oc  oc  oc  oc  oe  a cCQC  a oc  ac  QC  o a 

X oc  oc  a:  fX  oc  QC  on  oc  QC  ex  cc  o (-1  u o o a a o o o o o o I- z 

3 3 3 3 3 3 3 3 3 3 3 3 U.  li.  U.  li.  Uu  U.  U.  UL  U.  UL  U.  li.  li.  ID  UJ 


OO  C3<M<MDI<M  >0<M<M<M<M.0 
D h3  i-^IM  iD 'T  ID  ^ <M  fO  D 
Dl<M  <M<M<MDI<M  <m4"4>4- 4- 4- 


O 


228 
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